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I .  INTRODUCTION 


A.  General  Introduction 


In  order  to  put  this  report  in  context  a  survey  of  some  penetration 
models  will  be  given  first.  For  this  purpose  it  is  convenient  to  divide 
these  models  into  elementary,  intermediate  and  advanced  types.  We  will 
not  attempt  a  complete  survey  of  all  such  models.  Instead,  we  will 
refer  to  or  quote  other  surveys  where  possible  and  content  ourselves 
with  an  analysis  of  some  of  the  more  popular  models  currently  in  use. 

A  recent  general  survey  of  penetration  models  has  been  given  by  Backman 
and  Goldsmith1 . 

B.  Elementary  Models 

A  principle  feature  of  elementary  models  is  that  they  are  very  easy 
to  use  and  require  little  calculation.  They  are  usually  constructed  by 
choosing  certain  experimental  variables  which  are  felt  to  be  important. 
These  are  incorporated  together  with  some  adjustable  constants  into  an 
arbitrary  functional  form  which  serves  the  ad  hoc  purpose  of  producing 
at  least  qualitatively  correct  behavior.  The  constants  are  then  adjusted 
to  give  a  reasonable  or  even  "best"  fit  to  the  data  which  is  available 
in  certain  ranges  of  the  allowed  values  of  the  experimental  variables. 

The  final  result  is  an  interpolation  formula  which  is  useful  within 
these  measured  ranges. 

One  reason  for  choosing  a  particular  functional  form  rather  than 
standard  multivariate  interpolation  formulas  such  as  those  of  Newton  or 
Lagrange  might  be  ease  of  calculation,  although  this  reason  has  less 
force  in  our  day  of  convenient  electronic  calculators.  A  better  reason 
might  be  given  if  the  functional  form  chosen  is  reminiscent  of  physical 
laws  so  that  it  can  be  helpful  in  mentally  organizing  a  mass  of  data. 

We  must  bear  in  mind,  however,  that  such  functional  forms  are  not  based 
on  physical  laws  and  are  purely  empirical.  A  danger  can  arise  if  the 
user  begins  to  place  too  much  confidence  in  such  formulas  and  uses  them 
for  extrapolation  beyond  the  measured  ranges  of  the  experimental  vari¬ 
ables.  At  best  the  predictions  of  such  an  extrapolation  are  highly 
tentative  and  are  explicitly  described  as  such.  At  worst  such  pre¬ 
dictions  can  be  absurd  without  the  user  even  realizing  this  fact. 

Numerous  examples  of  elementary  or  empirical  models  exist  and  many 
of  them  have  been  surveyed  by  Zook2,  We  will  select  one  for  further 
analysis  since  it  is  both  ambitious  and  popular.  This  model  consists 


1M.  E.  Backman  and  W.  Goldsmith,  "The  Mechanics  of  Penetration  of 
Projectiles  Into  Targets Int.  J.  Eng.  Sci.  16_,  1-99  (1978). 

2J.  Zook,  "An  Analytical  Model  of  Kinetic  Energy  Projectile /Fragment 
Penetration, "  BRL  ME  2797,  1977. 
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of  the  two  Thor  equations3.  Five  experimental  variables  were  chosen: 
striking  mass  (mp) ,  projectile  presented  area  (A),  striking  velocity, 
that  is,  speed  (VQ)  and  obliquity  (0O) ,  and  plate  target  thickness  (T) . 
Five  adjustable  constants  for  each  target-projectile  material  combina¬ 
tion  were  incorporated  together  with  these  variables  into  each  of  two 
arbitrary  functional  forms,  one  for  the  residual  speed  (V^) 

vi  =  v0  -  10*1  e™)*2  (ro/3  (sec  V*4  (V*5  (1) 

and  another  for  the  residual  mass  (n^) 

b1  b  b  b  b 

=  mQ  -  10  (TA)  (mQ)  (sec  0q)  (Vq)  (2) 

remaining  after  perforation.  Here  we  use  numerical  subscripts  instead 
of  letters  to  anticipate  multiplate  target  applications.  For  example, 

Vj  is  the  speed  remaining  after  perforation  of  target  plate  number  one, 
which  of  course  becomes  the  striking  speed  for  target  plate  number  two 
with  V2  remaining  and  so  forth.  Similarly  we  have  mj  and  0,.  These 
equations  were  fitted  to  a  particular  data  set  for  steel  cylinders 
impacting  ten  different  target  plate  materials  and  the  constants  were 
adjusted  to  give  a  least  squares  type  of  best  fit. 

In  the  Thor  equations  dimensions  are  preserved  by  suitable  choices 
for  the  power  of  ten  factors  which  govern  the  values  of  a1  and  b-, .  In 
Equation  (1)  the  fitting  procedure  gave  positive  values  for  a2  and  a4 
in  all  ten  cases  while  a^  was  always  negative  and  a^  was  positive  but 
less  than  unity  for  eight  target  materials  and  negative  for  two, 
magnesium  and  aluminum.  Since  a2  and  a 4  are  positive  these  parts  of 
Equation  (1)  are  reminiscent  of  physical  behavior.  The  factor  (TA)  is, 
roughly  speaking,  the  volume  of  the  plug  which  is  often  punched  out  of 
thin  target  plates  by  blunt  projectiles  and  a2  >  0  implies  that  the 
larger  this  plug  is  the  smaller  the  residual  speed  will  be.  This  is 
reasonable  since  more  energy  must  be  spent  punching  out  and  accelerating 
a  larger  plug.  One  might  think  that  an  association  of  T  with  sec  0O 
should  also  have  been  made.  A  factoraL  =  T  sec  0O  is  the  line-of-sight 
thickness  of  the  plate  and  (T  sec  0O)°^  with  a^  >  0  would  imply  that  a 
larger  L  due  either  to  larger  T  or  larger  0n  means  a  smaller  V1 ,  which 
is  also  reasonable.  Unfortunately,  there  is  no  easy  way  to  split  the 
a2 

factor  (T)  as  given  to  remind  us  of  both  plugs  and  line-of-sight  thick¬ 
ness. 

The  authors  of  the  Thor  report  cited  above  went  on  to  assume  that 
for  compact  projectiles  the  area  (A)  is  proportional  to  the  two-thirds 

power  of  the  striking  mass  and  absorbed  (A)  ^  into  the  mass  factor, 


^Project  Thor  Technical  Report  No.  47,  1961. 
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changing  a3  to  and  of  course  aj  to  aj*.  When  this  is  done,  T  is 
left  free  to  be  associated  with  sec  0Q  and  equality  of  a.2  and  a4  would 
remind  us  of  line-of-sight  thickness.  Usually  and  a4  differ  by  about 
25  percent  so  such  an  interpretation  is  not  impossible.  The  new  value 
given  for  a3*  is  not  a3  +  0.67  as  we  might  expect.  Usually  it  is  more 
like  aj/3. 


If  we  set  Vj  in  Equation  (1)  equal  to  zero  and  interpret  VQ  as  the 
perforation  limit  speed  V0£  we  have  for  the  compact  projectile 


(*o)  °3  ow*"*5  =  (10)O1  (T)~2  (sec  V~4  • 


(3) 


If  a_  were  (-  .5)  and  a5  were  zero,  then  squaring  Equation  (3)  and 
multiplying  by  one-half  would  give 


1  21  ^ai  ^a2  ^a4 

7  »„  Vt  ■  7  (10>  m  (sec  V 


(4) 


This  result  reminds  us  of  a  threshold  energy  for  perforation  proportional 
to  the  square  of  something  which  resembles  a  line-of-sight  thickness. 

For  magnesium,  aluminum,  titanium  and  homogeneous  steel  target  plates, 
a  fit  to  the  data  set  gave  a3*  closer  to  (-  .4)  than  (-  .5).  In  addition 
ar  was  usually  close  to  zero,  especially  for  steel.  However,  for  copper 
and  lead  a,*  was  closer  to  (-  .2)  while  a5  was  about  (.8).  This  gives 

(mQ  VoSL)°'2  on  the  left  side  of  Equation  (3)  so  that  by  taking  the  fifth 
power,  we  have 


t»o  V  -  <10> 


5a4 

(sec  6q) 


(5) 


which  is  suggestive  of  a  threshold  momentum.  When  a3  is  negative  as  for 
magnesium  and  aluminum  we  obtain  a  threshold  value  for  a  variable  mQ 
where  x  >  2  and  physical  reminders  are  difficult  to  find. 

If  we  use  V0£  from  Equation  (4)  in  Equation  (1)  for  a  compact  pro¬ 
jectile  with  a^*  =  -  ,5  and  a,-  =  0  we  obtain 


V  =  V  -  10 
1  o 


3-  3. 

(mo)"-5  (T)~2  (sec  0Q)  4  =  Vq  -  V 


o£ 


(6) 


Thus  the  Thor  prediction  for  the  most  common  target  plates  is  a  roughly 
linear  relation  between  Vj  and  VQ.  This  is  in  contrast  to  the  hyperbolic 

relation  V  «  (V  2  -  V  2)1//2  so  often  proposed  using  simple  physical 

J.  0  OX/ 
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considerations  and  presumably  confirmed  by  numerous  experiments.  If  we 
put  Equation  (5)  in  Equation  (1)  for  a  compact  projectile  with  a,*  =  .2 
and  a^  =  . 8  5 


V  =  V  -  10  1  (m  )'-2  (T)  2  (sec  6  )&4  (V  ) ‘8  =  V  ‘8  (V  '2  -  V  !2) 

1  o  o  v  o  o  o  o  o£  J 

(7) 

which  also  contrasts  with  the  hyperbolic  relation.  The  latter  of  course 
is  derived  for  zero  obliquity  and  constant  striking  mass.  The  Thor 
attempt  to  include  non- zero  obliquity  and  mass  erosion  in  a  simplistic 
formula  can  of  course  lead  to  such  contrasts. 

These  relations  are  noted  here  to  show  that  Equation  (1)  is  not  only 
easy  to  use  but  can  be  reminiscent  of  some  sort  of  physical  behavior  and 
so  aid  in  mentally  organizing  a  mass  of  data,  provided  we  take  note  of 
such  relations.  If  we  do  not,  the  formula  can  be  as  confusing  as  the 
data  itself. 

In  Equation  (2)  bj  and  b^  were  always  found  to  be  positive.  Positive 
b2  implies  that  a  thicker  target  plate  or  larger  plug  will  cause  greater 
mass  loss  which  is  reasonable.  Positive  be;  implies  that  higher  velocity 
projectiles  will  suffer  more  erosion  which  is  not  in  accord  with  experi¬ 
ence  for  unbroken  projectiles  as  we  will  see  below.  However,  if  one 
looks  only  at  the  largest  piece  of  a  broken  projectile,  then  positive 
b$  can  describe  experiment.  In  most  cases  b^*  is  positive  and  not  much 
less  than  unity  which  suggests  the  not  unreasonable  form 


mj  =  mQ  [1  -  f  (Vo,  9o,  T)].  (8) 

In  the  single  case  of  titanium  b^*  =  -  .024.  This  is  inconsistent  with 

with  the  other  target  materials  used  and  could  lead  one  to  suspect  the 
data  and/or  the  fitting  procedure.  The  thought  that  titanium  is  somehow 
different  and  that  the  Thor  results  are  revealing  a  new  phenomenon  should 
be  resisted  since  we  are  dealing  with  an  empirical  formula  not  based  on 
any  physical  law. 

For  eight  target  materials  b^  was  found  to  be  positive  which  is 
reasonable  since  greater  mass  loss  is  expected  for  a  greater  line-of- 
sight  thickness.  For  two  target  materials  b^  was  found  to  be  negative 
with  b^  =  -.172  for  magnesium  and  b^  =  -  .361  for  aluminum.  This  leads 
to  the  surprising  prediction  that  a  greater  line-of-sight  thickness  due 

to  an  increase  in  obliquity  means  less  mass  loss.  In  addition,  it  is 

inconsistent  with  the  other  target  materials  used,  although  interestingly 
enough,  the  value  found  for  a^  for  these  same  two  materials  was  also 
inconsistent  with  the  other  materials.  Again  we  might  suspect  the  data 
and/or  the  fit.  Indeed,  at  least  for  aluminum  two  obvious  blunders  can 
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be  found  which  might  help  to  account  for  the  negative  value  of  b^.  In 
experiments  256  and  257  the  residual  mass  was  reported  to  be  three  or 
four  times  the  striking  mass.  If  a  simple  typographical  blunder  is  not 
responsible  then  possibly  a  failure  to  pick  out  the  remains  of  the  steel 
projectile  could  be  the  cause.  For  example,  if  a  large  piece  from  the 
aluminum  target  were  chosen  instead  of  the  projectile  or  if  a  piece  of 
the  target  adhering  to  the  projectile  were  included  in  the  weighing 
procedure,  perhaps  an  explanation  could  be  found,  This  type  of  mistake 
is  commonly  called  a  blunder  to  distinguish  it  from  an  experimental 
error  which  can  be  reduced  but  not  removed  entirely.  Again  we  should 
resist  the  temptation  of  thinking  that  new  properties  of  magnesium  and 
aluminum  might  be  revealed  and  await  discovery,  reminding  ourselves  that 
we  are  dealing  with  a  purely  empirical  formula.  An  argument  might  be 
made  that  target  pieces  adhering  to  the  projectile  should  be  recorded  as 
the  residual  mass  since  this  combination  can  inflict  more  damage  behind 
the  target.  If  we  accept  this  argument,  then  to  be  consistent  we  should 
always  record  the  masses  of  target  pieces.  This  is  desirable,  but  was 
never  done  in  the  Thor  collection. 

If  we  try  to  verify  an  increase  in  residual  mass  with  increasing 
obliquity  (negative  b4)  in  the  Thor  data  sets  for  magnesium  and  aluminum 
we  quickly  find  that  the  data  is  too  sparse  and  inaccurate  to  make  a 
judgment.  If  we  hold  three  factors  constant  for  magnesium,  namely, 

T  =  .75  inches,  mQ  =  240  grains  and  VQ  =  5,400  ft/sec,  we  find  mx  equal 
to  195.7,  165.4  and  230.7  grains  respectively  for  0Q  equal  to  0°,  60° 
and  70°.  It  is  difficult  to  believe  that  a  10°  increase  in  obliquity 
could  lead  to  a  forty  percent  increase  in  residual  mass,  and  similar 
behavior  is  not  found  elsewhere.  Such  data  is  usually  labeled  anomalous 
and  calls  for  a  repetition  of  the  experiment  before  inclusion  into  a 
formula  fit.  If  w*e  keep  T  =  1  inch,  mn  *  120  grains  and  VQ  =  5,000 
ft/sec,  we  find  m-j  =  104.5  and  119.0  (average  111.75)  grains  at  6p  =  0  , 
m,  =  103.7  grains* at  45°  and  105.2  and  106.5  (average  105.85)  grams  at 
60°.  Here  we  have  more  than  one  experiment  in  two  cases  at  least  and 
the  trend  is  not  inconsistent  with  decreasing  or  even  constant  mj  as  a 
function  of  60  in  view  of  the  experimental  errors  involved.  Other  data 
with  T,  m0  and  VQ  constant  was  not  obtained.  If  we  allow  VQ  to  vary  a 
bit  we  find  for  T  =  1  inch  and  mQ  =  240  grains  that  for  0Q  =  0°  and 
Vn  =  4,649  ft/sec,  m-,  =  214.7  grains;  for  0n  =  45°  and  VQ  =  4,763  ft/sec, 
m^  -  202,2  grains;  and  for  0Q  =  60°  and  VQ  =  5,000  ft/sec,  m^  =  185.0 
grains.  In  addition  we  find  for  T  =  2  inches  and  m0  =  240  grains  that 
for  0O  =  0°  and  V  =  5,000  ft/sec,  mx  =  213.3  and  219.0  (average  216.15) 
grains  while  for  0Q  =  45°  and  VQ  -  5,500  ft/sec,  m^  =  174.7  and  190.7 
(average  182.7)  grains.  Admittedly  this  information  is  complicated  by 
the  fact  that  increasing  V  should  also  decrease  mj ,  but  there  is  no 
evidence  to  support  an  increase  of  m^  with  increasing  0Q. 

The  residual  mass  data  for  aluminum  is  even  more  sparse  than  for 
magnesium  and,  as  we  have  noted,  contains  some  blunders.  About  the  only 
example  close  to  what  we  want  is  for  T  =  0.5  inches  and  mQ  =  240  grains 
where  we  find  that  nij  =  229.8  grains  for  0Q  =  45°  and  VQ  =  2,662  ft/sec 
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while  m,  =  222.1  grains  for  0Q  =  60°  and  VQ  =  2,638  ft/sec.  Here  the 
slight  decrease  in  V  should  tend  to  increase  m^,  yet  in  spite  of  this 
m^  decreases  slightly  as  0Q  increases.  Of  course  with  experimental 
errors  in  the  order  of  10  to  20  percent  for  mi  and  at  least  10  percent 
for  V0  it  is  difficult  to  say  anything  about  m^  as  a  function  of  0  from 
this  information.  u 

What  are  we  to  make  of  a  negative  value  for  or  a^  in  the  cases  of 
magnesium  and  aluminum?  The  fitting  error  was  no  worse~for  these 
materials  than  for  the  other  eight.  An  explanation  can  be  found  for  this 
anomaly  in  the  arbitrary  functional  form  of  the  model  in  which  a  simple 
product  of  variables  is  strung  together  and  given  adjustable  exponents 
fitted  to  a  data  set  which  is  unevenly  weighted  over  the  measured  ranges 
and  which  contains  large  experimental  errors  and  occasional  blunders. 

If  a  different  mathematical  form  had  been  chosen  to  fit  the  same  data 
set,  the  result  could  well  be  a  decrease  in  residual  mass  with  increasing 
obliquity.  The  same  result  might  be  obtained  if  the  same  mathematical 
form  were  fitted  to  a  different  (more  accurate,  more  evenly  weighted) 
data  set.  These  observations  serve  to  illustrate  the  recognized  but 
sometimes  ignored  fact  that  purely  empirical  models  add  nothing  to  our 
fundamental  understanding  and  can  give  erroneous  predictions  if  they  are 
used  for  extrapolation. 

For  example,  a  fit  of  Equation  (1)  to  another  data  set  might  give  a 
negative  value  for  a^,  implying  that  would  increase  with  increasing 
obliquity,  contrary  to  experience.  Or,  it  might  give  a  value  of  a^ 
greater  than  unity,  implying  a  decrease  in  Vj  as  VQ  increases  with  the 
absurd  prediction  of  embedment  or  ricochet  rather  than  perforation  at 
sufficiently  high  VQ .  This  kind  of  possibility  is  not  idle  speculation 
as  is  clear  from  a  recently  published  extension  of  the  Thor  procedure  to 
tungsten  fragments  penetrating  steel  and  aluminum  targets14.  There  it 
was  found  that  a  fit  of  Equation  (2)  to  the  data  for  steel  targets  gave 
a  value  of  b^  greater  than  unity,  implying  a  decrease  in  residual  mass 
with  increasing  striking  mass  with  complete  disappearance  of  m^  for 
sufficiently  large  m  .  Fortunately,  the  authors  point  out  this  particular 
danger.  Hopefully  users  of  their  formula  will  bear  this  in  mind  and 
explicitly  mention  this  behavior  if  they  are  bold  enough  to  use  this 
formula  for  extrapolation.  It  is  interesting  to  note  in  this  report  that 
b„  is  positive  for  aluminum  as  well  as  for  steel  in  contrast  to  the 
negative  value  found  in  the  Thor  case.  This  tends  to  confirm  the  obser¬ 
vations  made  above  on  the  nature  of  purely  empirical  formula  fits  to 
particular  data  sets. 

Thor  Report  No.  47  does  not  make  it  clear  exactly  how  the  residual 
speed  was  measured.  It  merely  states  that  most  of  the  data  was  provided 
by  BRL  but  that  other  sources  were  also  included  (no  reference  given). 

4£.  L.  Holloway 3  M,  B.  Danish  and  J.  A .  Matts*  " Penetration  Relations 
for  Tungsten  Alloy  Fragments  versus  Selected  Target  Materials 
ARBRL-TR-02087 3  1978.  (AD  #0016203) 
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Other  Thor  Reports  before  and  after  No.  47  are  also  of  little  help.  It 
is  very  likely  that  velocity  screens  were  used,  although  depth  of  pene¬ 
tration  into  a  catcher  material  is  also  a  possibility.  In  making  such 
measurements  for  oblique  impacts  it  is  common  practice  to  cock  the  target 
plate  at  the  desired  angle  but  to  leave  the  velocity  screens  or  witness 
plates  normal  to  the  axis  of  the  gun  which  is  used  to  launch  the  pro¬ 
jectiles.  This  is  done  not  only  because  of  space  limitations  or  con¬ 
venience,  but  also  because  the  deviation  of  the  projectile  from  its 
original  flight  path  is  not  known  ahead  of  time  and  would  require  addi¬ 
tional  experiments  and  equipment  to  determine  if  it  could  not  be  calcu¬ 
lated.  This  practice  can  lead  to  errors  in  the  reported  magnitude  of 
the  residual  velocity  since  only  one  component  of  this  velocity  is  being 
measured.  In  cases  where  the  obliquity  is  not  large  or  the  impact  speed 
is  not  greatly  different  from  the  exit  speed,  the  neglect  of  an  exit 
velocity  component  perpendicular  to  the  line  of  fire  is  not  serious 
since  the  deviation  is  small.  However,  for  large  obliquities  and  marginal 
perforations  this  neglect  can  lead  to  serious  errors.  This  practice  is 
illustrated  in  Figure  3  (p  40)  of  Thor  Report  No.  47  where  the  assumption 
of  normal  impact  on  a  second  target  is  explicitly  mentioned  (last  line  of 
p.  38). 

C.  Intermediate  Models 


Intermediate  models  are  harder  to  use  and  require  more  calculation, 
but  not  enough  to  preclude  a  description  of  the  main  features  of  practical 
problems  using  a  reasonable  expenditure  of  time  and  money.  They  are 
based  on  some  type  of  physical  law  which  is  assumed  to  govern  the  motion 
of  the  projectile  through  the  target  and  usually  they  treat  the  projectile 
as  a  point  mass  and  contain  other  simplifying  assumptions  as  well  as 
adjustable  constants.  They  come  under  the  general  heading  of  particle 
physics  and  as  such  can  be  called  semi -empirical  or  semi-theoretical. 
Although  they  require  more  effort  to  use  than  elementary  models,  they 
can  be  used  with  some  confidence  in  tentative  extrapolations  beyond  the 
measured  ranges  of  the  variables.  The  degree  of  confidence  we  assign  to 
such  extrapolations  will  depend  on  the  kind  of  physical  laws  that  are 
assumed  to  apply  and  the  extent  of  the  agreement  between  model  predic¬ 
tions  and  known  facts.  Ideally  such  models  will  predict  at  least 
qualitatively  correct  behavior  even  for  the  extreme  values  permitted  for 
the  variables. 

A  simple  example  of  this  type  of  model  is  contained  in  a  1963  paper 
by  Recht  and  Ipson5.  Using  a  result  of  Jameson  and  Williams6  that 
projectile  and  target  plug  residual  speeds  are  practically  the  same, 
these  authors  wrote  an  energy  balance  as  follows: 


5R .  F.  Recht  and  T.  W.  Ipson ,  " Ballistic  Perforation  dynamics , "  Trans¬ 
actions  of  the  ASME ,  p.  384,  1963. 

6fl.  L.  Jameson  and  J.  S.  Williams,  " Velocity  Losses  of  Cylindrical  Steel 
Projectiles  Perforating  Mild  Steel  Plates,"  BRL  R  1019,  1957. 

(AD  # 142447 ) 
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where  mn  is  the  plug  mass,  W  is  the  work  done  in  shearing  out  this  plug 
and  E  i§  the  energy  lost  in  other  ways.  The  authors  assume  that  E  can 
be  equated  to  the  energy  transfer  which  might  occur  in  a  free  collision 
between  mQ  and  nu  as  if  the  rest  of  the  target  were  not  present.  Since 
momentum  conservation  in  such  a  free  collision  is 


m  V  =(m  +  m)V1_ 

o  o  o  p  If 


(10) 


with  the  subscript  f  standing  for  free  collision,  we  have 


E  =  Ef  =  7  "o  (Vo2  •  Vlf>  ■  7  mo  Vo2 


,  /  m  V 

1  I  o  o 
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(11a) 


or 


(lib) 


We  have  put  Equation  (11a)  into  the  equivalent  form  (lib)  in  order  to 
point  out  that  the  authors  implicitly  assumed  in  their  paper  that  the 
factor  (2  m^  +  m^)/ (mQ  +  m^)  can  be  set  equal  to  unity  when  they  wrote 
their  Equation  (2).  Of  course  this  requires  mQ  to  be  zero  to  be  strictly 
correct  or  at  least  <<  nip  to  be  approximately  true.  In  case  mQ  =  nip 
(a  not  unreal  situation  in  some  cases),  this  factor  is  equal  to  1.5 
(not  too  far  from  unity)  and  it  approaches  2  if  mQ  >>  mp.  We  will  not 
make  this  approximation  here  and  will  note  the  consequences. 

At  the  perforation  limit,  V1  =  0  and  VQ  =  VQ^.  In  this  case  a 
combination  of  Equation  (11a)  and  Equation  (9)  gives  us 


o 


(12) 


which  is  the  authors'  Equation  (4)  except  for  the  second  power  of  the 
mass  ratio  instead  of  the  first  power  given  by  the  authors.  This  of 
course  does  not  affect  their  argument  based  on  high-speed  machining 
studies  that  W  can  be  taken  to  be  independent  of  VQ  >  V0£  to  a  good 
approximation,  since  the  dynamic  shear  stress  is  also  approximately 
constant. 
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Now  if  we  put  Equations  (11a)  and  (12)  in  Equation  (9)  and  solve 
for  Vj,  we  obtain 


which  is  the  authors'  Equation  (5)  except  for  the  3/2  power  of  the  mass 
ratio  instead  of  their  first  power.  This  is  interesting  because 
Equation  (13)  above  is  exactly  the  authors'  Equation  (14)  which  they 
obtained  by  postulating  an  extra  factor  of 

P-T 

\°  7 

to  account  for  perforation  of  thick  plates.  If  mp  <<  mQ  in  Equation 
(13),  then  we  have  the  simple  hyperbolic  relation  mentioned  in  the 
previous  subsection. 

Another  tradition  exists  in  this  category  of  models  which  analyze 
the  particle  dynamics  of  penetration.  It  consists  of  assuming  a  force 
field  in  which  the  projectile  moves  while  it  is  in  the  target,  solving 
Newton's  second  law,  and  comparing  with  experiment.  Newton  himself 
apparently  followed  this  procedure  in  his  theory  of  gravitation.  A 
learning  exercise  used  in  some  physics  texts  asks  the  student  to  assume 
various  forms  for  a  gravitational  force  field,  a  simple  one  being 

a  a?  at 

C  (n^)  1  (m2)  /(r)  ,  where  the  notation  is  obvious.  Then  the  student 

is  required  to  fit  the  adjustable  constants  C,  a^,  a^  and  a  3  to  Kepler's 
data  and  to  give  values  for  these  constants  as  well  as  their  average 
error,  for  example,  a$  =  2  ±  .002.  Explanations  for  this  error  in 
addition  to  experimental  uncertainties  are  then  elicited,  for  example, 
perturbations  by  the  other  planets,  and  so  forth.  This  method  is  so 
powerful  that  Newton's  contemporaries,  Robins  and  Euler,  soon  applied  it 
to  projectile  penetration  of  dense  media7.  They  assumed  a  constant  force 
field  opposing  the  motion  and  for  constant  projectile  mass  and  zero 
obliquity  their  equation 


m  dv/dt  =  -  a  (14) 

o 

is  easily  solved.  A  generalization  of  Equation  (14)  to  include  a 
speed-dependent  force  field  is 


7P/.  Johnson ,  Impact  Strength  of  Materials ,  N.Y.:  Crane ,  Russak3  1972 , 
cc.  9-10. 
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(15) 


2 

dv/dt  =-  (a+bv+cv) 

where  a,  b  and  c  are  constants.  This  equation  with  b  =  0  is  usually 
attributed  to  Poncelet  (c.1830),  while  the  same  equation  with  a  =  0  is 
generally  associated  with  Resal  (c.  1895).  Of  course  Equation  (15)  is 
also  easily  solved8. 


Robertson's  1943  report9  is  an  example  of  a  paper  which  solves 
Equation  (15)  with  b  =  0.  This  author  factors  the  constant  a  into  Aa' 
where  A  is  the  "area  of  impression"  of  the  projectile  and  a'  is  called 
a  "shatter  coefficient".  He  also  allows  A  to  be  a  function  of  the  depth 
of  penetration,  but  he  proposes  no  explicit  form  for  A  except  in  the 
case  of  a  cylinder  when  A  is  taken  to  be  a  constant  ttR^  for  a  non-deform¬ 
ing  cylinder  of  radius  R.  He  also  lets  C  =  A(l/2  y  pt)  where  pt  is  the 
density  of  the  target  material  and  y  is  a  dimensionless  constant  which 
he  calls  the  "inertial  coefficient".  He  then  uses  the  relation  dv/dt  = 
v  dv/dz  where  z  is  the  depth  of  penetration  to  integrate  Equation  (15) 
once  using  limits  VQ  and  Vj  for  v  and  0  and  T  for  z,  obtaining 


-cT/m  1/2 


2  2  cT/m 

where  Vq£  =  (a/c) (e  °  -  1). 

is  the  plug  mass  and  1/2  pt  (a/c)  =  (a'/y).  The  ratio  (a/c)  has  the 
dimensions  of  a  velocity  squared  while  (a'/y)  has  the  dimensions  of  a 
stress.  Here  the  same  type  of  hyperbolic  form  as  above  has  been  obtained 
directly  from  Newton's  law  with  an  analytical  rather  than  an  empirical 
form  for  the  limit  speed.  Of  course  this  form  contains  the  adjustable 
constants  a'  and  y. 

Lambert  and  coworkers10  suggest  a  generalization  of  Equations  (13) 
and  (16)  to  include  obliquity  and  rod  length,  namely, 


(V, 


(16) 


Here  2  cT  =  y 


Pt  A  T 


=  Y 


m  where  m 

n  n 


®W.  Gr'dbner  and  N.  Hofrez ter,  "Integraltafel.  Erster  Teil.  Unbestimmte 
Integrate  ". 

9 H .  P.  Robertson,  "The  Mechanics  of  Armor  Perforation,"  NDRC ,  Armor  and 
Ordnance  Report  No.  A-227  (OSRD  No.  2043),  1943. 

10e/.  P.  Lambert  and  G.  H.  Jonas,  "Towards  Standardization  in  Terminal 
Ballistics  Testing:  Velocity  Representation ,"  BRL  Report  1852,  1976; 

J.  P.  Lambert  and  B.  E.  Ringers,  "Standardization  of  Terminal  Ballistics 
Testing,  Data  Storage  and  Retrieval, "  ARBRL-TR-02066 ,  1978 ;  and  J .  P. 
Lambert ,  "The  Terminal  Ballistics  of  Certain  65  Gram  Long  Rod  Penetra— 
tors  Impacting  Steel  Armor  Plates,"  ARBRL-TR-02072,  1978. 
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(17) 


The  constant  p  is  2  plus  an  empirical  factor  which  depends  on  rod  diameter, 
obliquity  and  target  thickness.  VQ£  is  an  empirical  factor  related  to 
these  quantities  as  well  as  to  the  rod  mass  and  length.  The  mass  m' 
bears  a  strong  resemblance  to  the  plug  mass,  nip.  If  the  exponential 
factor  in  Equation  (16)  is  expanded  in  powers  bf  (nip/m0) ,  neglecting 
powers  higher  than  the  first,  it  becomes  1  -  (y/2)  Xnip/m0) .  If  the 
mass  ratio  in  Equation  (17)  is  similarly  expanded  we  nave  1  -  1/3  (m'/m0) 

where  m'  =  pt  (ttR2)  T  (sec  90)3^4  with  R  the  rod  radius.  If  sec  0Q  had 
been  chosen  to  the  first  power  we  would  have  a  line-of-sight  thickness 
and  m'  roughly  equal  to  hl.  A  choice  of  y  =  2/3  would  complete  the  first 
order  equivalence  of  m'  and  mp.  A  similar  expansion  of  the  mass  ratio 
in  Equation  (13)  gives  1  -  3/2  (mp/mQ). 

Zook2  has  continued  this  tradition  by  applying  the  solution  of  Equa¬ 
tion  (15)  to  that  portion  of  the  Thor  data  which  was  obtained  for  zero 
obliquity.  This  was  a  wise  restriction  because,  as  we  shall  see  below  in 
the  body  of  the  present  paper.  Equation  (15)  is  only  applicable  to  the 
case  of  zero  obliquity.  Like  Robertson,  Zook  factors  the  constant  a  in 
Equation  (15)  and  writes  a  =  Cj  A  Ht  where  A  is  the  cross-sectional  area 
of  the  projectile  and  H^.  is  the  Brinell  hardness  of  the  target  material. 
Similarly,  he  writes  c  =  Cj  A  so  that  his  is  half  of  Robertson's 
"inertial  coefficient",  y.  Unlike  Robertson  he  retains  the  constant  b 

and  proposes  the  form  b  =  C2  A  vlTJT- pT.  However,  he  made  no  attempt  to 

justify  this  geometric  mean  type  of  form,  nor  to  justify  the  use  of  a 
hardness  coefficient  of  any  type  much  less  the  Brinell  hardness  which  is, 
of  course,  only  one  of  the  many  such  parameters  which  appear  in  the 
literature.  Consequently,  a,  b  and  c  remain  adjustable  or  empirical 
constants  in  an  analytical  theory. 

Backman  and  Finnegan11  modified  Equation  (15)  to  be 


2  2 

mQ  dv/dt  =  -  A  (aj  +  Cj  cos  ip  v  )  (18) 

where  \p  is  the  angle  between  the  projectile's  current  line-of-f light 
(tangent to  the  trajectory)  and  the  current  normal  to  the  spherical 
surface  which  is  in  contact  with  the  deforming  target.  In  order  to 


11M.  E.  Backman  and  Finnegan ,  "Dynamics  of  the  Oblique  Impact  and 
Ricochet  of  Ron  Deforming  Spheres  Against  Thin  Plates , "  NWC  TP  5844, 
1976. 
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determine  the  projectile  motion  for  oblique  impacts  they  also  postulate 
the  following  rule.  The  force  opposing  the  motion  is  directed  through 
the  center  of  the  sphere  and  lies  antiparallel  to  the  current  contact 
surface  normal  mentioned  above.  Given  this  postulate,  one  can  then 
resolve  this  force  into  a  resistive  component  antiparallel  to  the  current 
line  of  flight  and  a  lift  component  perpendicular  to  this  line.  Torques 
tending  to  rotate  the  sphere  are  neglected  and  the  change  in  speed,  Av, 
is  calculated  for  time  step.  At,  from  the  resistive  force  component, 
while  the  change  in  direction  is  calculated  from  the  lift  component. 

Next,  the  change  of  position  during  the  following  time  interval  is  cal¬ 
culated  using  the  values  of  the  speed  and  curvature.  Finally,  a  critical 
penetration  depth,  P*,  for  plate  failure  is  assumed,  beyond  which  depth 
part  of  the  resistive  force  component  is  assumed  to  vanish  over  a  portion 
of  the  contact  surface  specified  by  an  assumed  angle  a*.  Thus,  the  model 
contains  four  adjustable  constants,  aj,  c^,  P*  and  a*  which  are  determined 
from  data  obtained  at  zero  obliquity.  The  constants  a^  and  c^  were  found 
by  a  best  fit  to  Equation  (13)  above  using  the  Recht-Ipson  power  of  unity. 
The  constants  P*  and  a*  were  chosen  to  match  observation  of  crater  for¬ 
mation  at  zero  obliquity.  Additional  postulates  about  a  kind  of  virtual 
origin  for  the  sphere  at  the  beginning  of  its  motion  also  enabled  the 
authors  to  calculate  crater  shapes  for  ricochet  and  embedment  from  the 
projection  of  one  side  of  the  tube  which  is  swept  out  by  the  sphere  along 
its  trajectory.  Since  A  is  a  point  at  the  moment  of  impact  with  zero 
force  pointing  along  the  normal  to  the  target  face,  an  initial  depth  of 
penetration  must  be  assumed  in  order  to  begin  the  calculation.  In 
addition,  the  outcome  will  depend  somewhat  on  the  size  of  the  time  steps 
which  are  chosen  and  very  definitely  on  the  spherical  cap  which  is 
taken  to  be  the  contact  surface  at  any  instant. 

An  advantage  of  this  model  is  its  ability  to  predict  exit  angles  as 
well  as  exit  speeds  for  oblique  impacts.  In  addition  to  perforation 
phenomena,  ricochet  and  embedment  phenomena  are  also  treated,  something 
which  is  rarely  if  ever  done  by  other  models  in  this  category.  An 
unsatisfactory  feature  of  this  model  is  the  somewhat  arbitrary  nature  of 
the  rules  which  have  been  selected  for  calculating  force  components  and 
crater  shapes.  Still,  it  appears  to  be  the  most  advanced  example  of  its 
kind  in  the  particle  mechanics  tradition  of  penetrator  models.  In  using 
a  numerical  time  and  space  step  method  for  calculating  the  projectile 
trajectory,  this  model  resembles  some  of  the  computer  code  techniques 
usually  reserved  for  solving  the  field  equations  which  are  referred  to 
in  the  next  subsection. 

D.  Advanced  Models 


Even  Newton's  law  of  gravitation  has  been  superceded  by  tensor  field 
equations  which  predict  the  existence  of  gravitational  waves  and  are 
currently  being  used  in  "black  hole"  theories  designed  to  explain  certain 
astrophysical  observations12.  Similarly,  a  large  number  of  investigators 

12L,  L._  Smarr  and  W.  H.  Press ,  "Our  Elastic  Spacetime:  Black  Holes  and 
Gravitational  Waves,"  American  Scientist  66J1)  Jan-Feb  1978,  p.  72. 
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have  abandoned  particle  physics  descriptions  and  turned  to  solutions  of 
the  field  equations  of  continuum  mechanics  in  an  effort  to  model  the 
penetration  process.  A  summary  of  such  efforts  can  be  found  in  a  paper 
by  Jonas  and  Zukas13  who  also  review  a  number  of  particle  models.  We 
will  content  ourselves  here  with  a  few  quotations  from  their  paper: 

"Computer  codes  are  invaluable  for  obtaining  a  qualitative 
picture  of  penetrator  and  target  deformation.  They  provide 
details  not  normally  obtainable  from  ballistic  experiments. 

Yet,  it  is  not  unfair  to  state  that,  on  the  whole,  they  have 
not  improved  our  understanding  of  penetration  phenomenology, 
except  possibly  in  the  hypervelocity  regime." 

With  few  exceptions 

"....the  present  state-of-the-art  in  computations  is  best 
summarized  by  a  paraphrase  of  Richard's  Law,  to  wit:  'One 
good  guess  is  worth  a  thousand  computer  runs'." 

"Aside  from  the  above  considerations,  codes  in  their  present 
form  cannot  be  relied  upon  for  quantitative  data  (except  by 
accident)  in  the  ordnance  velocity  regime." 

"It  is  necessary  to  add  here  that  despite  these  limitations, 
code  results  for  deformation  fields  often  bear  reasonable 
resemblance  to  those  found  experimentally." 

"Refinements  in  theory  will  be  to  little  avail  however  until 
characterization  of  materials  at  the  strain  rates  appropriate 
to  ballistic  impact  conditions  (up  to  10°  s"x)  is  achieved, 
for  errors  in  input  will  outweigh  any  gains  in  modeling." 

"The  last  and  most  difficult  hurdle  to  overcome  will  be 
determination  of  reliable  and  computationally  suitable  models 
for  the  onset  of  fracture  and  the  characterization  of  failed 
material. " 

"The  use  of  such  codes  is  neither  straightforward  nor  in¬ 
expensive.  A  typical  code  will  output  about  10'  words  of 
information  and  cost  upwards  of  $1000  per  run.  Of  necessity 
much  of  the  output  is  presented  in  plot  form.  Considerable 
experience  is  required  to  run  the  codes  and  frequently 
manual  intervention  is  called  for.  In  no  way  can  present  codes 
be  treated  as  'black  boxes'." 


13J.  H.  Jonas  and  J.  A.  Zukas ,  "Armor  Penetration:  Theory  and  Experiment," 
USA  BEL,  invited  paper  at  the  14th  Annual  Meeting  of  the  Society  of 
Engineering  Science ,  Lehigh  University,  PA,  1977,  Also  published  as 
(' Mechanics  of  Penetration:  Analysis  and  Experiment,"  Int.  J.  Eng.  Sci. 
16,  879 ,  1978. 
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Although  Newton's  gravitational  law  has  been  superceded  for  some 
purposes  it  is  still  extremely  useful  for  most  terrestial  and  solar 
system  phenomena.  Similarly,  while  continuum  models  may  be  useful  for 
some  purposes  and  may  eventually  become  routine,  our  present  requirements 
for  further  fundamental  understanding  as  well  as  our  needs  for  relatively 
simple,  inexpensive  methods  of  making  quantitative  predictions  for  at 
least  the  main  features  of  practical  problems  prompt  us  to  investigate 
further  the  possibility  of  constructing  a  better  particle  model.  The 
purpose  of  the  present  paper  is  to  construct  such  a  model  so  we  can 
describe  the  penetration  process  in  a  simple  yet  rational  manner.  This 
should  have  the  double  effect  of  increasing  our  understanding  and  provid¬ 
ing  a  useful  tool  for  making  predictions. 
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II.  LIMITATIONS  OF  A  SPEED - DEPENDENT  FORCE  FIELD 


Even  the  most  general  force  field  usually  considered  (Equation  (15) 
above)  is  severely  limited.  We  can  note  this  by  writing  it  in  vector  form 
as  follows: 

?=  -  (a  +  b  v  +  c  vz)  v  =  -  (i+b  +  cv)^  (19) 

where  v  =  v/v  is  a  unit  vector  in  the  current  direction  of  the  velocity 
v  which  has  magnitude  v.  If  r  is  a  position  vector  from  the  coordinate 


origin  to  the  center  of  mass  of  the  projectile,  then  in  Cartesian 
coordinates  r  =  x  i  +  y  j  +  z  k  where  i,  j  and  fc  are  constant  unit  vec- 

j  r  ->■  — v  ,  c  •  c  .  c  .  ,  f .2  .2  .2.1/2 

tors  and  of  course  v=r=xi+yj  +zk  with  v  =  (x  +  y  +  z  ) 

Here  the  Newtonian  dot  convention  for  time  derivative  is  being  used.  In 

most  discussions  of  this  force  field  in  the  past,  the  fact  that  ?  is 


being  taken  antiparallel  to  v  is  mentioned  at  least  implicitly  by  the 
statement  that  only  zero  (or  normal)  obliquity  is  being  considered. 
Attempts  are  then  sometimes  made  to  "patch  up"  the  result  in  order  to 
make  application  to  impacts  at  non-zero  obliquity.  These  attempts  are, 

of  course,  in  vain  since  a  direct  consequence  of  ?  being  anti  parallel 


to  ^  at  the  beginning  of  the  motion  is  that  the  motion  is  necessarily 
rectilinear.  This  is  immediately  obvious  because,  if  there  is  no  force 

component  perpendicular  to  v  =  vQ  at  time  zero,  there  never  will  be  such 
a  force  component  and  the  projectile  must  continue  to  move  along  its 
original  straight-line  path.  In  short,  curvilinear  motion  is  impossible 
in  such  a  force  field.  Since  curvilinear  motion,  including  extreme 
cases  such  as  ricochet,  is  observed  experimentally,  then  a  different 
force  field  must  be  used  if  we  hope  to  describe  the  projectile  motion. 


This  observation  about  rectilinear  motion  can  also  be  demonstrated 
by  solving  Newton's  second  law  for  a  constant  mass  particle  with  Equation 
(19)  as  the  force.  We  may  write 


[m  -4-+(-+b  +  cv)]v=0  (20) 

1  o  dt  Vv  J 

where  for  convenience  we  have  factored  the  equation  into  an  operator 

(in  brackets)  and  an  operand,  If  we  take  the  scalar  or  dot  product 

of  this  equation  with  v,  using  ^  .  v  =  v,  we  obtain  the  usual  separable 
scalar  equation 


m 


dv 
o  dt 


(a  +  bv  +  cv)  =  0 


(21) 
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As  noted  above  in  the  introduction,  the  solution  can  be  looked  up  in 
standard  integral  tables.  If  our  interest  is  in  obtaining  v  =  v(s)  where 
s  is  the  displacement  or  distance  traveled  along  the  trajectory,  we  can 
use  the  relation  dv/dt  =  v  dv/ds,  since  v  =  ds/dt.  As  we  have  seen, 
Robertson  has  done  this  for  the  special  case  b  =  0,  while  Zook  has  done 
it  for  the  general  case.  Our  interest  here  is  in  obtaining  v  =  v(t) 
instead.  If  we  define  the  discriminant 

q  =  4  ac  -  b2  (22) 

then  for  q  <  0, 

V  =  IE  (h+  '  A  h"  e_Bt )/ Cl  -  A  e~Bt)  (23) 

+  r — +  _ 

where  h  =  -  b  ±  /-q,  A  =  (2  c  vq  -  h  ) / (2  c  vq  -  h  )  and  B  =  /Iq/m  , 
using  the  initial  conditions  v  =  v  ,  s  =  0  at  t  =  0,  with 

s  -  (f)  £n  C^A-V)  +  (^)  t  .  (24) 

Alternative  solutions  for  this  case  which  involve  hyperbolic  tangents 
will  not  be  considered  here  because  of  the  physical  requirement  that  v 
decrease  with  increasing  time. 

2  2 

For  q=0,  a+bv+cv  =  c  (V  +  v)  where  V 


V  =  [(V  +  vo)_1  +  (^-)  t]'1  - 


s  =  (— )  £n  [1  +  (^-)(V  +  v  )  t] 

o 

For  q  >  0, 

v  =  [/q  tan  (C  -  Dt)  -  b]/(2c)  (27) 


V  (25) 

-  Vt.  (26) 


s  =  [2  mQ  £n  (cos  (C  -  Dt)/cos  (C)}  -  b  t]/(2c)  (28) 

where  C  =  arc  tan  [(2c  vq  +  b)/*^]  and  D  =  /q/ (2  m^) .  Our  main  interest 
here  is  to  note  from  Equations  (23),  (25)  and  (27)  that  v  =  v(t)  can  be 
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found  explicitly. 

If  we  write  Equation  (20)  in  Cartesian  component  form  we  have,  for 
example, 

["o3F*  (^b  *  cv)1  *'  °  C29> 

r  sur.zss.s-s  r  wsr-s-"1;; 

and  we  obtain,  using  Equation  (15)  or  (21J , 


dx  _  d£  _  _ 


_L  r£  +  b  +  c  v)  dt  =  ^ 
mo 


By  integrating  once,  we  have 


i_  .  -i-  =  —  »  exp  [-2-  f(t)1  =  f 

7  v  V  m  n 


......  ft  fa  .  v-)  dt.  This  completes  our  demonstration, 

where  f(t)  -  -  J  +  u  ^  VJ 

since  Equation  (31)  says  that  the  ve '““^^itially’wLr^he  con,"" 
ti*  t  are  always  in  the  rat  to  t  at  they  ^  rectiUnea 

ponents  were  v  ,  v  ana  vqz  ^ 

and  the  projectile  will  continue  along  the  line  of  vQ.  is  imp 

that  s  =  r,  the  position  vector  so  that 

needed  to  describe  the  motion.  A  simple  coordinate  r 

any  one  component  serve  as  well  as  another. 

Since  v  is  known  explicitly  as  a  function^  time^  -^als^find 
f(t)  explicitly  as  a  function  of  t  m  .  Equation 

tables  or  more  simply  by  ^^g^^usiir^iUai  “conditions  ‘ 
f31)  If  we  integrate  Equation  (31),  using 
z  =  y  =  z  =  s  =  r  =  0  at  t  =  0,  we  have 


x  _  _y_  = 


where  s  is  given  by  Equations  (24),  (26)  or  (28) • 
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vrtlWai!de.,milCertainly  P°ssible  «  “btain  explicit  solution,  for 

SB* 

(28)^edu^rRobertSOT’sSCaseedestinh  H° -"^u  that  Ecluations  (2?)  and 
b  =  0,  q  =  4  ac  0  ^  the  previous  section  for 

U  o,  c  >  0).  m  particular  they  become 


V 

£ 

c  tan  , 

tan  ^  ( v 

•/“  )  -  */ac  t/m 

1 

\  ° 

1 a/  0 

and 


(33) 


S  =  (m°/C)  tn  [C0S  |ta"_1  (“ o#)-  ^t/no|/cos  jta„-1(vo1/|)|], 

(34) 

At  the  perforation  limit  vq  -  v^,  we  require  v  =  0  when  s  .  T,  the  target 

iciness.  From  Equation  (33)  we  find  that  this  occurs  at  t„  =  (m  //^c) 

tan  (vq£  ST/ a).  if  we  use  this  value  of  t  in  Equation  (34)  and  °et 
s  -  T,  we  find 


oi 


/a/c  tan 


!»•■'  r--)l 


(35) 


which  I“iehe,^as  T  *  0  or  n,Q  w  -  and  increases  without  limit  as  T  ■*  « 


tan 


'  jcos-1 


-cT/m 


cos  <  cos 


(e'CT/"j| 


2cT/m 

-  1  =  e  °  -  1, 


(36) 


then  Equation  (35)  is  obviously  the  same  as  Robertson’s  v  ,  given 
above  after  Equation  (16).  When  perforation  occurs  for  v  >  v  ,  then 
1  -  tj.  v  =  yi  and  s  =  T.  From  Equation  (34)  we  can  find°t  or\he 

argument  of  the  cosine  involving  and  use  this  in  Equation  (33)  to 
obtain 
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Vl2  =  |  tan2  j  cos'^e  '  °  cos  J  tan"1^  ^  J  J 


•2cT/mQ  2 


/ 

rr  \  ) 

1 

Wt  - 

1 

\ u 

J 

-2cT/m 


rv2 

O 


vot] 


which  is  Robertson's  result  in  Equation  (Ib^above.  In  Equation  (37)  we 
have  used  Equation  (36)  and  the  identity  sec  x  1  t 

^  for  rectilinear  motion  during  zero  obliquity  impacts  the  force 


Even  for  rectilinear 

*  - 1-  c.  -  *  *  m2 phr™ung 

io  ^STih^iti-  when  J  vanishes  in  a  thick  enough  target.  Thus 

ricochet  can  be  described  since  the  s=  0,  t  <  ”  can 

0  <  s  <  T  for  t  <  »,  and  the  ricochet  “nditions,  v  ;  requires 

he  met._  However,  embed^nt^cannot  ^described^inc  ?  ^  ^  y  , 

If*aU=“o  because“a  -  a(v0)  vanishes,  then  embedment  can  be  described. 

However,  embedment  will  occur  only  for  a  particular  vQ  or  vo's  (depen  ing 
on  the  functional  form  of  a(»0)  )  instead  of  over  a  continuous  range  o 
v  as  observed.  We  could  define  a(vD)  to  be  zero  over  a  continuous 
range  of  V  ,  but  this  seems  somewhat  artificial.  If  a  a 

F  -  a.  -  b°v  -  c  v2  and  again  embedment  cannot  be  described.  Neither^  ^ 

unphysical  behavior  which  is  not  obsen  d.  J  only  if  a  -  0. 

requires  F  0,  v  -*•  u,  s  d  cT/m 


-Physical  behavior  which  is  not  ^  described  only  l(  a  .  o, 

requires  F  0,  v  -*•  u,  s  d  T/ 


Ill  c 

so  that  smav  .|ln(l‘5  v„t>  °r  V, 


-  1).  Robertson’s 
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case  (b  =  0)  then  implies  v^n  =  0  which  agrees  with  Equation  (35)  for 

a  =  0.  Since  finite  values  for  v  n  are  observed,  then  Robertson's  case 

o£ 

should  not  be  used.  A  ricochct/embedment  limit  requires  F  ->  0,  v  ->  0, 
s  -*  0  as  t  ->  00  and  cannot  be  described  since  the  requirement  a  >  0  pre¬ 
vents  F  ->  0.  If  we  define  a(v  )  to  be  zero  for  all  v  <  v  _  ,  then  we 

o  o  ol 

eliminate  all  possibility  of  ricochet  and  have  only  two  regions,  embed¬ 
ment  for  0  ^  v  ^  v  n  and  perforation  for  v  >  v  „  ,  Ricochet  or  rebound 
o  oi  1  o  oi 

after  penetration  for  low  speed,  zero  obliquity  impact  has  not  been  well 
explored  experimentally,  but  seems  possible. 

o 

In  summary  the  force  F  =  [-(a  +  b  v  +  c  v") ]  can  describe  only 
rectilinear  motion  and  so  is  limited  to  zero  obliquity  impacts.  Although 
perforation  can  always  be  described  for  thin  enough  targets,  and  ricochet 
can  be  described  if  a  >  0,  embedment  over  a  continuous  range  of  vQ  can 
only  be  described  by  defining  a(vn)  to  be  zero  over  this  range. 
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III.  FORCES  DEPENDENT  ON  DEPTH  OF  PENETRATION  AS  WELL  AS  SPEED 


As  we  have  seen,  the  constant  a  in  Equation  (15)  has  been  modified 
to  be  proportional  to  some  type  of  projectile  area  multiplied  by  a 
"shatter"  or  "hardness"  coefficient.  The  area  dependence  is  reasonable. 
In  addition  there  is  a  more  physical  interpretation  for  the  constant 
coefficient  of  proportionality.  Johnson7  points  ou£  the  usefulness  of 
dividing  impact  regimes  according  to  the  ratio  (p  v  /a)  which  is  approx¬ 
imately  equal  to  twice  the  crater  volume  for  simple  cases.  Here  o’  is 
the  mean  flow  stress  and  p  the  density.  This  suggests  adopting  the 
form 


a  =  a?  o  A  (38) 

rr*  _  _  _ 

so  that  Recht  and  Ipson's  W  =  /  a.  o  A  ds  =  a.  o  A  T,  which  is  constant 

o  1  z 

if  o’  and  A  are  constant.  They  argue  from  machining  studies  that  a  is  a 

constant,  while  A  is  a  constant  at  least  for  a  non-deforming  blunt 
cylinder  hitting  end-on  at  zero  obliquity.  For  most  metals,  except  soft 
ones  like  tin  and  lead,  their  plastic  shear  strength  is  not  very  depen¬ 
dent  on  shearing  speed14,  so  o  constant  is  a  good  approximation. 

However,  for  many  projectiles  A  depends  somewhat  on  depth  of  pene¬ 
tration,  although  a  constant  average  presented  area  like  that  in  Equa¬ 
tion  (38)  may  be  quite  adequate  to  represent  the  shearing  force.  For  a 
sphere  with  radius  R,  the  cross-sectional  area  at  the  original  target 

plate  impact  plane  is  -rr [R2  -  (R  -  s)2]  =  ir[2  R  s  -  s2]  where  s  is  the 

depth  of  penetration.  This  becomes  a  constant  it  R  for  s  >  R  which 
applies  for  the  rest  of  its  motion  through  a  thick  target.  For  a  thin 
target  it  may  be  more  appropriate  to  use  an  average  value  proportional 

to  tt  R2  for  A  in  Equation  (38).  For  a  cylinder  with  a  conical  nose  the 

2  2 

cross-sectional  area  at  the  impact  plane  is  tt  tan  <t>  s  where  1S  the^ 
constant  half-angle  of  the  cone  vertex.  This  becomes  the  constant  tt  R 
when  the  cylindrical  portion  of  the  projectile  (of  radius  R)  enters  the 
target. 

In  addition  to  shearing  forces  there  are  also  frictional  forces 
which  depend  on  the  contact  area  between  projectile  and  target  and  which 
may  be  given  the  simple  form 


F 


f 


*f  P  Ac 


(39) 


14P.  W.  Bridgman ,  Studies  in  Large  Plastic  Flow  and  Fraoture,  N.Y. , 
McGrow  Hitt ,  1982,  p.  284. 
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where  Uf  is  a  constant  coefficient  of  sliding  friction  and  P  is  the 
pressure.  In  order  to  write  this  simple  form  we  must  assume  that  the 
stress  is  isotropic  as  for  a  fluid,  an  approximation  which  should  improve 
with  higher  speeds.  If  we  retain  all  components  of  the  general  stress 
tensor  we  will  obtain  a  much  more  complicated  form  which  is  probably 
not  necessary  for  our  purpose.  Here  we  will  take  P  to  have  a  constant 
average  value  during  the  motion.  However,  we  will  take  into  account  the 
variation  of  A  with  depth  of  penetration  s.  For  a  sphere  the  area  of 

o 

the  contact  surface  is  2  n  R  s,  which  becomes  2  tt  if  for  s  >  R.  For  a 

i  ?  1  /  ? 

conical  nose  at  zero  obliquity  the  contact  area  is  tt  R  (R^  +  s  ) 
where  R  is  the  radius  of  the  cone  base  as  well  as  the  radius  of  the 
cylindrical  body.  When  the  body  of  this  projectile  enters 

the  target  the  nose  contact  area  becomes  the  constant  rr  R  (R^  +  h^)^^ 
where  h  is  the  cone  height  and  an  additional  cylindrical  area  2  tt  R  s  of 
the  body  comes  into  contact  with  the  target.  In  general  A  might  be 

reasonably  represented  by  the  form  Aq=  Aj  +  Cj  s  +  C  s2  with  ,  C  and 

C2  constant  and  usually  with  C2  <  C  . 

For  oblique  impact  and  for  irregular  projectile  shapes  like  those 
coming  from  naturally  fragmenting  warheads  or  armor  spall  and  debris, 
the  areas  important  for  shear  and  friction  forces  will  be  more  compli¬ 
cated,  although  a  quadratic  form  in  s  is  probably  still  adequate  for  an 
average  contact  area.  Projectile  deformation  will  increase  the  contact 
area  as  s  increases.  Tumbling  and  yawed  impact  present  additional 
problems  while  projectile  erosion  further  complicates  matters.  Frosion 
will  be  considered  in  a  later  section  of  this  report. 

For  constant  mass  projectiles  we  will  adopt  a  force  field  which 

depends  on  depth  of  penetration  in  a  simple  way  namely  -  (a  +  d.  s  +  d0 
2  o  1  2 

s  ),  where  a^,  d1  and  d9  are  constants.  This  form  is  suggested  by  the 

shear  stress,  friction  and  deformation  considerations  given  above.  It 
seems  to  be  a  reasonable  way  to  represent  the  "plug  shear”,  ’’penetrat  i  on” 
and  "deformation”  processes  described  by  Recht  and  Ipson  in  their  1974 
paper15.  The  ”bulging  or  dishing”  response  of  a  thin  target  plate  also 
depends  on  how  far  the  projectile  tip  has  progressed  beyond  the  original 
impact  plane.  An  initial  elastic  response  dependent  on  the  first  power 
of  s  is  expected,  followed  rapidly  by  an  inelastic  response  dependent 
on  higher  powers  of  s,  much  as  in  the  form  we  have  chosen.  However,  the 
"breeching  or  petaling"  mode  of  plate  failure  described  by  these  authors 
is  not  expected  to  depend  on  depth  of  penetration.  Instead  it  seems  to 
be  dependent  on  speed  and  will  be  included  in  the  speed-dependent  part 
of  our  force  field.  It  is  of  more  importance  for  thin  plates  and  high 


15/?.  F.  Recht  and  T.  W.  Ipsony  n Ballistic  Penetration  Resistance  and  Its 
Measurement y  fr  Proceedings  of  the  First  International  Symposium  on 
BallisticSy  Orlandoy  Floriday  1974 . 
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obliquities  since  the  plugging  mode  of  failure  is  more  common  near  zero 
obliquity.  For  example,  thin  aluminum  plates  fail  by  petaling  when 

>  30°  while  thin  steel  plate  failure  has  been  described  as  "hinging"11. 

When  a  body  moves  through  a  fluid  its  motion  is  opposed  by  a  viscous 
or  internal  friction  force  which  is  proportional  to  the  first  power  of 
the  speed.  This  is  in  addition  to  the  sliding  friction  force  mentioned 
above  which  depends  only  on  the  contact  area  and  pressure.  For  a  sphere 
of  radius  R  Stokes'  force  is  6ttRijv,  where  y  is  the  coefficient  of 
viscosity.  Thus  the  constant  b  in  Equation  (15)  can  be  given  a  simple 
physical  interpretation.  The  viscosity  of  solids  is  a  well-known 
physical  concept16  and  for  common  metals  it  has  been  measured  at  high 
speeds1 7 . 


For  a  body  moving  through  a  fluid  another  force  proportional  to  the 
square  of  the  speed  called  the  drag  force  is  also  considered.  Since  we 
expect  metal-metal  projectile-target  penetration  to  be  a  non-linear 
process  we  should  also  consider  a  force  dependent  on  v2  as  well  as  v. 
This  could  include  such  non-linear  processes  as  conversion  of  kinetic 
energy  into  heat  energy  (or  even  radiant  electromagnetic  and  mechanical 
energy).  For  example,  Johnson7  estimates  the  temperature  rise  in  a 
penetration  to  be  v2/(2  Cn  J)  where  J  is  Joule's  mechanical  equivalent 


of  heat.  Using  =  0.1  cal/(gm-deg)  for  the  specific  heat  of  steel  or 

aluminum,  we  calculate  a  temperature  rise  of  about  100°C  when  v  is  0.3 
km/sec.  In  addition,  a  v2  term  could  represent  the  energy  transferred 
to  pieces  which  are  separated  from  the  target  such  as  the  plug  sheared 
out  at  low  obliquity,  the  petaling  or  hinging  mentioned  above  and  the 
energy  used  in  accelerating  ejecta  from  the  impact  side  of  a  target 
plate  as  well  as  assorted  bits  of  spall  and  other  debris  emerging  from 
the  rear.  This  is  the  usual  inertial  reaction  interpretation. 


Our  complete  force  field  is  then  of  the  foim 


F  =  -  fa  +  b,  v  +  b„  v2  +  d,  s  +  d^  s2] .  (40) 

‘01  z  i  z 

Admittedly  this  is  an  oversimplification  which  contains  for  example  no 
cross  product  terms  in  s  and  v  and  no  explicit  functions  of  pressure, 
temperature,  time  or  other  factors  which  might  be  included.  However, 
we  have  given  simple  physical  interpretations  for  the  coefficients 


16 E .  A.  Evans,  "Viscosity  of  Solids,  "  in  American  Institute  of  PJnjsics 
Handbook,  N.Y.,  McGraw-Hill,  1957,  c.  2h. 

17 S.  K.  Godunov,  A.  A .  Deribas  and  V.  I.  Mali,  "Influence  of  Material 
Viscosity  on  the  Jet  Formation  Process  During  Collisions  of  Metal 
Plates,"  Fiz.  Goreniya  i  Vzryva  If  (1)  3,  1975,  and  V.  I.  Mali,  V.  V. 
P ai  and  A.  J.  Skovpin,  " Investigation  of  the  Breakdown  of  Flat  Jets,  " 
Ibid.,  10  (15),  775,  1974. 
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a  ,  b.  ,  bn,  cL  and  dn  which  should  eventually  enable  us  to  estimate  their 
o  1  2  1  2 

values  from  experimental  data  obtained  in  studies  other  than  penetration 
studies.  If  we  can  use  such  values  and  successfully  predict  the  outcome 
of  penetration  experiments,  then  we  have  gained  some  fundamental  insight 
into  the  penetration  process,  since  these  coefficients  will  no  longer 
be  purely  empirical  constants.  If  we  cannot  do  this,  we  will  have  to 
be  content  with  calling  them  empirical  constants.  Still,  if  we  can 
adjust  them  to  known  data  and  successfully  predict  the  unknown,  we  will 
in  any  case  have  a  relatively  simple  method  for  solving  practical 
problems  without  a  large  expenditure  of  time  and  money. 

A  more  general  form  than  Equation  (40)  which  includes  higher  powers 
of  s  and  v  as  well  as  cross  product  terms  was  attributed  to  Thompson18 
by  Murff  and  Coyle19,  although  in  practice  they  used  Equation  (40)  with 
^2  =  ^2  =  ^  anc^  f°und  it  adequate  to  describe  embedment  in  semi-inf inite 

soil  targets.  Their  suggestion  of  an  infinite  series  of  terms  (s11  vm) 
did  not  include  suggestions  for  physically  interpreting  these  terms  and 
so  would  reduce  the  coefficients  of  these  terms  to  empirical  constants. 


18L.  G.  Thompson,  et  at. ,  "The  Effect  of  Soil  Parameters  on  Earth  Pene¬ 
tration  of  Projectiles,"  Technical  Report  to  Sandia  Laboratories, 

Texas  ASM  Research  Foundation,  July  1969 . 

1 V.  D.  Murff  and  H.  M.  Coyle,  " Prediction  Method  for  Projectile  Penetra¬ 
tion,"  J .  Soil  Mech .  Foundations  Division,  Proc.  ASCE  99,  1033,  1973 . 
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IV.  SPECIAL  AND  GENERAL  SOLUTIONS  DESCRIBING  THE  MOTION 


As  we  have  seen  in  Section  II,  Equation  (15)  as  it  stands  permits 
only  rectilinear  motion.  However,  we  can  write  it  in  a  component  form 
which  will  permit  curvilinear  motion,  provided  we  require  the  coefficients 
in  our  force  field  to  be  anisotropic.  This  is  a  reasonable  thing  to  do 
in  any  case  even  for  a  semi-infinite  target  since  the  forces  parallel 
and  perpendicular  to  the  impact  surface  at  the  point  of  impact  must  be 
different.  Thus,  even  for  small  obliquities  the  trajectory  will  be 
curved.  This  is  because  the  impact  surface  is  an  interface  between 
materials  of  different  properties,  commonly  metal  and  air  in  cases  of 
Army  interest.  If  the  impact  surface  is  planar  and  effectively  of 
infinite  extent  but  the  target  thickness  T  is  finite,  this  will  have  a 
strong  anisotropic  influence  too.  The  message  that  a  projectile  has 
struck  the  front  face  of  a  target  plate  is  rapidly  transmitted  by  shock 
wave  to  the  rear  surface  where  transmission  and  reflection  occur.  Be¬ 
cause  of  this  the  resistance  offered  by  the  plate  perpendicular  to  the 
surface  plane  is  different  from  that  parallel  to  this  plane.  In  the 
special  case  where  d^  =  =  0  in  Equation  (40),  Newton's  law  becomes 


m  x  +  a  +  b,  x  +  b0  (x)  =0 

o  ox  lx  2x 


(41) 


and 


m  z  +  a  +  b.  z  +  b-  (z)  =  0 

o  oz  lz  2z 


(42) 


where  the  y-component  is  ignorable  because  of  the  symmetry  of  an 
oblique  impact  on  a  plate  of  thickness  z  =  T.  If  we  allowed  aQx  =  aQZ> 

b  =  b  and  b„  =  b„  we  would  be  right  back  to  Equation  (15)  with  its 
lx  lz  2x  2z 

rectilinear  motion  and  other  limitations.  Now,  however,  we  can  find 
explicit  solutions  for  x(t)  and  z(t)  which  have  the  same  form  as  Equa¬ 
tions  (24),  (26)  and  (28)  but  which  are  not  in  the  same  ratio  for  all 

time  since  a  #  a  ,  b,  *  b,  and  b„  *  b~  .  Thus  curvilinear  motion 
ox  oz  lx  lz  2x  2z 

can  be  described.  However,  both  ricochet  and  embedment  cannot  be 

described  without  artificial  definitions  for  a  .  We  can  describe 

O  Zr 

either  perforation  and  ricochet  or  perforation  and  embedment  but  not 
perforation,  ricochet  and  embedment,  all  three  as  a  continuous  function 
of  vq  for  given  obliquity. 

Now  consider  the  special  case  in  which  -  d^  =  0  in  Equation  (40)  , 
leading  to 


m  x  +  a  +  b..  x  +  d.  x  =  0 
o  ox  lx  lx 


and 


m  z  +  a  +  b,  z  +  d.  z  =  0 
o  oz  lz  lz 


(43) 

(44) 
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for  a  plate  target.  Not  only  is  curvilinear  motion  possible  because  of 
the  anisotropy  of  the  force  field  coefficients,  but  ricochet,  embedment 
and  perforation  can  all  be  described  as  a  single  function  of  v^.  The 

solutions  for  each  Equation  (43)  and  (44)  are  well  known  since  we 
recognize  this  differential  equation  form  as  that  of  the  damped  harmonic 
oscillator  with  constant  imposed  force  a  .  This  force  we  have  already 

interpreted  as  that  needed  to  remove  enough  material  from  the  target  to 
make  way  for  the  passage  of  the  projectile.  If  s  now  represents  either 
x  or  z  (or  y  also  in  the  general  case),  we  can  write  the  solution  for 
each  component  as  the  superposition 

Ys+  1  Ys"  1 

s  =  Als  e  +  Bls  e  +  5S  (45) 


where  6  =  -a  /d  ,  y  +  =  -  a  ±  6  ,  cx  =  b.  /  (2  m  )  and  6  = 

2s  os  is  s  s  s  s  Is  o  s 

Cots  '  d^s/mQ)  .  If  we  take  the  time  derivatives  of  the  component 

displacements  and  impose  the  initial  conditions  x  =  z  =  0,  x  =  v  and 

ox 

z  =  v  at  t  =  0,  we  have 
oz 


for  each  component.  The  usual  cases  for  real,  zero  and  imaginary  3^ 

will  be  studied  in  the  next  section.  For  the  moment  we  will  merely 
point  out  that  ricochet,  embedment  and  perforation  can  all  be  described 
as  a  single  function  of  v  . 

As  we  have  seen,  part  of  the  constant  a  consists  of  the  term 

os 

a^s  o'  which  should  depend  on  the  target  thickness  since  a  thicker 

target  requires  a  greater  force  to  shear  out  more  material.  However, 
this  increase  does  not  continue  indefinitely  and  for  semi - inf inite 
targets  the  influence  of  the  rear  surface  interface  is  never  felt  so 
the  force  is  T-inaependent .  Such  behavior  can  be  represented  by  using 
components  of  the  line-of-sight  thickness  L  =  T/cos  0  ,  namely, 
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L  =  L  sin  6  =  T  tan  9  and  L  =  cos  0  =  T  in  a  form  like  a2g  o  Ag  = 

x  O  O  L  u 

f(L  )  =(als/a3s)  Cl  -  e’a3s  Ls) ,  (s  =  X,  z)  .  Thus  for  thin  targets 
f  (LS)  j,  a  5  L  while  for  semi-infinite  targets  (Ls  ->■  “) ,  fd-s)  *  (als/a3s) 

r-JM  the  form 


COTlbLclllL  i - r  .  -x  T 

Mf  PA  v  =  C  V  .  Thus  we  can  write  aQs  =  *v-s,+  cls  vqs  -  aJs  ^ 

+  r  v  CSfor  thin  targets.  For  any  target  we  can  write  the  z-component 
Is  os 

of  the  force  as 

Fz  »  -  If(Lz)  *  cu  voz  +  ^lz  dlz  Zl  (48) 

.  t  _  n  Tn  addition,  the  influence  of  target  thick - 
which Lls=no  h^r°een  in'cluded  in  the  force  governing  the  motion.  This 

form  also  enables  us  to  require  the  "the  ri^ochet'Ibed- 

case  of  embedment,  including  the  two  lt«p"|  cases,  ..  Umit 

ment  limit^nd^he^erforation/embedment^mi.  ^  becomes 

o  oRL 


0  =  -  [fCT)  *  c  v  ] 


while  at  the  perforation  limit  (when  v0  ■  voPL)  z  =  0,  z  =  T,  and 


0  =  -  [f  (T)  *  c,z  vozpL  *  du  T], 


These  equations  may  be  solved  simultaneously  to  obtain  f(T)  -  dizTvozRL/ 

(v  PT  "  vozRL)  and  Clz  “  '  dlz  T/<:VozPL  "  VozRL5  S°  tl,at  f<T>  3nd  lz 
bo^  depend  ™  d,,.  The  constant  term  in  Equation  (45)  then  becomes 


«z  ■  -  aoz/dlz  -  -  [f(T)  *  Clz  Vozl/dlz 


=  T  <voz  -  V0ZRL)/(V0ZPL  '  VozRL) 


r  ,  ,  Since  the  physical  requirement  of 

which  is  negative  for  vqz  <  vqzRL.  ' 

finite  s  and  s  requires 


Y  <  y  <  0 

S  S 
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in  Equation  (45),  at  least  for  their  real  parts,  we  see  that  z  tends 
eventually  to  be  negative  in  the  ricochet  regime.  When  z  returns  to 
zero,  the  problem  ceases  and  the  exit  velocity  component  z  will  be 
negative.  At  the  ricochet  limit  6^  =  0  and  z-^0,  z  0  as  t  For 

VozRL  <  Voz  <  VozPL>  T  >  5Z  >  0  and  embedment  occurs  since  z _  <  T. 


For  v  >  v  rkT 
oz  ozPL 


max 


,  6  >  T,  and  perforation  occurs. 


have 


If  we  substitute  [f(T)  +  c,  v  ]  =  -  d  5  in  Equation  (48),  we 

i  Z  O  Z  I  z  z 


Fz  =  -  [bU  1  4  dlz  <z  •  6z>] 


(53) 


where  6^  is  seen  to  play  the  role  of  the  eventual  value  of  z  as  t  be- 

comes  arbitrarily  large  and  F  -*  0,  z  -*  0.  Because  of  this  6  will  be 
called  a  final  position  parameter.  z 

Similarly,  for  the  x-component 

[f(L  )  +  c 


F  = 
x 


,  v  +  b,  x  +  d-  x] 
lx  ox  lx  lx  J 


(54) 


We  also  require  ->  0,  x  0  at  the  ricochet  limit,  so 


°  =  '  [f(V  4  clx  VxRL  4  dlx  XRL> 


and 


(55) 


°  =  ■  [f(Lx)  4  clx  voxPL  4  dlx  XPI,1 


(56) 


at  the  perforation  limit  where  x  also  vanishes.  The  limiting  values 
XRL  an<^  XPL  are  no *  usuaHy  known.  From  these  relations  we  find 


f(V  ^lx  ^XRL  VoxPL 


"VPL  v oxRLJ '  v  voxPL  voxRLJ 


clx  =  -  dlx  (XPL  -  xRL’/f''oxPL 


v  .)  which  may  be  substituted  into 

OXKL 


Equation  (54)  to  obtain 

Fx  =  -  tblx  *  +  d 


lx  (X  -  V1 


(57) 


where 


[XPL  VoxRL  '  XRL  VoxPL  "  (xPL  "  XRL]  Vox]/(voxPL  '  VOXRL)  (58) 


plays  a  role  similar  to  6^.  Of  course  6^  in  Equation  (51)  can  be  cal¬ 
culated  from  experimental  data  for  a  given  6^  to  within  the  accuracy 

possible  for  these  types  of  limit  measurements.  However,  the  limit 

values  in  Equation  (58)  are  not  usually  measured  and  it  is  necessary  to 

estimate  6  from  its  definition.  This  can  also  be  done  for  6  . 

x  z 
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In  summary,  we  have  introduced  a  soluble  set  of  equations  (Equations 

(43)  and  (44))  and  pointed  out  their  ability  to  describe  curvilinear 
motion,  including  the  extremes  of  ricochet,  embedment  and  perforation. 

A  discussion  of  the  constant  force  term  a  (composed  of  shear  and 
frictional  parts)  led  us  to  the  general  position  parameter  8^  given  in 

Equation  (51)  and  its  x-analog  in  Equation  (58) .  In  this  country  it  is 
customary  to  determine  the  perforation  limit  speed  v0ZpL  for  given  target 

thickness,  and  from  Equation  (51)  we  see  that  8^  =  T  when  vqz  =  vozpL' 

From  Equation  (45)  with  s  =  z,  we  also  see  that  z  =  6^  =  T  (since  t  -*■  °°)  . 

In  Germany2®  it  is  common  to  determine  that  thickness  of  target  plate 
which  will  just  stop  the  projectile  for  given  striking  speed.  This  limit 
thickness,  T  ,  of  course  corresponds  to  the  limit  speed,  vozpL>  so  a8ain> 

from  Equations  (45)  and  (51)  z  =  (6^,)^  =  T  =  when  vqz  =  vozpL*  f°r 

an  effectively  semi-infinite  target  neither  vQzpL  nor  T£  is  ever  reached 

experimentally.  Thus,  even  for  the  highest  striking  velocity  vqz  <  vozp^ 

and  by  Equations  (45)  and  (51)  we  have  zro  =  (6^  <  T  =  T^.  In  other 

words,  the  maximum  depth  of  penetration  in  a  semi-infinite  target  is  less 
than  the  limit  thickness.  Experimentally20  it  is  found  that  the  ratio 
T„/(6  )  decreases  as  v_  increases  and  seems  to  approach  1.5  as  vq 

approaches  the  longitudinal  sound  speed  in  the  target  material  (at  least 
for  steel  targets  and  short  steel  rods).  Because  of  spalling,  dishing, 

petal ling  and  other  finite  target  effects  as  well  as  projectile  erosion 

and  breakup,  T„  is  roughly  an  exponential  function  of  vq  in  the  ordnance 
range . 

If  we  allow  b2s  and/or  d2s  to  be  non- zero  in  the  components  of 

Equation  (40),  then  we  have  a  set  of  non-linear  equations.  The  general 

closed- form  solution  to  any  one  of  these  equations  is  not  known  at 
present  and  will  be  the  subject  of  future  work.  If  need  be,  we  can 
always  write  approximate  solutions  in  series  form.  However,  we  will, 
first  try  to  answer  the  question  of  how  well  we  can  describe  the  avail¬ 
able  data  for  ordnance  velocity  penetrations  using  Equations  (43)  and 

(44) .  These  are  the  simplest  forms  we  can  expect  to  use  to  describe 
curvilinear  motion  resulting  in  perforation,  embedment  or  ricochet 
depending  in  a  continuous  manner  on  v^  and 


20V.  Hohler,  E.  Schneider  and  A.  Stilp ,  ", Endballistik  von  Splittem  I-V, 
Beriahte  E  6/73,  E4/74 ,  E4/75,  E9/75  and  E  7/75  (1972-1975)  Emst- 
Maah-Instituty  Freiburg. 


37 


V.  ILLUSTRATIONS  OF  MOTIONS  IN  THE  SIMPLEST  FORCE  FIELD 

In  this  section  we  wish  to  illustrate  some  of  the  capabilities  of 
the  force  with  components  Fs  =  -  (as  +  s  +  dg  s) .  We* will  not  give 

an  exhaustive  analysis.  Rather  we  will  show  the  type  of  agreement  with 
experiment  which  might  be  expected  by  showing  this  agreement  in  a  few 
simple  cases.  For  convenience  we  reproduce  the  z -component  of  the 
solution  for  the  equation  of  motion  as  well  as  its  time  derivative, 
dropping  numerical  subscripts: 


z  = 


(Y_*  1  ~  UV02  *  6Z  V’  °  2  -  <VOZ  *  4Z  V*6'2  1  *  5Z 


z  z 


(59) 


z  = 


(vz  -  Yz  > 


+ 

[(voz  +  6Z  Yz’)  YZ+  e  Z  -  (voz  +  Vz  +  )yz~ 


Y  t 

-  (v  +  6  y  )y  eZ  ] 

(60) 


where  y  =  -  a  ±  @  a  =  ,  0 

z  z  t  z  2m  z 


a  ^  -  d  / m  and 
z  z  o 


sz  -  -  T  (vozRL  -  \,z)/(VozPL  -  VozRJ.)  =  ‘  °„z/dz' 

As  is  well  known,  can  be  real,  zero  or  imaginary.  If  B_  is  real  or 

^2  z 

zero,  then  d^  <  mQ  .  For  d^  =  0  we  have  a  solution  analogous  to  that 

wnicn  describes  a  body  falling  through  a  viscous  medium  in  the  earth Ts 
gravitational  field,  aoz  =  However  we  have  already  seen  that  such 

a  force  field  dependent  only  on  the  velocity  is  not  adequate  for  describ¬ 
ing  an  embedment  region  between  ricochet  and  perforation  regions  in  a 
plot  of  exit  speed  versus  striking  speed.  Experimentally,  velocity 
decreases  with  time  which  requires 


Yz  <  Yz  <0  (62) 

W 

so  that  z  and  z  will  always  be  finite  in  Equations  (59)  and  (60).  The 
same  is  true  of  course  for  the  other  components.  Equation  (62)  implies 

az  ^  dz/m0  >  °»  since  dz  <  0  would  make  Bz  >  and  >  0,  leading 

to  a  non-physical  velocity  increase  with  time.  The  special  case 
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m  a  or  B 
o  z  z 


=  0  requires  the  degenerate  solution 


a  6  )  t  -  6  ] 
z  z  z 


z  =  e 


r  -  (v  *  a  6  J  “z  • 
oz  oz  z  z  ^ 


By  rearranging  Equations  (59)  and  (60)  we  can  find 


Y„  t 


"  _  A  11/fv  +6  Y  ) 


e  "  =  iz  -  Yz  i*  ’  VJ/  v ¥ oz  'z  'z 


e  Z  =  [z  -  Yz+  (z  ‘  6zn/(voz  +  6z  Yz  ).  ( 

If  we  divide  Equation  (65)  by  Equation  (66)  and  solve  for  t,  we  find 


)1  fv  +  &  Y  )  ) 


lz  -  Yz  IZ  -  ozJj  .voz  •  ~z 


(Y  +  -  Yz”)  |  [*  -  ^z+  (z  -  6z)]  (Voz  +  6z  Yz  ^ 


In  case  3. 


=  0  we  find  from  Equations  (63)  and  (64) 


+  (z  -  sz)]/(Voz  ‘  az  V 


t  =  [V  (Z  -  6  )  +  6  z]/[(z  +  a 

1  OZ  z  z 


u  -  «z)>  <y01  -  az  V1 


=  1-  tn  ((v„  -  a.  «,)/<;  -  5z),] 


Now  let  us  illustrate  the  app 


analogous  x-component  s 


thP  auulication  of  this  solution  (and  the 
olut ion)1 by  using  some  of  Backman  and  nnnegans 
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data11.  Figures  1,  2  and  3  reproduce  some  of  their  data  for  a  1  gram 
hard  steel  (SAE  52100)  sphere  penetrating  a  6061  aluminum  alloy  plate 
of  thickness  T  =  9,53mm.  The  sphere  diameter  was  6.35mm.  As  before, 
we  will  take  our  z  axis  anti-parallel  to  the  normal  to  the  impact 
surface  and  our  origin  at  the  point  of  impact  with  the  motion  occuring 
in  the  x-z  plane.  We  note  that  our  choice  of  axes  makes  L  =  0  =  v 

y  oy 

so  6  =  0  and  y  =  y  =  F  =  0  for  all  t.  By  convention  the  obliquity  6 

j  y  o 

is  the  angle  measured  from  the  impact  face  normal  to  the  negative  of  the 
striking  velocity  vector.  In  order  to  conform  to  the  conventions  of 
Backman  and  Finnegan  we  will  measure  the  exit  angle  G  =  arc  tan  (x/z) 
from  the  negative  of  the  impact  face  normal,  that  is,  from  the  z  axis, 
to  the  exit  velocity  vector.  Thus  for  perforations  we  will  have 
0  <  0  <  90°,  for  ricochets  we  will  have  0  >  90°,  while  for  embedments 

G  is  undefined.  An  upper  limit  of  180°  for  0  is  suggested  by  the  avail¬ 

able  data  near  the  ricochet/embedment  limit,  but  the  data  is  not  good 

enough  to  be  sure.  There  are  also  indications  that  0  -*  180°  -  G  as 

o 

VQ  4  0  as  we  might  expect  since  we  are  approaching  the  case  of  elastic 

rebound  (no  penetration)  for  which  the  angle  of  reflection  is  equal  to 
the  angle  of  incidence.  Again  experimental  evidence  is  lacking  near 
this  limit  to  confirm  or  deny  such  behavior. 

At  the  other  extreme  of  very  large  striking  speed  (v  °°)  we  expect 

the  exit  speed  v  -+  v^  and  the  exit  angle  G  -►  0  .  This  expectation  has 

been  verified  by  experiments  with  very  thin  (aluminum  foil)  targets  and 

ordnance  projectiles  which  move  almost  as  if  the  target  were  not  there. 

In  both  extremes  (vq  0  and  vq  °°)  the  time  spent  penetrating  the 

target  is  very  small.  That  is,  the  ricochet  time,  tR,  which  is  the 

time  it  takes  for  z  to  return  to  zero,  and  the  perforation  time,  t  , 

which  is  the  time  it  takes  for  z  to  reach  the  value  T,  both  approach 
zero  for  these  extreme  striking  speeds.  It  is  clear  from  Equation  (60) 
and  its  x-analog  that^the  solution  behaves  properly  at  these  two 
extremes,  namely, x  / 1  z  |->  v  /v  so  0  =  arc  tan  (x  /z  )  0  =  arc  tan 

r  P  OX  OZ  p  p  O 

(vox/vqz)  for  very  large  vq,  while  0  180°  -  0q  for  very  small  vq  (small 

zn<  0  and  xn  >  0)  .  For  intermediate  values  of  v  which  lead  to  ricochet, 
R  R  o 

0  >  180°  is  a  possibility  (zR  <  0,  xR  <  0),  but  there  seems  to  be  no 

published  evidence  for  this  case.  For  zero  obliquity  impacts,  0  =  0  = 

sin  G  ,  x  =  x  =  0,  v  =  v  sin  9  =0  and  L  =  L  sin  0  =  0  for 
o  ox  o  o  x  o  1 

-a„  L 

SY  Y 

(1  -  e  )  =  0]  so  Fy  =  0  in  Equation  (54).  As  0^  -*  90°,  grazing 

x  o 

incidence  turns  into  a  near  mass  and  no  penetration  occurs. 
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Figure  3.  Measured  and  calculated  exit  angles  for  steel  spheres  vs  aluminum  plates  at 
45°  obliquity. 


For  real  Bz,  z  can  vanish  only  twice.  The  first  time  can  be  found 
by  letting  £  =  0  in  either  Equation  (67)  or  Equation  (69)  and  then 
z  =  zmnx  Can  f°und  by  using  this  value  of  t  in  either  Equation  (59) 

or  Equation  (63).  The  second  zero  of  z  occurs  as  t  -►  <*>.  It  is  of 
course  possible  for  t  ->  00  to  be  the  only  time  that  z  vanishes.  This 
occurs  at  perforation  limits  or  for  embedment  in  semi -infinite  targets. 

£ 

At  the  ricochet  limit  z  has  already  vanished  once  at  the  turning 

point  z  =  z  <  T.  For  this  limit  it  vanishes  again  as  z  -*  0  and 
max 

t  00 .  From  Equations  (67)  or  (69)  we  see  that  this  requires  6^  =  0 

which  is  insured  by  Equations  (51)  or  (61)  when  v  =  v  .  At  the 

oz  ozRL 

perforation  limit  v_  =  vozPL  and  62  =  T  in  E9uaticm  (61)  so  that  z 

vanishing  and  z  approaching  T  makes  t  -*  00  here  too.  For  embedment  in 

a  finite  target,  £  -*■  0  and  z->6  <  T  makes  t  ->  °°,  For  a  semi-infinite 

z 

target  which  guarantees  embedment  even  for  the  largest  vq  we  can  realize, 

z  =  as  t  ->  °°.  Since  v  D  increases  without  limit  as  T  increases 
max  z  oz  r  L 

without  limit  in  Equation  (51),  6  =  z  depends  on  v  . 

z  max  r  oz 

Figure  4  illustrates  schematically  some  of  the  trajectories  which 
can  occur  during  target  penetration.  These  trajectories  become  effec¬ 
tively  straight  lines  when  the  projectile  exits  from  the  target  into  a 
low  resistance  medium  like  the  air.  Alternative  trajectories  are 
illustrated  by  solid  or  dashed  lines  while  dotted  lines  indicate  pro¬ 
jections  of  trajectories  if  the  motion  had  continued  or  if  a  different 
distribution  of  target  material  had  been  present.  For  example,  if  there 
had  been  more  target  material  in  the  upper  left  of  Figure  4(a),  then  the 
x  coordinate  might  reach  a  maximum  in  finite  time  and  eventually  turn 
negative.  The  solid  line  in  Figure  4(a)  shows  a  typical  ricochet  while 
the  dashed  line  illustrates  a  possibility  which  does  not  seem  to  have 
been  documented  in  the  literature  as  mentioned  above.  In  the  dashed 
case  illustrated  the  exit  velocity  is  anti -paral lei  to  the  striking 
velocity  and  0  =  180°  +  6^.  Smaller  or  larger  values  for  0  are  possible 

with  the  latter  involving  a  trajectory  which  crosses  itself.  This 
crossing  could  occur  inside  or  outside  the  target.  Figure  4(b)  illus¬ 
trates  possible  trajectories  when  v  =  v  ...  ,  the  ricochet  limit.  Here 

o  oRL 

Vox  =  VoxRL  am*  ^x  =  XRL  *n  Equation  (56)  so  x  ->  6^  and  x  ■+  0  makes 

t  ->■  “  in  the  x-component  analogs  of  Equations  (67)  or  (69).  If  x 

remains  positive  and  approaches  zero  only  as  t  +  then  x  .  is  a 

RL 

maximum  value  which  x  approaches  only  as  t  00  and  the  exit  angle  0O1  = 

180°  as  illustrated  by  the  solid  line  in  Figure  4(b).  If  x  turns  nega¬ 
tive  in  finite  time  and  then  approaches  zero  as  t  -*■  ",  then  x  <  x 

RL  max 

and  0j^  >  180°  as  shown  by  the  dashed  line  in  Figure  4(b)  which  is  the 
the  limiting  case  of  a  dashed  line  in  Figure  4(a).  Figure  4(c) 


44 


illustrates  by  a  solid  line  an  embedment  trajectory  for  which  v  < 

oRL 

vo  <  v0pl  anc*  z  approaches  a  maximum  value  as  t  ->■  °°,  z  ^  0.  Of  course 
z  +  0  as  t  +  “  even  if  the  final  z  value  £  z  ,  (solid).  The  dashed  lines 

ITlciX 

in  Figure  4(c)  illustrate  such  possibilities.  At  the  perforation  limit 
the  available  data  suggest  that  the  exit  angle  is  close  to  zero.  The 
solid  line  in  Figure  4(d)  illustrates  this  and,  if  true,  then  xm  is  a 

maximum.  Since  x^^  0  as  t  co,  the  other  possibility  is  that  Xp^  < 

Xmax  ^  x  a^rea<^ir  vanished  once  in  finite  time.  This  is  shown  by 

one  dashed  curve  in  Figure  4(d).  Of  course  z  ->  0  and  z  ->  6  =  T  when 

z 

Vo  =  VoPL  wil1  make  t  ~y  00  in  Equations  (67)  or  (69).  Finally,  a  typical 

perforation  trajectory  with  the  exit  angle  0  <  0^  is  shown  by  the  solid 

line  in  Figure  4(e),  As  v  °°,  of  course,  0  ->  0  .  Another  possibility 

o  o 

shown  by  a  dashed  line  is  that  x  vanish  while  z  <  T  so  that  x  occurs 

max 

inside  the  target  and  the  exit  angle_0  <0.  A  third  possibility  in 
which  a  point  of  inflection  (x  =  0,  x  >  0)  occurs  inside  the  target  in 
finite  time  could  lead  to  0  >  0^  and  is  illustrated  by  another  dashed 

line  in  Figure  4(e).  A  similar  type  of  behavior  is  possible  for  ricochets 
also  but  was  not  illustrated  in  Figure  4(a)  to  avoid  crowding  and  confu¬ 
sion.  In  short,  a  variety  of  trajectories  can  be  described  by  Equation 
(59)  and  its  x-analog  although  not  all  of  them  may  be  realizable  for 
particular  projectile  and  target  material  combinations.  For  targets 
other  than  plates  of  homogeneous  material  or  for  projectiles  with 
special  shapes  or  motions  at  impact,  it  may  be  necessary  to  invoke  un¬ 
usual  trajectories  or  to  include  x- ,  y-  and  z-components  of  motion  for 
an  adequate  description.  We  also  note  that  for  real  3^  >  0  with  s  =  x, 

z  (or  y  if  needed)  all  combinations  such  as  3  =  0,  3  >  0  or  3  =3  = 

x  z  X  z 

0  or  B  >  0,  8  >  0,  etc.  are  allowed, 

x  z 

In  summary,  there  are  some  constraints  imposed  on  the  parameters 
an?.  >  and  d^  (for  s  =  x,  z)  from  general  considerations  of  plate  target 

penetration  with  Bs  >  0.  From  Equation  (62)  and  its  x  analog  we  have 

b  >  0,  b  >  0,  0  <  d  <  m  [b  /(2  m  )]"  and  0  <  d  <  m  [b  /(2  m  )]2, 

A  ^  xoxo  zozo 

while  from  Equations  (49)  and  (50)  we  have  a  =  d  V  /(v  t  -  v  ) 

J  z  z  ozRL  v  ozPL  ozRL' 

since  f(T)  «  a  T  for  thin  targets,  while  c  =  -  d  T/(v  -  v  „,) 
z  z  z  ozPL  ozRL 

in  any  case.  Similarly  we  can  relate  c  to  d  and  a  to  d  since 

xx  xx 

f(L„)  «  a..  T  tan  0_  for  thin  targets  from  Equations  (55)  and  (56).  We 

A  A  U 

can  also  find  cross  constraints  between  the  x  and  z  parameters  by  equat¬ 
ing  Equation  (67)  (or  Equation  (69))  to  its  x  analog  since  time 
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(like  mass)  is  common  to  all  components.  For  example,  at  the  time  of 

perforation,  x  =  v  sin  0  ,  z  =  v  cos  0  and  z  =  T.  Since  we  know 
PP  P  PP  P  P 
v  ,  0^  and  T  experimentally  and  can  estimate  x^  experimentally,  we  can 

obtain  cross  constraints  for  as  many  values  of  vq  (and  0  )  as  are  needed. 


The  same  can  be  done  for  ricochets  with  x  =  vR  sin  0R, 


z  =  vR  cos  0R, 


z  =  0  and  x  estimated.  In  addition  (at  least  for  some  projectile-target 
combinations)  other  conditions  must  be  satisfied  in  the  ricochet  region. 


For  example,  from  Figure  2  it  is  clear  that  vR  =  ^ xR2  +  zR2  goes  through 

a  maximum  as  a  function  of  v  .  Thus  at  some  value  of  v  ,  dv  /dv  -  0 
2  °  o  k  o 

and  v  /dv  1  <  0,  which  imposes  additional  cross  contraints.  Simi- 
R  o 

larly,  from  Figure  3  (and  similar  data)  it  seems  likely  that  0R  =  arc 
tan  (x  /z  )  goes  through  a  minimum,  so  that  perhaps  d  0  /dv  =  0  and 

~  K  K  r\  u 

d  0„/dv_2  >  0  at  some  value  of  v_ . 

K  O  u 


For  the  zero  obliquity  case  shown  in  Figure  1  there  is  no  x-component. 
There  is  also  no  ricochet  limit  shown  although  one  probably  exists. 

Table  I  contains  estimates  made  from  Figures  11  to  17  of  Backman  and 
Finnegan  with  estimates  of  error  ±  . 05mm/ps  on  the  speeds.  For  this 


Table  I.  Limit  Speeds  vs  0Q  (1  gm  Steel  Spheres  vs  6061  ML  plate) 


0 

0° 

O 

O 

45° 

60° 

o 

cos  0 

l 

.94 

.71 

.5 

0 

VoRL  (mm/ys] 

? 

.5 

.71 

1.06 

vozRL(mm/ys) 

? 

.47 

.5 

.53 

voPL 

.71 

.73 

.9 

1.2 

.71 

.69 

.65 

.6 

\r  -v 

? 

.22 

.15 

.07 

ozPL  ozRL 

kind  of  error  we  can  say  that  vqzRL  is  approximately  .5mm/ps  over  the 
whole  range,  including  perhaps  0q  =  0°.  However,  vozpL  clearly  decreases 

If  we  take  a  =  [v 

7 


as  0  increases  (T  =  constant) 
o 


ozRL^VozPL  VozRL^^z 


=  [ . 50/ ( .71  -  .50)J  dz  =  2.38  dz  and  cz  =  -  dz  T/(vQzpL  -  vqzRL) 

=  -  d  (9 . 53) / ( . 71  -  .50)  =  -  45.38  d7  for  0^  =  0°  with 

Z  Z  O 

0  <  d  <  .25  for  m.  =  1  gm,  then  we  only  have  to  estimate  b?  in 
z  z  u  “ 
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order  to  constrain  d  and  so  a  and  c  .  Gudunov  and  co-workers17  have 

z  z  z 

measured  the  kinematic  coefficient  of  viscosity,  v,  for  copper,  aluminum, 
lead  and  steel.  They  give  v  =  1.0  m7/s  for  aluminum,  2.5  rrr/s  for 
duraluminum  and  5.5  m^/s  for  steel.  The  dynamic  coefficient  of  viscosity 
is  p  =  pv.  For  example,  for  aluminum  with  p  =  2.7  gm/cm^,  p  =  2.7  x  10^ 
gm/(cm-sec).  If  a  Stokes  type  of  law  were  to  hold  for  a  non-deforming  sphere 

penetrating  aluminum  then  b  =  6-rrRy  =  6tt  Cm)  (2.7  x  IQ4  gm/(cm-sec)) 

6  ^  ^  ^ 

=  . 16  x  10  gm/sec  =  .16  gm/ps  or  2.5  times  this  for  duraluminum.  In 

this  way  we  can  estimate  bz  from  an  experiment  of  a  different  type.  If 
we  obtain  reasonable  agreement  with  penetration  experiments,  then  we  will 
be  encouraged  to  believe  that  the  constant  bz  is  more  than  just  a  fitting 
parameter.  If  we  take  b2  =  .10  gm/ps  and  dz  =  .0023  gm/(ps)  ,  we  obtain 
the  solid  line  in  Figure  I  calculated  for  v0  increments  of  0.1  km/s  (or 
0.1  mm/ps).  The  X’s  are  the  data  points  of  Backman  and  Finnegan  (their 

Figure  11).  The  model  allows  a  ricochet  region  for  0  <  vQ  <  .50  mm/ps 

with  an  embedment  region  for  .50  <  vQ  ^  .71  mm/ps.  In  the  perforation 
region  there  is  rough  agreement  with  experiment  but  the  curvature  of  the 
calculated  line  is  less  then  that  of  a  line  fitted  through  the  experi¬ 
mental  points.  This  is  probably  due  to  our  neglect  of  the  force  term 
t>2z  vz^*  N°  attempt  was  made  to  optimize  the  parameters.  Since  0O  =  0, 
vox  ~  0  and  Lx  =  T  tan  0O  =  0  so  5X  =  x  =  x  =  0  for  all  time  for  any  x- 

parameters.  The  x-parameter  values  in  Figure  1  are  dummy  values  to 

satisfy  the  program  input  requirements  of  the  9830A  Hewlitt-Packard 
Computer  (see  Appendix) . 

In  Figures  2  and  3  the  data  of  Backman  and  Finnegan  (their  Figures 
14  and  15)  for  the  same  parameters  as  in  Figure  1  except  0Q  =  45°  is 
shown  by  X's.  The  calculated  solid  curves  were  obtained  by  using  a  = 

.002  gm/(ps)2,  bx  =  .3  gm/ps,  cx  =  -  .4  gm/ps,  dx  =  .0225  gm/(ps)"  with 
bz  and  d7  the  same  as  for  0^  =  0°,  but  with  a^  and  c^  determined  using 
vozPL  =  *65  km/s  instead  of  .71  km/s  from  Tabfe  I.  fto  attempt  has  been 
made  to  optimize  the  agreement  which  presumably  could  be  improved  with 
some  effort.  It  might  be  wiser  to  apply  this  theory  to  a  greater  variety 
of  data  first  to  learn  more  about  the  trends  shown  by  the  parameters. 

This  will  be  done  in  future  work  before  we  return  to  optimizing  data 
fits. 

From  Table  I  we  noted  that  v  =  v  cos  0  =  .5  was  essentially 

ozRL  oRL  o  J 

constant,  implying  that  v  T  =  .5/cos  0  .  We  saw  however  that  v  - 

oRL  o  ozPL 

VoPL  COS  9o  decreased  as  0Q  increased,  implying  that  is  some 

other  function  of  9  .  If  we  postulate  that  v  n.  and  v  riT  differ  by  a 

o  r  oPL  oRL  7 

constant  amount,  namely,  v  n  =  . 21  +  v  OT  =  . 21  +  .5/cos  6  so  that 

oPL  oRL  o 

VozPL  =  +  -21  cos  0^,  we  can  reproduce  the  values  given  in  Table  I. 

The  parameter  a  now  becomes  a  =  d  v  ^T/(v  -  v  _ )  =  d 

z  z  z  ozRL  ozPL  ozRL  z 

(2.38/cos  0  ) .  In  our  calculation  for  Figures  2  and  3,  we  assumed  that 
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d  had  the  same  value  for  0q  =  45°  as  it  did  for  90  ~  0  which  made  az 

increase  as  0  increased.  If  az  were  to  be  independent  of  0Q  we  would 

have  to  assume  dz  proportional  to  cos  0Q.  Actually,  az  is  piupun-i^.a, 

to  the  z-components  of  the  projectile  presented  area  and  the  contact 
area.  Consequently  az  ought  to  be  proportional  to  cos  0Q  instead 

(cos  0  )~*  as  assumed.  This  implies  that  dz  ought  to  be  proportiona 
to  (cos  0  )2.  A  similar  argument  applies  to  cz  which  is  also  propor¬ 
tional  toV  In  future  work  we  will  explore  the  merits  of  including 
the  dependence  of  these  parameters  on  obliquity. 

If  e  is  imaginary,  than  the  motion  has  a  periodic  form  with  decay- 

Slau™.^:  executed  before  the  motion  ceasedwhen t  grew 


arbitrarily  large.  For  imaginary  3z  ±  %  -1-  f  ~z'  "o  z 

H  >  m  a  2  >  0,  for  example,  we  have 
~z  o  z 


ri  ,, 


i  6  ) 
z  z 


)  t  -  6  cos  0)  t 
z  c  z  J 


+  5^ 


-a  t  ' 

7  -  e  Z  v  cos  u> 
z  e  oz 


,t  -  f  f,  voz  -  5z  (“z2  *  “z2))  Sin  V 


with  analogous  relations  for  the  x-component.  As  before  from  Equations 
(70)  and  (71)  we  can  find 


e  sin  a)zt  =  [<5Z  z  +  vqz  (z  -  «Z)1/X 


e  Z  cos  u2t  =  ~  t(voz  az  V  z 


(z  -  6  )  (az  VQZ  -  6z  (az2  +  «/)>]/) 
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where 


X  =  —  [v  2  +  6  2  (a  2  +  u  2)]  =  —  [v  2  +  5  2  (— Z-)] 

0)  OZ  z  Z  Z  J  J  /i\  L  *7  '  j 

Z 


0)  OZ  z  TT1 

Z  O 


(74) 


or  by  dividing  Equation  (72)  by  Equation  (73) 


tan  u^t  =  u»z  [6z  z  +  vqz  (z  -  6  )]/[v_  -  a,  6J  z 


OZ  z  z 


+  (z  -  6_)  {a  V  -  6  fa  2  +  CJ  2)  11 

-  Z  OZ  Z  '  Z  z  '  '  J 


(75) 


The  times  at  which  turning  points  of  the  motion  occur  can  be  found 
by  setting  z  =  0  in  Equation  (75)  with  the  successive  values  of  w  t 
occurring  modulo  tt .  Jhe  perforation  limit  can  be  represented  by  the 

condition  that  the  first  maximum  of  z  is  equal  to  the  target  thickness 
(zmax  i  =  T)  while  the  ricochet  limit  can  be  represented  by  z^x  ^  <T, 

z  .  i  =  0.  Since  oscillatory  embedment  is  a  motion  allowed  by  the 
mm  1  J  ^ 

general  model,  provided  d  >  m  a  ,  z  <  T  and  z  .  >  0,  it  might 

r  z  o  z  max  1  mm  1  ’  b 

be  interesting  to  design  an  experiment  to  look  for  this  event  since  there 
does  not  seem  to  be  any  information  on  the  subject  in  the  published 
literature.  However,  such  an  experiment  would  require  a  motionless 
target  (very  massive)  since  this  model  takes  no  account  of  target  motion 
(a  reasonable  approximation  for  many  targets  during  the  short  times  in¬ 
volved  for  perforations  or  ricochets).  In  this  paper  we  will  not 
further  explore  the  possibility  that  one  or  more  3S  (s  =  x,  y,  z)  may 
be  imaginary.  Of  course  mixed  6  (one  component  6  real  or  zero,  another 
imaginary)  are  allowed  in  general. 

In  summary,  we  have  shown  how  to  apply  the  solutions  of  Equations 
(43)  and  (44)  to  a  particular  case  of  plate  target  penetration  for  which 
experimental  exit  speeds  and  angles  were  available  in  both  ricochet  and 
perforation  regimes.  The  agreement  between  theory  and  experiment  can 
certainly  be  improved,  but  the  overall  agreement  shows  promise.  An 
important  improvement  in  agreement  is  expected  when  force  terms  involving 
the  square  of  velocity  components  are  included  (as  they  will  be  in  future 
work) .  This  should  take  the  burden  of  agreement  in  the  perforation  region 
off  the  first  power  terms  and  allow  them  to  aid  the  constant  and  penetra¬ 
tion  depth  terms  in  improving  agreement  in  the  ricochet  regime.  The 
reader  is  urged  to  use  the  computer  code  in  the  Appendix  with  somewhat 
different  parameter  values  in  order  to  get  a  feel  for  how  this  might  be 
possible.  We  have  also  noted  the  considerable  variety  of  trajectories 
which  are  allowed  by  the  theory,  many  of  which  might  be  realized  in 
practice  with  a  suitable  choice  of  materials. 
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VI.  MASS  LOSS,  EXPERIMENT  AND  EXISTING  THEORY 


A.  Comments  on  the  Experimental  Data  Base 

The  largest  data  set  which  relates  to  the  question  of  mass  loss  by 
a  penetrator  consists  of  the  Thor  collection3  for  steel  fragments  to¬ 
gether  with  its  extension  to  tungsten  fragments4.  In  the  Thor  collec¬ 
tion  residual  mass  was  frequently  not  recorded,  but  when  it  was,  it  was 
never  equal  to  the  striking  mass.  In  the  extension  to  tungsten  fragments, 
the  residual  mass  after  perforation  was  almost  always  recorded  and  was 
always  less  than  the  striking  mass.  In  other  words,  if  our  aim  is  to 
model  at  least  the  important  features  of  penetration  our  model  must  be 
amended  to  include  non-constant  mass. 

"With  low  striking  velocities,  the  loss  in  weight  of  a  fragment 
during  perforation  is  small  and  is  usually  ignored.  In  such  cases,  the 
residual  velocity,  alone,  serves  as  a  good  measure  of  the  resistance  of 
the  target  to  perforation  and  the  capacity  of  the  residual  fragment  for 
perforating  another  target.  As  striking  velocity  increases,  the  break¬ 
up  of  the  fragment  becomes  more  and  more  pronounced  until  this  aspect 
of  the  impact  must  be  taken  into  account.  The  residual  weight  of  the 
fragment  must  be  known  as  well  as  the  residual  velocity  before  any 
reasonable  estimate  can  be  made  of  the  capacity  of  the  fragment  for 
perforating  another  target."3 

There  are  other  data  collections  concerning  rod  penetrators  and  these 
also  record  mass  loss  as  a  frequent  phenomenon  of  practical  interest31 . 

In  addition  many  individual  reports  record  mass  loss. 

There  appear  to  be  two  basic  types  of  mass  "loss".  The  first  can  be 
described  as  erosion  or  wearing  away  of  the  projectile  occuring  mostly 
at  the  tip  or  leading  edge.  This  might  occur  with  or  without  prior 
distortion  of  the  tip.  The  second  can  be  described  as  fracture  or 
breakup  of  the  projectile  into  pieces  of  more  or  less  comparable  size, 
differing,  say,  by  no  more  than  an  order  of  magnitude.  This  phenomenon 
can  be  called  a  mass  "loss"  only  if  one  chooses  to  disregard  the  smaller 
pieces  of  the  broken  projectile  and  concentrates  on  the  largest  surviv¬ 
ing  piece.  This  point  of  view  was  adopted  by  those  who  collected  the 
Thor  and  extended  Thor  data  sets.  An  attempt  was  made  to  justify  this 
procedure: 

"In  many  experimental  cases,  the  weight  of  the  largest  piece  of 
residual  fragment  approximates  the  total  weight  of  fragment  perforating 
the  barrier  target.  At  any  rate,  the  capacity  of  a  fragment  to  perforate 
a  primary  target  beyond  an  initial  barrier  can  be  conservatively  esti¬ 
mated  by  considering  only  the  largest  piece  of  fragment  which  perforates 


2V.  P.  Lambert ,  "The  Terminal  Ballistics  of  Certain  65  Gram  Long  Rod 
Penetrators  Impacting  Steel  Armor  Plate,"  ARBRL-TR-02072,  May  1978. 
(AD  M0577S7) 
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the  barrier.  This  approach  is  justified  whenever  the  hypothetical 
primary  target  is  one  for  which  damage  from  the  impact  of  small,  slow- 
moving  particles  is  not  anticipated,  i.e.,  damage  to  such  a  target  will 
essentially  be  that  caused  by  the  largest,  fastest  particle  that  impacts 
on  it."3 

The  authors  admit  however  the  desirability  of  observing  more  than 
the  largest  piece  so  that  one  can  estimate  damage  to  all  types  of 
interior  components.  However,  they  seem  to  ignore  the  possibility  that 
the  projectile  might  breakup  into  pieces  of  comparable  size.  Suppose, 
for  example,  the  largest  piece  were  0.4  mQ,  while  three  other  pieces 
half  this  size  but  with  comparable  residual  speed  accompanied  the  largest 
piece  in  a  debris  cone.  Then  the  target  damage  might  be  much  worse  than 
that  estimated  by  considering  only  the  largest  piece.  Add  to  this  the 
target  plug,  spall  and  other  debris  and  the  practical  need  for  recording 
more  than  the  largest  piece  becomes  clear.  From  a  theoretical  point  of 
view  the  need  for  recording  the  number  and  sizes  of  as  many  such  pieces 
as  can  be  practically  measured  is  also  clear.  Unfortunately  this  was 
not  done  even  in  the  extension  of  the  Thor  data  to  tungsten  fragments 
which  used  a  similar  wallboard  collection  method  and  recorded  only  the 
largest  piece.  However,  x-rays  were  used  to  measure  residual  speed  in 
some  of  these  later  experiments  so  there  may  be  some  unreported  fragment 
mass  distribution  information  still  recoverable. 

If  one  is  to  compare  a  physical  theory  with  experiment  and  not  merely 

fit  an  arbitrary  function,  at  least  three  things  must  be  known  about 

residual  penetrator  masses,  (1)  we  must  know  whether  erosion  and/or 
breakup  has  occurred;  (2)  if  erosion  has  occurred,  we  must  know  its 
extent  and  if  breakup  has  occurred  we  must  know  the  masses  of  all  the 
pieces;  and  (3)  we  must  know  whether  erosion  and/or  breakup  has  been  cans 
by  the  target  or  by  the  observation  technique.  We  are  interested  in 

events  caused  by  the  target.  If  erosion  and/or  breakup  is  caused  by  use 

of  a  recovery  medium,  this  leads  to  falsification  of  the  data  we  need. 

One  might  imagine  firing  into  the  recovery  medium  alone  in  order  to 
correct  for  this  falsification.  Such  a  procedure  might  seem  to  resemble 
the  calibration  procedure  which  is  usually  carried  out  when  such  a  medium 
is  used  to  determine  residual  speeds.  However,  this  is  not  satisfactory 
for  mass  determinations.  For  example,  it  is  entirely  possible  that 
incipient  cracks  produced  by  the  target/projectile  encounter  will  be 
converted  into  projectile  breakup  by  the  recovery  medium,  although  the 
projectile  emerging  from  the  target  was  intact.  To  show  that  this  did 
not  happen  would  require  a  pre-conditioning  of  the  projectile  comparable 
to  firing  it  through  a  target.  For  observation  of  residual  masses  direct 
methods  such  as  x-rays  seem  preferable  in  spite  of  their  difficulties. 
Some  experiments  of  this  type  have  been  performed22.  However,  the  Thor 
data  collection  and  its  extension  do  not  contain  enough  information  for 


Herr  and  C.  Grabarek >  "Ballistic  Performance  and  Beyond  Armor  Data 
for  Rods  Impacting  Steel  Armor  Plates BRL  MR  2575 ,  1976 . 

(AD  #B009979L) 
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comparison  with  a  physical  theory  in  spite  of  the  large  number  of  firings 
involved.  In  general,  they  do  not  tell  us  the  three  pieces  of  informa¬ 
tion  enumerated  above.  In  addition,  they  contain  very  few  cases  in 
which  only  one  variable  was  changed  while  all  others  were  held  constant, 
something  which  is  essential  if  fundamental  comparisons  are  to  be  made. 

B.  Comments  on  Previous  Models 

Some  of  the  advantages  and  disadvantages  of  purely  empirical  models 
such  as  the  second  Thor  equation  have  been  pointed  out  and  will  be 
elaborated  upon  below.  Mass  erosion  and  breakup  is  now  being  incorporated 
into  some  advanced  continuum  models  with  of  course  an  increase  in  the 
complexities  which  already  exist  for  such  computer-dependent  models. 

In  the  class  of  theoretical  but  simple  intermediate  models  there  are  at 
least  three  which  have  been  proposed  to  describe  mass  erosion,  although 
there  seem  to  be  none  for  breakup. 

The  earliest  intermediate  erosion  model  was  given  by  Tate  .  He 
adapted  a  modified  hydrodynamic  theory  first  proposed  by  Hill,  Mott  and 
Pack2*4,  later  published  by  Bishop,  Hill  and  Mott^3  and  still  later 
modified  by  Eichelberger26.  If  we  assume  that  the  deceleration  is 
negligible  during  some  part  of  the  penetration  then  Bernoulli's  steady- 
state  equation  will  hold  approximately  for  the  pressure  at  the  stagna¬ 
tion  point 

P  =  -jPtu^  +  P  =  jP  (v  "  +  Y  (76) 

where  R  and  Y  are  strength  parameters  or  stresses  above  which  the  target 
and  projectile  respectively  behave  as  fluids.  Here  u  is  the  penetration 
speed  and  v  is  the  speed  of  the  rear  of  the  penetrator  for  zero  obliquity 
impact  on  a  target  plate,  while  p^.  is  the  density  of  the  target  and  p  is 
the  density  of  the  projectile.  Since  u  <  v,  we  have 


u  = 


(77) 


where  y  = 


and  A  =  —  (1 


y2)  (R  -  Y). 


Since  v  and  u  are  the 


23A.  Tate ,  "A  Theory  for  the  Deceleration  of  Long  Hods  After  Impact ," 

J.  Mech.  Phys.  Solids  35,  387 ,  1967,  and  "Further  Results  in  The  Theory 
of  Long  Rod  Penetration,"  Ibid.  If,  141,  1969. 

24 Unpublished  Ministry  of  Supply  Report,  1945. 

2SR.  F.  Bishop,  R.  Hill  and  N.  F.  Mott,  Proa  Phys  Soc  5f,  147,  1945. 

26i?.  J.  Eichelberger ,  J.  Appl.  Phys.  27,  68,  1956. 
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speeds  of  the  back  and  front  end  of  the  rod  respectively,  then  the 
decrease  in  length  with  time  is 

l  =  -  (v  -  u)  .  (78) 

The  equation  of  motion  of  the  back  end  of  the  rod  was  taken  to  be 

P  l  v  =  -  Y,  (79) 

assuming  only  a  constant  deceleration  supported  by  the  strength  of  the 
rod.  Equation  (79)^  not  strictly  speaking  compatible  with  Equation 
(76)  which  assumes  v  =  0.  It  is  also  for  the  speed  of  the  rear  of  the 
rod,  not  its  center  of  mass.  Since  dt  =  -  p  l  dv/Y  from  Equation  (79), 
we  can  eliminate  dt  from  Equation  (72)  and  write 


d£ 

l 


=  (v  -  u)  dv. 


(80) 


If  we  use  Equation  (77)  in  Equation  (80)  to  eliminate  u,  and  integrate 
from  initial  values  of  l  =  L  and  v  =  V  we  obtain 


+  A  - 


(81) 


Here  l  cannot  vanish  for  finite  V.  Tate  distinguished  several  cases. 

If  R  >  Y,  the  rod  behaves  as  a  fluid  until  penetration  ceases  and  u  =  0 

in  Equation  (77)  or  v  =  u  .  The  length  of  the  rod  at  this 

c  f  l-u“ 

point,  namely  is  found  by  using  v in  Equation  (81).  Beyond  this 

point  the  hydrodynamic  theory  does  not  apply.  Instead  the  bulk  of  the 

rod  is  taken  to  be  a  rigid  body,  and  v  in  Equations  (78)  and  (79)  is 

taken  to  be  the  speed  of  the  center  of  mass  with,  of  course,  u  =  0. 

The  center  of  mass  continues  to  move  forward  but  the  tip  merely  erodes 

in  place.  Now  dt  =  -  d£/v  from  Equation  (78)  can  be  eliminated  in 

Equation  (79)  to  obtain  Equation  (80)  with  u  =  0.  If  we  integrate  from 

initial  values  SL  =  £  and  v  =  v  to  final  values  and  v  =  0  we  obtain 

c  c  f 

the  final  length 
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£f  =  ic  exp  [-  (R-Y)/Y]  .  (82) 

If  the  impact  speed  V  <  vc>  the  tip  merely  deforms  without  penetra¬ 
tion  (u  =  0)  and  an  integration  of  Equation  (80)  from  initial  values 
£  =  L  and  v  =  V  to  £  =  £^  and  v  =  0  gives 

If  =  L  exp  [-  (y  pV2)/Y]  (83) 

for  the  final  length  of  the  undeformed  portion  of  the  rod  which  depends 
exponentially  on  the  ratio  of  the  initial  kinetic  energy  density 

j  0\‘L  to  the  rod  strength  Y.  The  lack  of  any  penetration  is  an  un¬ 
reasonable  feature  and  holds  only  for  weak  rods  versus  strong  targets. 
Tate  points  out  that  his  model  only  applies  to  perfectly  plastic  rod 
material  for  which  the  plastic  wave  speed  C  =  0.  He  compares  Equation 
(83)  to  Taylor's27  result  for  the  deformation  of  a  rod  striking  a  rigid 
target,  namely, 


£f  =  L  exp  [-  (^r)  v  (V  +  2C)  ]  (84) 

which  reduces  to  Equation  (83)  when  C  =  0. 

If  R  <  Y,  then  the  rod  behaves  as  a  fluid  until  u  =  =  V-A/(l-u2) 

=  /2  (Y-R) /p^  from  Equation  (77).  The  length  of  the  rod  at  this  point, 

£  ,  is  found^by  using  in  Equation  (81).  Beyond  this  point  the  rod 
c  c 

continues  to  penetrate  as  a  rigid  body  with  Since  the  length 

does  not  change,  l  =  0  in  Equation  (78)  and  u  =  v  throughout  this  second 
stage  with  v  interpreted  as  the  speed  of  the  center  of  mass  or  of  either 

end  and  Y  =  y  pt  v2  +  R  from  Equation  (76).  Then  Equation  (77)  becomes 

p  £  v  =  -  (i  p  v2  +  R)  (85) 

which  is  simply  Equation  (15)  with  b  =  0,  mQ  =  p  A,  a  =  AR  and 

i  • 

c  =  —  p  A  with  A  the  cross-sectional  area.  If  we  take  v  =  z  and 
,  2  t 

v  =  v  dv/dz  in  Equation  (85)  and  integrate  from  z  =  0  and  v  =  vp  to  final 

values  z  =  T  and  v  =  0  we  obtain  Robertson's  result.  Equation  (16),  with 
=  0,  namely 

27G.  I.  Taylor ,  Proa.  Boy.  Soo.  287,  1948. 


55 


TAc  =  -J  in  [1  +  Pt  vc2)/R]  =  ~  Zn  fl  +  p2  (Ip  v  2)/R]  (86) 

where  T  is  now  the  extra  depth  of  penetration  in  a  semi-infinite  target 

and  depends  logarithmically  on  the  ratio  of  the  rod's  energy  density  at 

the  beginning  of  its  rigid  body  motion  to  the  strength  of  the  target. 

2  2 
If  we  use  v  =  2  (Y-R)/p  in  Equation  (86)  we  obtain  T  =  ( £  / p  )  C,n  (Y/R). 
^  U  c 

If  R  =  Y,  then  A  =  0  and  u  =  v/(l  +  p)  in  Equation  (77)  while 

Equation  (80)  becomes  dZ/Z  =  (p/Y)  (|~— -)  v  dv  so 


(Sl/L)  =  exp  [-  W2  -  v2)]  (87) 


and  the  final  length  Z  -  Z^  when  v  =  0  depends  exponentially  on  (y  p  V^)/Y 
and  as  before  can  never  be  zero  for  finite  V.  The  depth  of  penetration 

z  =  /tu  dT  =  T— j; —  /tv  dt.  If  we  use  Equation  (79)  to  eliminate 
o  1  M  o 

dt  =  -  p  H  dv/Y  and  use  Equation  (87)  for  l  we  obtain 


z  =  (L/p)  1  -  exp  2  TTl'  '+  V)  (X'2  *  y2)}  •  (88) 

—  -J 


If  p^  =  p  so  p  =  1,  then  A  =  0  and  the  above  formulas  simplify 
further,  although  closed  form  solutions  for  z  (v)  cannot  be  found  for 
the  general  case  R  ^  Y  unless  R  =  nY  with  n  integer. 

A  model  similar  to  Tate*s  has  been  proposed  by  Walters  and  Majerus28 
These  authors  give  a  unified  theory  for  zero  obliquity  penetration  by 
shaped  charge  jets  or  kinetic  energy  projectiles.  In  particular,  in 
their  Appendix  II  for  steady-state  jet  penetration  neglecting  viscosity 
they  obtain  in  their  Equation  ( 1 1 - 1 j 


[a 


(v  -  uH  -  [Qj  + 


(89) 


28fv\  P.  Walters  and  J.  N.  Majerus ,  rr Impact  Models  for  Penetration  and 
Hole  Growth 3  "  ARBRL-7R- 02069 ,  May  1978.  (AD  UA0S6294) 
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which  is  the  same  as  Equation  (76)  above  provided  a  =  2Y  and  o  =  2R 

and  the  jet  is  not  stretching  so  that  the  speed  of  its  rear  is  equal  to 
the  speed  of  the  interface  of  its  tip  with  a  postulated  interaction 
region.  In  the  special  case  R  =  Y  above  we  found  u  =  v/(l  +  u) .  If  zq 

is  the  standoff  distance  and  we  use  the  form  v  =  (z  +  z)/t,  then 

o 

u  =  dz/dt  =  v/(l  +  n)  =  (1/(1  +  p))  (z  +  z)/t  so 

1 

z  =  zo  [(t/tQ)1+P  -  1]  (90) 

which  is  the  formula  obtained  by  Allison  and  Vitali29  and  discussed  by 
DiPersio,  Simon  and  Merendino30. 

In  their # Appendix  IIIA,  Walters  and  Majerus  consider  penetration  by 
a  rigid  rod  i  =  0,  u  =  v  (the  second  stage  of  the  case  R  <  Y  considered 
above)  and  their  Equation  ( 1 1 1 - 1)  with  the  height  of  their  interaction 
region  H  =  1  becomes 

2  p  £  v  =  -  (p.  v  +  o  )  (91) 

c  t  t 

which  is  Equation  (85)  provided  a  =  2  R.  In  their  Appendix  III-B  they 

consider  penetration  by  a  non-stretching  jet  (or  rod)  with  constant  speed 
v  and  obtain  as  their  Equation  (1 1 1-9)  for  u 


u  =  a'  +  b'  u  +  c'  u^  (92) 

'l  J-J 

with  pH  a"  =  (a  -  at)  +  p  v*\  pH  b'  =  -  (r^  p'  +  2  p  v)  and  p  He'  = 

(p  -  Pt)  where  A  =  ir  rQ  is  the  cross  sectional  area  of  the  cylindrical 

rod  and  p'  is  the  target  viscosity.  Since  p  H  A  is  the  mass  in  the 

interaction  region  moving  with  velocity  u,  Equation  (92)  is  the  equation 
of  motion  for  this  region  or  for  the  entire  rod  if  H  =  £.  The  coefficient 
b'  is  negative,  but  the  coefficient  c'  is  negative  only  if  p^  >  p,  while 
the  coefficient  a'  is  negative  only  if  >  o  +  p  v  .  Equation  (92)  has 

the  same  form  as  Equation  (15)  with  a'  =  -  a/mQ>  b'=  -  b/mQ  and  c'=  -c/mQ 


29F.  E.  Allison  and  R.  Vitali,  "A  New  Method  of  Computing  Penetration 
Variables  for  Shaped  Charge  Jets,"  BRL  1184,  1963.  (AD  #400483) 

30I?.  DiPersio,  J.  Simon  and  A.  Merendino,  "Penetration  of  Shaped  Charge 
Jets  Into  Metallic  Targets,"  BRL  Report  1296,  1965.  (AD  #476717) 
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(with  u  playing  the  role  of  v)  and  its  solutions  are  formally  the  same 

so  we  can  find  u  (t)  and  J u  u  dt  in  closed  form  for  any  value  of  q  = 

2 

4  a'  c'  -  fb")  in  Equation  (22),  Then  from  Equation  (78)  we  can  find 


£ 


L 


v  t 


+ 


JL  u  dt 


(93) 


in  closed  form.  For  example,  for  q  =  0,  the  simplest  case,  we  have  by 
analogy  with  Equation  (26) 


o 


(94) 


which  can  be  used  in  Equation  (93)  to  obtain  £  (t).  Similarly  equations 
analogous  to  Equations  (24)  and  (28)  could  be  used  for  q  <  0  or  q  >  0 
respectively.  Since  v  is  not  really  a  constant  in  practical  cases, 
especially  toward  the  end  of  the  motion  in  a  semi -inf inite  target,  we 
can  expect  Equation  (93)  to  agree  with  experiment  only  approximately  and 
only  during  the  early  part  of  the  motion.  Equation  (93)  requires  £  to 
vanish  in  finite  time  for  any  q,  contrary  to  experiments  with  semi-infinite 
targets,  but  the  agreement  is  reasonable  during  the  initial  stages  as 
noted  by  Walters  and  Majerus.  Presumably  numerical  solutions  of  their 
more  general  equations  will  give  predictions  in  better  agreement  with 
experiment  at  later  times. 


In  their  Appendix  IV  Walters  and  Majerus  consider  the  case  of  a  non¬ 
eroding  rod  (£  =  0)  for  which  they  obtain  Equation  (IS)  for  u  =  v  with 
a  =  A  a  ,  b  =  2tt  H  y  and  c  =  A  p  where  a  ,  u  and  p  are  the  yield 

1  L  L  t  t  t 

stress,  viscosity  and  density  respectively  for  the  target.  The  solutions 
have  already  been  discussed  above  and  by  Zook2. 

Another  erosion  model  for  zero  obliquity  impact  has  been  proposed  by 
Recht31.  This  model  was  previously  described  at  the  Fourteenth  Annual 
Meeting  of  the  Society  of  Engineering  Sciences  at  Lehigh  University  in 
1977  and  earlier  at  the  Workshop  on  Mechanics  of  Impact  and  Penetration, 
December  1976,  Aberdeen  Proving  Ground,  Md.  An  outline  of  his  model 
can  also  be  found  on  pp  25-29  of  the  Backman  and  Goldsmith  review  article 
cited  above1.  When  the  projectile  speed  exceeds  an  assumed  constant 
plastic  wave  speed  C,  a  shockwave  develops  near  the  tip  of  the  projectile, 
eroding  it  away.  When  the  projectile  speed  falls  below  C,  deformation 
of  the  tip  takes  place,  shortening  it  because  of  lateral  spreading. 

When  this  spreading  exceeds  some  critical  amount  (postulated  to  be  25% 
in  Recht's  application  to  a  case  of  a  steel  cylinder  impacting  a  steel 

31i?.  F.  Recht,  " Taylor  Ballistic  Impact  Modeling  Applied  to  Deformation 
and  Mass  Loss  Determinations ,"  Int .  J.  of  Eng .  Sciences,  16,  809,  1978 . 
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plate),  then  the  material  which  has  spread  beyond  this  distance  is 
assumed  to  be  sheared  off  as  the  motion  proceeds. 


In  arriving  at  Equation  (41)  of  his  1978  paper,  Recht32  considered 
that  the  main  projectile  changes  its  mass  from  m  to  m  +  dm  and  changes 
the  speed  of  its  center  of  mass  from  v  to  v  +  dv  during  the  time  interval 
dt.  Here  dm  and  dv  are  negative.  The  eroded  mass  (-dm)  was  assumed  to 
have  the  plate  plug  speed,  v^.  More  properly,  at  the  end  of  the 


interval  dt  the  plug  speed  is  vp  +  dvp,  although  this  will  not  matter  as 

we  will  see  shortly.  If  one  makes  the  same  free  collision  assumption 
that  Recht  and  Ipson  made  in  their  1963  paper  cited  above5  with  mp 

already  formed  and  constant  throughout  the  motion,  then  the  momenta 
before  and  after  this  collision  would  be  equal,  that  is, 


mv  +  m^  vp  =  (m+dm)  (v+dv)  +  mp(vp+dvp)  +  (“dro)  (Vp+dVp)  .  (95) 

If  we  neglect  the  second  order  terms  dm  dv  and  dm  dvp,  we  find 

m  dv  =  (-dm)  (v-vp)  -  mp  dvp  (96) 

and  if  we  let  (F  dt)  =  m  dv  with  F  =  -°A,  we  obtain  Recht' s  Equation 
(41)  mentioned  above,  namely, 

a  A  dt  +  (-dm)  (v-vn)  =  mn  dvn  .  (97) 

More  properly,  when  dm  ^  0,  Newton's  law  requires 

m  dv  +  v  dm  =  F  dt  =  -  a  A  dt  (98) 


which  gives  us  from  Equation  (96) 

o  A  dt  +  v  dm  =  m  dv  (99) 

P  P  F 

instead  of  Equation  (97) .  There  seem  to  be  at  least  three  shortcomings 
in  Recht 's  model.  First  of  all,  he  makes  the  unrealistic  assumption  of 
a  free  collision  in  Equation  (95).  In  his  earlier  paper  with  Ipson  he 
was  able  to  compensate  for  this  by  equating  the  catchall  energy  E  to  the 
change  in  energy  of  the  mass  mQ  during  the  assumed  free  collision,  re¬ 
quiring  that  the  exit  speed  be  equal  to  the  speed  VJf  after  the  free 


3 ^Private  communication  to  J.  Dehn,  14  Dec  1978. 
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collision,  thus  producing  an  expression  for  the  shearing  work  W.  In  his 
erosion  model  above  he  does  not  employ  energy  balances  to  compensate  for 
the  unrealistic  assumption  of  a  free  collision,  so  the  error  might  be 
uncorrected.  However,  as  we  have  noted,  he  does  ignore  the  theoretical 
requirement  that  (F  dt)  =  m  dv+vdm,  so  that  this  error  might  compensate 
for  the  free  collision  assumption.  It  might  also  make  it  worse.  Thirdly, 
it  would  be  desirable  to  use  more  than  a  simple  constant  force  F  =  -aA. 
Even  his  earlier  model  with  Ipson  produced  a  result  similar  in  form  to 

1  9 

that  of  Robertson  who  used  F=  -  (a"  A)  -  (-  y  p  A)  v" .  We  have  also 

seen  the  desirability  of  including  terms  dependent  on  v  and  on  the  depth 
of  penetration,  s. 

Many  physicists  take  the  attitude  that  compensating  errors  are 
perfectly  acceptable  in  a  theory  as  long  as  they  result  in  a  definite 
prescription  for  calculations  which  agree  with  experiment.  There  is 
something  to  be  said  for  this  point  of  view  just  as  there  is  something 
to  be  said  for  purely  empirical  interpolation  models  which  also  agree 
with  experiment  at  least  for  limited  conditions.  However,  there  is  also 
something  to  be  said  for  the  view  that  we  should  not  be  satisfied  with 
a  fortuitous  compensation  of  errors  if  we  expect  to  make  any  progress  in 
achieving  a  broader  understanding  of  the  physical  world.  Dirac  is  one 
of  the  foremost  proponents  of  this  latter  view. 

"My  early  research  work,  in  the  early  1920 !s  was  based  on  Bohr  orbits, 
and  was  completely  unsuccessful.  I  was  taking  the  Bohr  orbits  as 
physically  real  and  trying  to  build  up  a  mathematics  for  them....  One 
sees  now  how  futile  such  work  was....  The  Bohr  orbits  were  an  unsound 
physical  concept  and  should  not  be  used  as  the  basis  for  a  theory.”33 

Dirac  goes  on  to  relate  how  his  dissatisfaction  with  the  Klein-Gordon 
equation  led  him  to  develop  his  own  relativistic  equation  for  the  elec¬ 
tron  from  which  the  concept  of  spin  naturally  followed.  He  also  de¬ 
scribes  more  recent  efforts  by  physicists  to  solve  the  relativistic  wave 
equation  for  interacting  particles,  efforts  which  have  led  to  divergent 
integrals  which  are  discarded  by  a  renormalization  of  fundamental 
physical  constants.  Although  most  physicists  are  satisfied  with  this 
since  their  calculations  agree  with  experiment,  Dirac  gives  cogent 
reasons  for  not  being  content.  He  then  says, 

”For  these  reasons  I  find  the  present  quantum  electrodynamics  quite 
unsatisfactory.  One  ought  not  to  be  complacent  about  its  faults.  The 
agreement  with  observation  is  presumably  a  coincidence,  just  like  the 
original  calculation  of  the  hydrogen  spectrum  with  Bohr  orbits.  Such 
coincidences  are  no  reason  for  turning  a  blind  eye  to  the  faults  of  a 
theory,”33 

33P.  A.  M.  Dirac,  " The  Mathematical  Foundations  of  Quantum  Theory,”  in 

Mathematical  Foundations  of  Quantum  Theory,  Ed .  by  A.  Ft.  Marlow,  N.Y. 

Academic  Press,  1978 . 
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VII.  A  NEW  MODEL  OF  MASS  EROSION 


A.  Development  of  the  Theory 

Here  erosion  is  being  thought  of  as  a  continuous  process  which  can 
be  described  by  differential  equations.  We  might  expect  at  least  the 
dissipative  terms  in  our  force  field  to  depend  on  the  instantaneous 
values  of  the  mass  (m)  and  /or  its  rate  of  change  (m) .  A  dissipative 
term  accounts  for  kinetic  energy  loss  by  a  decrease  in  velocity  in 
the  constant  mass  case.  When  the  mass  itself  is  decreasing,  additional 
energy  loss  should  also  be  described  by  such  a  term.  Let  us  presume 
that  this  is  true  and  proceed  to  make  particular  assumptions  about 
our  dissipative  terms. 

The  viscous  loss  coefficient  b|  has  the  dimensions  of  mass  per  unit 
time  so  a  simple  modification  of  this  term  might  be 

-(bls-alsm)s  (100) 

where  e  is  a  positive,  dimensionless  constant.  When  the  mass  change 

rate  is  negative  (m  <  0),  then  the  dissipative  coefficient  b^  is 

effectively  increased,  while  for  m  >  0  mass  gains  tend  to  offset 
velocity  losses  and  the  coefficient  b^g  is  effectively  decreased. 

When  m=0  we  have  the  constant  mass  caseand  return  to  the  original 
form  of  this  term.  Examples  involving  m<0  are  target  penetrators  or 
space  vehicle  re-entry  shields.  Examples  involving  m>0  might  be 
rain  or  sleet  in  the  atmosphere  or  crystallites  precipitating  from 
a  chemical  solution. 

2 

If  we  wish  to  modify  the  coefficient  of  the  v  term  we  might 
consider  a  form  like 


"  b2s  CVm)  *^01) 

This  gives  us  a  reduction  in  retardation  when  (mo/m)<l  and  an 

increase  for  (m  /m)>l  with  a  return  to  our  original  form  when 
o  2 

(m  /m)=l,  the  constant  mass  case.  We  will  not  treat  the  v  term  in 

this  paper. 

•  •  • 

Since  d(m  s)/dt  =  m  s  +  m  s  is  the  kinematic  expression  when 
m=m(t),  we  must  solve  the  equation 

*  2  2 

m  s  +  m  s  +  (bjs  -  m)  s  +  b2g  (mQ/m)  s'  +  ag  +  dls  s  +  d2s  s  =0  (102) 
where  m(t)  and  s(t)  (s=x,  y  or  z)  must  both  be  determined. 
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Let  *-*  (t)  be  an  arbitrary  transformation  of  the  independent  variable 
so  that  s  =  <j>  ds/d<f>  and  s  =  *  2  (d2s/d*  2)  +  }  (d^  /d<f>)  (ds/dcf))  .  Then 

Equation  (102)  becomes 


where 


B 


Is 


B2s  =  b2s  (mo/m)  ^ 


(103) 


(104) 

(105) 


(106) 


Now  let  us  choose  <}>  so  that  at  least  M  and  B,  are  constants  and 

Is 

explore  the  implications  of  this  choice.  One  implication  is  evident 
from  the  form  of  Equation  (103),  namely,  that  whatever  solutions  we 
found  for  z  (t)  in  Equation  (44)  can  be  used  to  write  s(*)  here. 
Another  implication  is  that  Equation  (104)  becomes 


M  -  mi})*'  moi0"  (104a) 

where  *q  is  the  value  of  *  at  t=0.  We  can  also  use  Equation  (104a) 
to  obtain 

dm/d*  =  i  [-2M/*3]  ^d*/d*j  = -2m  d*/d*  (107) 

and  use  this  to  eliminate  (md*/d*)  in  Equation  (105)  to  obtain 


Is 


■K-  "  ]  *  *  [»!,-  (Els‘  i)*0]  V 


(108) 


Now  for  the  constant  mass  case  (rn^o) ,  Bls=bls  *q  from  Equation  (108). 
However  in  this  case  we  require  B^b^,  so  *^1.  This  implies 
M=mo  in  Equation^ (1 04a).  If  we  put  Equation  (104a)  in  Equation  (106) 
we  find  B2s=b2s  ^ which  is  not  constant.  If  we  had  used  (m/m  ) 

instead  of  (m/m)  in  Equation  CL01) ,  then  B^b^,  also  a  constant. 
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teed  to  solve 


Consequently ,  Equation  (108,  is  the  second  relation  we  •»"- 

for  «t)  and  so  m(t)  and  s(t)  Since  m=mo/+  from  Equatton  (104a,  and 
;=.2  (i  /i2,  di/d*  from  Equation  (107),  Equation  (105,  becomes 


!)(■•«)  *  i"  * 

which  can  be  integrated  to  obtain 

*-[(1-bls/Bls)“P{-Bla  (+)/.[2”o(Els-i)]}*bls/'is]"1  tU0) 

where  we  have  used  iQ=l  and  let  *0  -  V  °  £“  '  EqUati°n 

(110)  can  be  integrated  once  more  to  obtain 


(bls/Bls)  41  +  [2  mo(£ls  2)  /Bls](1_  bls/Bls)[1  6XP  ^  ^ 


(111) 


[2mo(els"l)  /Blsj}J=t^  ipftsideQf  . 

which  gives  us  a  relation  beJ”®®^14>f^tion  of  4>  and  becomes  algebraic 
Equation  (111)  is  a  transcenaenu-  *  constant  mass  case 

only  in  certain  simple  cases  such  as  t>ls  »ls  t 

for  which  *=t)  or  bls=0  for  which 

/  •  \  .  It  .71  /I  m  \  t  I  (112) 


(%/2  mo)  ln  L1  +(mo/2  mo)  M 


nno.  The  general  form.  Equation  (111),  gives 
thentransformation°which  as  required  to  make  M  and  Bls_in  Equation 
(103,  constant.  By  using  Equation  (110,  in  Equation  U«1  « 
also  obtain 

m  =  m0((l-bls/Blsj  exp| -*/  (l  ”0  (Eis‘  2  )/Bls])  *  bls/Bls)  ( 

If  b  >  B,  ,  then  m  =  m  1  .  the  constant  mass  case.  If  bls  =  0.  tnen 
Is  IS  0 


*/mD  =  m0  [l  *  (m0/2  »0)  t]  2. 


(114) 
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In  summary,  instead  of  assuming  a  form  for  m(t),  we  have  made 
reasonable  assumptions  about  modifications  to  the  coefficients  in  our 
force  terms,  we  have  transformed  the  independent  variable  in  order  to 

obtain  another  equation  with  constant  coefficients  neglecting  the  v^ 
term  and  have  derived  a  relation  for  m  [<j>(t)J  •  This  enables  us  to 
write  closed  form  expressions  for  the  displacement  components 
s  [<!>(£)]  >  provided  we  require 


blx  /  ('lx  \)  =  bly  /(ely"!)=  blz  /(elz  '  \) 

and 


(115) 


B 


lx 


/ 


so  that  m 
Equations 


b,  =  B,  /  b. 
ix  ly  iy 


-hz' 


iz 


(116) 


and  4>  will  be  the  same  for  all  components  in  Equation  (108)  . 
(115)  and  (116)  then  guarantee  that 


Blx/(Elx-l)‘Bly/(ely-|)-Blz/(El2-?)  Cl  17) 

so  that  <J>  and  m  in  Equation  (111)  and  Equation  (113)  will  be  the  same 

in  all  three  component  Equations  (103)  neglecting  the  vz  term.  Such 
a  transformation  might  be  called  a  constant  mass  frame  transformation, 
since,  in  the  new  time  frame  <p,  M=m^  is  constant  and  is  acted  upon  by 

force  terms  containing  constant  coefficients  one  of  which  multiplies 
the  transformed  speed  ds/d<j>.  The  utility  of  such  a  transformation  can 
only  be  judged  by  comparison  with  experiment. 

First  let  us  examine  the  qualitative  behavior  of  our  solution  for 
m.  If  we  use  Equation  (111)  and  Equation  (113)  in  Equation  (107),  we 
obtain 


“  '  "o  exP  {  -♦  ’  [2  no  (els-  \)  /Bls]  )  018) 

where  we  have  also  used  Equation  (108).  We  see  that  Els=  \  corresponds 

to  the  constant  mass  case  since  it  makes  m  =  o  for  all  <j>  =  t  >  0 
(bjs  =  B|S)  in  Equation  (111)  .  Here  we  are  interested  in  the  case 

m<0  for  which  it  is  convenient  to  write  Equation  (108)  as  B  =  b  + 

1  •  Is  Is 

0ls-5>  (-V- 
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When  £js  >  -j,  then  B^s  >  b^s  >  0  and  the  right  side  of  Equation  (111) 

is  positive  and  increases  without  limit.  If  <|>  is  also  positive,  then 
both  terms  on  the  left  of  Equation  (111)  are  positive  with  the  second 
term  approaching  a  maximum  value  while  the  first  increases  without  limit 
so  that  for  late  times  <f>  becomes  approximately  proportional  to  t. 

Negative  <(>  is  not  allowed  since  both  terms  on  the  left  of  Equation  (111) 
would  then  be  negative.  For  <p  >  0,  Equation  (118)  says  that  |m|  decreases 
as  <p  (or  t)  increases.  Consequently  a  plot  of  m  versus  t  will  be  con¬ 
cave  upward  with  m  decreasing  to  the  limiting  value  mQ  S/Bj s) ^  as  t 

increases  without  limit.  In  this  case,  m  will  truly  vanish  only  when 


b.  =  0. 
Is 


When  e  <  j,  then  Blg  7  0.  If  Blg  >  0  because  bjs  >  |ej  -  j\  (-  mQ) , 

then  b,  >  B,  >0  and  Equation  (111)  has  a  solution  only  for  the  finite 
Is  1  c 


0  <  t  <  t 


MAX 


Is 

with  <j>  >  0. 


Negative  <j>  is  not  allowed  since  the  first  term 
on  the  left  of  Equation  (111)  would  be  negative  and  eventually  dominate 


the  second  term.  At  t^^,  dt/d<j> 


(m/mo) 


1/2 


from  Equations  (111), 


(110)  and  (104a)  and  vanishes  for 


mx 


[-  2  m  (e 


Is 


-  i^Bisi 


*n  [1/(1  -  Bls/bls)]  (119) 


or 


“MAX 


[2  mo  <Els 


-  - 


(bls/Bls’ 


In  [1/(1  -  Bls/bls)l).  (120) 


From  Equation  (118),  |m|  increases  as  4>  (or  t)  increases  so  that  a  plot 
of  m  versus  t  is  concave  downward  with  m  vanishing  at  t  =  t^^.  ^ 

B.  =0  because  b=  |e,  -  i|  (-  m  ) ,  then,  since 

Is  Is  Is  2'  o 


£im 


B 


Is 


->0 


=  4>/  [2  mQ/(e 


Is 


(121) 


Equation  (111)  reduces  to 


so  =  [rn  /(-m  )]  and  =  2  when  m  vanishes.  In  this  case, 

•  # 

from  Equation  (118),  m  =  mQ  and  m  decreases  linearly  with  time  as  follows 

m  =  mn  [1  -  (-  mn/mn)t]  .  (123) 

If  Bis  <  0  because  -  ~|  (-  mQ)  )  then  <b  must  again  be  positive 

and  a  solution  exists  only  for  a  finite  time  as  before  expressed  by 

Equation  (120).  Now  however  from  Equation  (118)  we  see  that  |m|  decreases 

as  <J>  (or  t)  increases  so  that  a  plot  of  m  versus  t  is  concave  upward  with  m 

vanishing  at  t  =  t..AV. 

MAX 

In  summary,  for  m  <  0,  we  can  obtain  two  types  of  mass  loss  with 
time,  depending  on  whether  is  greater  or  smaller  than  0.5.  If 

cls  >  0.5,  then  m  will  not  vanish  even  in  a  semi -infinite  target  provided 

^  0.  However,  for  sufficiently  small  b/b  it  might  appear  to 

vanish.  Erosion  decreases  as  the  speed  decreases.  If  e  <  0.5,  then 

m  must  vanish  in  finite  time,  presuming  the  target  is  thick  enough,  with 
the  erosion  rate  increasing,  decreasing  or  remaining  constant  during  the 
motion.  Thus  the  model  offers  a  variety  of  possibilities.  The  case 
eis  =  0,  that  is  no  modification  of  the  coefficient  b^  in  Equations  (100) 

and  (102),  is  a  particular  case  of  e ^  <  0.5  and  so  cannot  be  used  to 

describe  penetration  in  a  semi- inf inite  target  since  m  does  not  always 
vanish  in  such  a  target. 

A  similar  analysis  can  be  made  for  m  >  0. 

B.  A  Classical  Application:  The  Oscillator  with  Diminishing  Mass 

A  well-known  classical  problem  is  the  constant  mass  harmonic  oscilla¬ 
tor  in  a  forcing  field  and  a  wide  variety  of  electrical  and  mechanical 
problems  have  been  solved  by  applying  the  solution  of  the  basic  differen¬ 
tial  equation.  Suppose  the  oscillator  mass  is  not  constant  but  changes 
during  the  motion.  How  does  the  solution  change?  An  answer  has  been 
given  in  the  previous  section  which  we  will  now  apply  to  a  particular 
case  which  can  be  easily  tested  in  the  laboratory. 

Consider  a  liquid-filled  vessel  of  constant  interior  cross-sectional 

area  A,  height  h  ,  and  mass  m^  =  m  +  m  suspended  in  the  earth’s 
o  lo  v  o 

gravitational  field  by  a  spring  of  negligible  mass  and  friction  and  with 

spring  constant  k.  Here  m  is  the  constant  mass  of  the  vessel  while  m 

v  o 

is  the  initial  mass  of  the  liquid  so  m^.  is  the  initial  total  mass.  At 

time  zero  the  vessel  is  released  from  its  position  a  distance  S  below 
its  equilibrium  position  and  simultaneously  a  hole  of  cross-sectional 


66 


area  a  <  A  is  opened  in  the  bottom  of  the  vessel.  The  mass  of  the  liquid 
at  any  time  is  m  =  p  A  h  where  p  is  its  constant  density  and  h  is  the 

height  of  the  liquid  at  time  t  >  0.  The.speedof  recession  of  the  liquid 
surface  relative  to  the  vessel  is  v  =  -  h  =  -  m/ (pA) .  From  Bernoulli's 
principle  the  exit  speed,  v£,  is  related  to  v  by 


1  P  (v  2  -  v2)  =4  P  v  2  [1  -  (a/A)2]  =  p  gh  =  p  g  (m/p  A)  (124) 

2  fc  2  t 

since  v  =  (a/A)  v„  by  the  equation  of  continuity.  Since  the  total  mass 
E 

is  irij,  =  mv  +  m  so  m^  =  m=  -pAv;:-pa  v^,  then 


m  =  - 


V 


2  p  g  A  (a/A)  1/2 

[1  -  (a/A)2] 


V 


m  2g  (a/A)  j/2 

5 - 2“  m1/Z.  (125) 

hQ[l  -  (a/A)2] 


If  we  integrate  this  equation  we  obtain 


nvj,  -  my  =  m  =  mc 


1  -  2  - 


i_  in. 

"•  V  h0 


(a/A)2  (p  A  hQ) 2 
[1  -  (a/A)2] 


(126) 


which  has  the  same  form  as  Equation  (114)  above  with 


m  *  - 
o 


2  g  (a/A)2  (p  A  hQ)2 
hQ  [1  -  (a/A)2] 


m 

iO 


i 


2  g  (a/A)' 


-  2  K  m 


hQ  [1  -  (a/A)  ] 


(127) 


which  defines  the  constant  K.  From  Equation  (127)  we  can  determine 

(-  m  ),  knowing  mo>  g,  a,  A  and  hQ.  A  specialized  form  of  Equation  (102) 

can  be  used  to  describe  the  system  we  are  considering,  namely,  an 
undamped  (bJs  =  b2s  =  0),  harmonic  (dj  =  k,  d2g  =  0)  oscillator  with 

diminishing  mass  (m^,  =  m  <  0)  in  the  earth's  gravitational  nelu 


(as  mg) .  Thus  we  have 
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(128) 


mT  s  +  (i  -  Ej)  nir  s  +  k  s  =  (-  jyr). 

If  we  let  4>  =  4> (t)  be  an  arbitrary  transformation  of  the  independent 
variable  t  as  before,  we  obtain 


itl 


d^s 

d<fi2 


m  *  d4>  ds  ,  . 

‘"T  T  d4>  d<}>  +  ^  ~  V 


nu 


ds 

d(J) 


+  k  s  =  -  rm  g  (129) 


or 


"jg) 


(130) 


where  Equations  (104)  and  (105)  for  the  constants  M  and  B  are  in  this 
case  M  =  ^  =  mTo  0q2  =  =  ry  +  my  and  Bj  =  [jy  dj/d*  +  (1  -  y) 

mTJ  *  =  "  f£i  '  J-1  mTo’ by  Rations  (105)  and  (108).  In  this  case  we  have 
from  Equations  (111)  and  (127) 


*  =  (1/K)  £n  (1  -  Kt)"1 
and  by  Equations  (114),  (126)  and  (127) 

ry  -  ly  =  m  =  mQ  e'2  K4>  =  mQ  (1  -  Kt)2 
then  Equation  (130)  becomes 


.2 

M  d  s 
M  — -  + 

d<T 


ds 

d$ 


+  k  s 


'  8  ("V  +  mo 


e’2  K*). 


(131) 


(132) 


(133) 


The  general  solution  is  the  sum  of  the  solution  of  the  homogeneous 
equation  (right  side  of  Equation  (133)  set  equal  to  zero)  (given  by 
Equation  (45)  with  6=0)  and  a  particular  solution.  In  other  words 
our  solution  has  the  form 


s 


+ 

Y  <!> 

C  e  + 


C2  e 


Y  $ 


(C3  *  C4 


e-2K*) 


(134) 
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If  we  substitute  Equation  (134)  into  Equation  (133)  and  require  the 
coefficients  of  the  exponentials  to  vanish  we  find 


+ 

Y 


-  k/M  =  - 


a  ± 


i  0) 


(135) 


C3  =  -  my  g/k 


(136) 


c4  =  -  mo  g/[4  M  K2-  2  K  Bj  +  k]  (137) 

with  C  and  C?  to  be  determined  from  initial  conditions. 

If  we  hang  the  empty  vessel  on  the  spring  it  will  extend  a  distance 
(-  C3)  which  determines  k  in  Equation  (136)  once  we  measure  and  g. 

If  we  fill  it  with  a  liquid  mass  mQ  and  keep  the  hole  closed  so  a  =  K  = 

0  it  will  extend  an  additional  distance  mQ  g/k  which  gives  us  a  second 

measure  of  k  once  we  measure  mQ.  If  we  open  the  exit  hole  and  observe 

the  relaxation  of  the  system,  the  motion  is  described  by  the  particular 
solution 


s  =  C3  +  C4  e'2U  =  -  (gMniyA  +  m0  C1  '  Kt)2/(4  M  **  ~  2  K  B1  +  101 ' 

(138) 

At  the  end  of  the  relaxation  time,  1/K,  or  <t>  =  00  the  vessel  is  empty 
and  motionless  at  s  =  -  g  THy/k.  If  we  measure  s  at  various  times  during 

the  motion,  we  can  determine  -  y)  (-  mQ)  =  "  y)  2  m0 

that  is,  determine  since  we  know  m^  and  K.  Now  we  know  and  C^. 

If  we  refill  the  vessel  as  before  and  pull  it  down  an  additional 
distance  (-  S),  then  release  it  as  we  simultaneously  open  the  exit  hole, 
an  oscillatory  motion  is  superimposed  on  the  relaxation  motion.  The 
initial  conditions  are 

So  =  C1  +  C2  +  (c3  +  c4}  =  '  s  +  (C3  +  C4)  (139) 

and  s  =  4>  (ds/d4>)  =  (ds/d<M  =  0,  which  gives 

o  o  o  o 
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0  =  c  v  +  r 

1  Y  +  C7  Y  -  2  KC  , 

4 

Fr°"  EquatiOTs  (139)  and  (140)  we  find 

C1  =  (2  K  C4  +  Y-  S)/(v*  _  y-) 

C2  ’  <2  t  C4  *  /  S)/(y*  .  v-, 

S°  Equati°"  (134)  becomes 


C140) 


(141) 

(142) 


+  C3  +  C4  e 


4 

-2  K<j> 


This  describes  the^Tic^  unt  ilTh3^  ^  ^  giVe"  b/  Ecluati°n  (137) 

‘7ld  1S  *«  «  t  .  msTc^r  ff  ?*  8one.  Agai„,  wbln  ^ 

an  the  motion  ceases.  For  t  <  l /v  "V  g-  k  h>'  Equation  (143) 

described  by  t  <  UK  the  decay  of  the  amplitude  is 


(143; 


where 


^  =  (1  -  Kt)a/E 


a 

K 


(c 


1 


')  (-  in  ) 
o 


2K~r\n 


(mo  *  ”\  )  <C1  ■  -rl 


(144) 


(145) 


so  observing  the  amplitude  oo  .  ^ 

measure  of  e,.  If  ^  would  also  givc  a 

ecays  Unearly  with  time  a,  a  constant  rate  ,(  ,,  ‘""PlitUdo 

swfif  ZStH  '.bu"  sr»™w,ih  P '  d;P ' 

tion  (143)  o  =  °*  ^en  0  =  t  and  from  Equa- 


where  u  =  4c/m 

u  To 

mass  with  constant 


lo 


describing  undamped  harmonic  motion  ( 
*■  "‘To  A  measurement  of 


(146) 

constant 

2tfv  =  2tt/t 
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where  v  is  the  frequency  and  t  is  the  period  in  real  time  will  give 
another  measure  of  k. 


When  the  hole  is  opened  at  t  =  <f>  =  0,  id  is  a  constant  only  in  the 
transformed  time  <J>.  We  can  determine  the  turning  points  of  the  motion 


by  setting  s  =  <J>  ds/d<(>  =  0  (or  ds/d(j>  =  0  since  <J>  = 


(1-Kt) 


0) .  In  the 


s-<J>  plane  the  period  =  2tt/u)  where  to  =  yk/mTo  -  [Bj/  (2  m?o)  ] c  is 


constant  and  successive  minima  occur  for  <J>  =  n  x,  =  n  (27r/oj)  while  real 

n  <p 

time  successive  minima  occur  at 


tn  =  (1  -  e'Kb/K  -  0  -  e-nK2l,/“)/K 


(147) 


with  n  =  0,  1,  2  ....  using  Equation  (131).  As  can  be  seen  the 

difference  t  ,  -  t  becomes  smaller  as  time  increases.  As  expected, 
n+1  n 

the  oscillation  becomes  more  rapid  as  the  mass  decreases. 

This  latter  point  was  checked  experimentally  by  suspending  a  16  gm 
vessel  (waxed  cardboard  carton)  of  cross  sectional  area  A  =50  cm  and 
height  h  =  9.8  cm  from  a  spring  with  constant  k  =  1.5  x  10  dyne/cm 
determined  by  hanging  various  weights  upon  it  and  observing  its  extension. 
This  vessel  was  filled  with  490  gm  of  soapy  water  (to  improve  the  flow) 
and  a  real  time  period  determined  by  recording  the  time  at  which  every 
fortieth  minimum  occurred  and  averaging  over  ten  such  measurements.  The 
time  was  measured  to  the  nearest  0.01  second  using  the  crystal  controlled 
timer  built  into  the  Hewl itt-Packard  Model  55  hand  held  electronic 
calculator.  The  memory  was  used  to  record  and  average  the  times,  giving 
an  average  period  of  0.413  ±  .002  sec  with  the  actual  difference  between 
high  and  low  values  being  0.2/40  =  0.005  sec.  If  =  A/m^  then 

k  =  mTo  t »  2  =  506  (2n/ . 413) 2  =  506  (231.451)  =  1.17  x  105  which  is  22% 
lower  then  the  value  measured  by  extension.  The  difference  is  due  to 


friction  which  we  have  been  neglecting.  Actually  w  =  ^k/mTo  -  [b1/2mTo]  , 


implying  that  b^  the  friction  coefficient. 


is  2  mTo 


8.16  x  10^  gm/sec  using  k  =  1,5  x  10^  dyne/ cm  so  that  b-  is  not 

negligible  and  must  be  known  if  we  wish  to  calculate  the  observed  decay 
in  amplitude.  For  this  system  with  m^  =  506  gm, oscil lations  could  be 

observed  for  more  then  five  minutes,  but  the  system  appeared  stationary 
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after  7  minutes.  A  solid  brass  weight  of  500  gm  was  observed  to 
oscillate  for  about  10  minutes  on  the  same  spring,  indicating  a  damping 
contribution  from  the  liquid  motion.  Still,  we  don't  need  to  know 

in  order  to  check  the  fact  that  is  constant  for  a  constant  mass.  In 

this  case  =  2u/(.413)  =  15.21  sec-1, 

A  circular  hole  was  opened  in  the  bottom  center  of  the  vessel  and 
the  drain  time  (until  dripping  began)  was  found  to  be  1/K  =  250  sec, 
so  K  =  .004  sec'l.  Thus  the  effective  diameter  was  0.19  cm  and  the 
effective  area  of  the  hole  was  0.028  cm2  which  is  only  35%  of  the  actual 
area  of  0.079  cm2  (actual  diameter  0.3175  cm)  because  of  the  well-known 
vena  contracta  effect.  The  filled  vessel  was  pulled  down  a  distance 
S  =  1  cm  and  released  at  time  zero  while  the  hole  was  opened  simulta¬ 
neously.  As  before,  the  time  at  which  every  fortieth  minimum  occurred 
was  measured,  giving  the  values  in  Table  II. 

Table  II.  Times  (sec)  for  Minima  of  Decreasing 
Mass  Oscillator 


n 

0 

40 

80 

i 

20 

160 

200 

240 

280 

320 

360 

t 

n 

(obs) 

0 

15.98 

31.33 

45 

.77 

59 

.38 

71.82 

83 

.90 

95 . 36 

106.59 

117.41 

t 

n 

_tn-40 

- 

15.98 

15.35 

14 

.44 

13 

.61 

12.44 

12 

.08 

11.46 

11.23 

10.82 

-t  ,n)/40 
n  n-4U 

- 

.400 

.384 

.361 

340 

.311 

362 

.286 

.281 

.270 

t 

n 

(calc) 

0 

16.25 

31.44 

45 

.64 

58 

.93 

71.34 

82 

.96 

93.81 

103.96 

113.45 

if 

we  extrapolated  (t 

-  t 

n  - 

•  40>/40 

to  1 

Lime  zc 

>ro 

we 

would 

obtain 

0.415 

which  is  the  time  for  one  vibration  with  the  constant  mass  system.  The 
observed  increase  in  vibration  frequency  agrees  at  least  qualitatively 
with  the  model.  If  we  wish  to  calculate  u)  in  Equation  (147),  namely. 


w  =Vk/mTo  -  [<b!  +  mo)}/(2  mTo)]2  (149) 

3 

we  can  use  b,  =  8.16  x  10  gm/sec  as  above.  Since  (-  m  )  =  2  K  m  = 

1  6  ^0^0 

3.92  gm/sec,  must  be  considerably  larger  than  1/2  to  have  much  effect 

on  the  value  of  uk  In  order  to  measure  e  we  would  need  precise  equip- 

ment  since  k  >>  2  K  B  in  Equation  (133)  for  f  <  10  .  If  we  assume 

1  -l  ^ 

=  100,  then  w  =  15  sec  in  Equation  (149)  which  is  not  much  less 
than  u)j.  =  2tt/(.413)  =  15.21  sec  b  If  we  use  w  =  15  sec  ^  and  K  =  .004 
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sec  *  in  Equation  (147)  we  can  calculate  the  t  values  given  in  the  last 

line  of  Table  II.  These  are  in  reasonable  agreement  with  the  observed 
values,  although  the  agreement  becomes  worse  as  the  motion  damps  out. 

A  slightly  different  choice  for  could  have  made  the  observed  and 

calculated  values  agree  initially.  If  the  spring  were  weaker  or  the 
mass  heavier,  better  agreement  might  be  expected  for  a  longer  time. 

Even  the  classical  damped  harmonic  oscillator  with  constant  mass 
requires  a  theoretically  infinite  time  for  its  amplitude  to  decay  to 
zero  and  so  does  not  perfectly  describe  real  mechanical  (or  electrical) 
oscillators  which  cease  their  motion  in  finite  time  because  of  factors 
not  adequately  represented  by  the  simple  equation  used.  Similarly,  we 
cannot  expect  a  perfect  description  of  a  variable  mass  oscillator  by  a 
simple  amendment  of  the  constant  mass  equation. 

We  have  given  this  illustration  to  show  that  the  proposed  amendment 
to  the  classical  equation,  namely  the  added  damping  factor  (-  m) ,  can 

provide  a  reasonable  description  in  at  least  one  simple  case  up  to  360 
oscillations.  Our  intention  here  is  to  use  this  amendment  to  describe 
an  eroding  penetrator  for  which  the  motion  should  be  completed  in  less 
than  one  cycle.  Thus,  our  expectation  is  that  the  model  will  be  good 
enough  for  the  practical  calculations  we  have  in  mind.  Of  course  only 
comparisons  with  experimental  projectile-target  data  will  tell  us  how 
close  this  agreement  will  be.  We  will  make  such  comparisons  in  the  next 
subsection. 

C.  Applications  to  Target  Penetration 

Tate04  has  given  time-dependent  data  for  the  penetration  ana  erosion 
of  dural  and  aluminum  rods  striking  a  semi-infinite  polyethylene  target 
at  zero  obliquity.  The  5.45  gm  rods  were  6.35mm  m  diameter  and  63.5mm 
long  while  the  striking  speed  was  1.646  km/s.  Figure  5a  gives  the 
position  of  the  tip  of  the  dural  rod  as  a  function  of  time.  Tate's 
data  also  gave  the  position  of  the  rear  of  the  rod  as  a  function  of  time, 
but  this  is  not  shown  to  avoid  crowding.  Instead  the  difference  between 
the  tip  and  rear  which  is  the  remaining  length  is  shown  in  Figure  5b. 

For  the  dural  rod  the  final  depth  of  penetration  of  the  tip  was  about 
155mm  while  the  remaining  length  was  18mm,  Figures  6a  and  6b  show 
similar  data  for  an  aluminum  rod  which  had  a  final  depth  of  penetration 
of  about  120mm  and  a  final  length  close  to  zero  (unobservable).  Tate's 
comment  was,  "The  whole  of  the  aluminum  rod  is  used  up  in  the  penetration 
process". 


34A.  Tate,  "A  Theory  for  the  Deceleration  of  Long  Rods  After  Impact," 
J.  Mech.  Rhys.  Solids  IS,  387,  1967. 
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Figure  S. 


Dural  rod  penetrating 
calculated  values. 


polyethylene: 


measured  and 
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Figure  6.  Aluminum  rod  penetrating  polyethylene:  measured 
and  calculated  values. 
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0  for 


Let  us  apply  our  model  to  Tate's  data.  We  recall  that  d?7  = 

a  cylindrical  rod.  In  addition  we  will  neglect  the  force  term 
(assuming  b^_  v  <<  B,_).  The  displacement  of  the  center  of  mass  is 

(omitting  the  subscript  z) 

z  =  e  ac^  [{ vo  +  6  (-  a  -  3)}  -  {v^  +  6  (-  a  +  3)1  c  +  5 

(ISO) 

from  Equation  (46)  with  4>(t)  instead  of  t  given  by  liquation  (111).  If 
the  rod  density  and  cross  section  are  constant  during  the  motion  we  have 
from  Equation  (113) 


m/n’o  =  ^/f0  =  [(1  -  bj/Bj)  exp  {-  <f>/[2  mQ  -  *)/B1  | }  +  b]  /  ] 2 .  (1 51 ) 

We  recall  that  B  =  b  +  (e  -  L  (-  m  )  and  will  take  c.  >  i  so  B,  >  0. 

111^0  I  2  l 

If  we  assume  that  the  viscosity  of  polyethylene  is  v  =  0.5  m2/s  =  .005 

2 

cm  /ps  as  was  done  by  Walters  and  Majerus28  based  on  analogous  Russian 
measurements20 

b  ~  6ttRPv  =  6tt  (.635/2)  (2.7)  (.005)  -  .08  gm/us  .  (152) 

As  before  in  our  application  to  the  constant  mass  data  of  Backman  and 
Finnegan  we  will  take  b]  to  be  less  than  the  estimated  value,  namely,  b} 

.04  gm/us.  From  Figures  5  and  6  we  can  estimate  m  =  -  .07  gm/us  and  the 

o 

final  position  of  the  center  of  mass  6  =  155-18/2  =  146mm  for  the  dural 
rod  or  120-0/2  =  120mm  for  the  aluminum  rod.  The  depth  of  penetration 
of  the  tip  will  be  P  =  z  +  i/2  in  either  case.  For  the  dural  rod  we 

-  5  2 

will  take  =  1  and  d^  =  10  gm/ (us)  ,  while  for  the  aluminum  rod, 

=  4  and  a =  3.5  x  10  ^  gm/(us)^. 

The  solid  lines  in  Figures  5  and  6  give  the  calculated  results  which 
compare  reasonably  well  with  experiment.  The  overestimate  of  P  at  early 
times  may  be  due  to  our  neglect  of  the  v^  force  term  which  would  be  more 
important  at  retarding  penetration  during  the  early  stages  when  the 
speed  is  high.  It  is  of  course  desirable  to  include  this  term  (though 
not  the  d?  s^  term)  and  future  work  will  discuss  methods  of  doing  this. 
For  the  present  we  will  compare  the  solution  for  the  force  -  (aQ  + 

v  +  d  s)  given  by  Equation  (150)  with  the  solution  for  the  force 
A  2 

-  (bj  v  +  b^  v  )  given  by 
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z  =  [Jin  {(1  +  vQ)  exp  (bj  t/mQ)  -  b2  v^bj}  -  4>/mQ]  (153) 


where  we  must  take  aQ  =  0  to  describe  embedment  as  explained  above. 
From  Equation  (153)  we  see  that  as  t  -*■  °°,  $  *  °° 


m^ 

2  *  2MAX  =  *"  (1  *  b2  W'  (154) 

From  Equation  (154)  we  see  that  2.im  z,,,„  =  m_  v_/b,  so  that 

,  _  MAA  O  O  1 

b2° 

b  <  m  v  /(z„.v)  observed  to  keep  b~  >  0.  For  the  dural  rod  for  example, 

bj  <  5.45  gm  x  1 . 646mm/ys/ 146mm  =  .06  gm/ys  while  for  the  aluminum  rod, 

b  <  5.45  gm  x  1 . 646mm/ys/120mm  =  .075  gm/ys.  In  either  case  bj  =  .04 

gm/ys  as  above  is  small  enough.  If  we  use  this  value  for  b^  and  the 

same  values  for  as  before,  namely,  =  1  for  dural  and  =  4  for 

aluminum,  then  (since  mQ  is  the  same)  it  is  clear  from  Equation  (151) 

that  the  calculation  of  i  will  not  be  different  than  before  since  it  is 
independent  of  a^,  d ^  or  b2*  The  same  is  true  for  4>  as  is  clear  from 

Equation  (111).  Only  z  and  therefore  P  will  change.  The  dashed  line  in 
Figure  5  gives  the  result  calculated  from  Equation  (153)  when  we  use 
b^  =  .03  gm/mm,  while  the  dashed  line  in  Figure  6  does  the  same  for 

aluminum  with  b2  =  .052  gm/mm.  There  is  little  to  recommend  one  force 

over  the  other  for  embedment . 

The  values  of  the  parameters  ,  mQ,  mQ,  vq  and  6  =  -  aQ/dj  were 
fixed  by  experiment,  that  is,  a,  vq  and  6  were  fixed  in  Equation  (150) 
as  were  4  or  mQ  and  b^  in  Equation  (151).  The  parameter  (  and  so 
Bj)  was  adjusted  in  Equations  (111)  and  (151)  to  give  reasonable  agree¬ 
ment  with  experimental  SL  vs.  t  values.  The  parameter  d^  (and  so  B)  was 

adjusted  to  give  reasonable  agreement  for  z  (or  P)  vs.  t.  Similarly 
b2  was  adjusted  in  Equation  (153).  In  neither  case  was  there  an  effort 

made  to  obtain  a  best  fit  since  our  purpose  here  is  merely  to  illustrate 
the  potential  of  the  method.  Future  work  should  include  fitting  a 
variety  of  data  and  codification  of  the  parameters. 
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Studies  of  penetration  into  a  semi-infinite  target  are  not  commonly 
made  by  those  who  study  KE  penetrators  in  this  country  (in  contrast  to 
those  who  study  shaped  charge  jet  penetrators).  usually  projectiles  are 
fired  through  a  plate  target  of  fixed  thickness  and  their  residual  speed 
and  mass  is  determined  as  a  function  of  striking  speed.  For  example, 

Herr  and  Grabarek22  have  observed  that  for  both  steel  and  tungsten  rods 
striking  steel  armor  plates  end  on  at  zero  obliquity  the  residual  mass 
increases  as  the  striking  speed  increases  in  the  ordnance  range.  This 
is  illustrated  in  their  Figure  16  for  a  7.78  gm  steel  rod  with  a  length 
of  50mm  and  a  diameter  of  5mm  perforating  a  6.35mm  thick  steel  plate. 

This  behavior  is  entirely  reasonable  for  erosion  since  for  higher  v  . 

-  o 

less  time  is  spent  in  the  target.  If  the  initial  erosion  rate  is  not  a 
strong  function  of  vq  and  declines  with  time,  then  less  erosion  will 

occur  as  less  time  is  spent  in  the  target. 

Let  us  use  the  force  -  (a^  +  b^  z  +  z)  to  describe  their  data. 

We  can  determine  the  amount  of  time  or  transformed  time  <j>  spent  in  the 

target  by  setting  z  =  T  in  Equation  (150)  for  a  given  v  >  v  =  0.64 

o  oPL 

mm/ps  measured  by  Herr  and ^Grabarek .  We  recall  that  at  the  perforation 
limit  speed  z  =  6  =  T  and  z  =  0  as  4>  «>.  The  ricochet  limit  speed  was 

not  measured  but  was  taken  to  be  v_ni  =  0.55mm/ps.  In  addition  b  was 

Ut\L,  i 

estimated  to  be  ~  6irRpv  =  6tt  (.635/2)  x  7.8  x  5.5  x  10“^  -2.5  gm/ys 
using  the  Russian20  value  for  steel  of  v  =  5.5  M2/s  =  .055  cm2/ps. 

Instead  of  cutting  b  by  half  as  we  did  for  aluminum  targets,  it  seems 

more  appropriate  to  divide  by  five  and  take  b  =  0.5  gm/'ps.  This  may  be 

due  to  the  approximate  nature  of  our  model  (in  particular  to  our  neglect 
2 

of  v  terms)  or  perhaps  partly  due  to  a  lack  of  Stokes1  law  transfer  from 
the  explosive  plate  measurement  technique  used  by  the  Russians  and  long  rod 
penetration  of  steel  plates.  Such  questions  require  much  more  investi¬ 
gation  and  will  be  left  for  future  work.  The  actual  values  used  for  the 

parameters  were  bj  =  0.5  gm/ys,  =  7.78  x  lO-'”*  gm/(ys)2,  t-  =  1  and 

mo  =  -  0.4  gm/ps.  Once  <p  is  determined  for  a  given  v^,  then  Equation 

(151)  is  used  to  calculate  the  residual  mass  and  the  derivative  of 
Equation  (150)  is  used  to  calculate  the  residual  speed,  namely. 

•  • 

z  =  <(>  dz/d(p 

=  V'm’  hi’  ^v0-6  (a+3P  (a-6)e  ^  +  {vq-6  (a-g) }  (a+8)e 

(155) 
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where  we  have  used  Equation  (104a)  for  <J>.  The  results  of  these  calcu¬ 
lations  are  compared  with  the  data  of  Herr  and  Grabarek  in  Figures  7a 
and  7b.  Table  III  also  shows  how  <J>  decreases  as  vq  increases  in  agree¬ 
ment  with  the  idea  that  less  time  in  the  target  means  less  erosion. 

Real  time  t  behaves  in  a  similar  way  and  can  be  found  from  Equation  (111) 

although  there  is  no  need  to  do  this  in  order  to  carry  out  the  calcula¬ 
tion. 

Since  we  are  dealing  with  a  case  of  rectilinear  motion  in  the 
perforation  regime  only, we  can  also  describe  the  data  by  using  the  force 

-(aQ  +  bj  z  +  b^  z2)  provided  we  require  aQ  to  vanish  when  vq  =  v^^- 

Let  us  take  a  =  a,  T  (v  /v  -  1)  with  a,  small.  Now  we  determine  4> 

o  1  0  oFL  i 

for  given  v^  by  letting  z  =  T  in 

z  =  (^)  In  [o  -  Ae'EV(l  -  A)  *  *]  .  (1S6) 

the  analog  of  Equation  (24)  where  (-q)  =  b^  -  4aQ  b^,  h  =  -  b^  ±  /-q, 

A  =  (2b.,  v^  -  h+)/(2b0  vn-  h')  and  B  =  /^q/mft.  Once  4>  is  determined 

then  Equation  (151)  is  used  to  calculate  the  residual  mass  and  the 

derivative  of  Equation  (156)  is  used  to  calculate  the  residual  speed, 
namely, 

z  =  *  dz/d*  =  ^  (4b  )  (h+  -  A  h"  e~BVd  -  A  e'B4>).  (157) 

O 

If  we  choose  =  ,001  gm/ (ps)  ,  b-  =  0.3  gm/us,  ^2  =  ~  ^ 

and  mQ  =  -  0.4  gm/ps  we  obtain  the  dashed  lines  in  Figure  7.  The 

calculation  for  the  residual  mass  is  so  close  to  the  previous  calcula¬ 
tion  that  the  solid  and  dashed  curves  are  indistinguishable.  Again  no 
attempt  has  been  made  to  optimize  agreement  with  experiment.  We  note 
that  an  even  lower  value  of  b seems  appropriate  when  a  v2  term  is  in¬ 
cluded.  If  we  limit  ourselves,  to  zero  obliquity  perforations  as  studied 
by  Tate  or  by  Herr  and  Grabarek,  there  is  little  to  recommend  the  ^ 
component  force  -(aQ+  b  z  +  dx  zj  over  the  force  -(aQ+  b1  z  +  b2  z  ) . 

However,  if  we  wish  to  include  richochet,  embedment  and  perforation  in 
one  unified  theory,  then  the  former  force  is  preferable  as  we  have 
seen. 
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(mm  //ts ) 


Figure  7,  Residual  speed  and  mass  versus  striking  speed  for  a 
long  steel  rod  penetrating  a  steel  plate:  measured 
and  calculated  values. 
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Table  III.  Calculated  Values  for  Rod  Mass  and  Speed 


Using  the 

Force 

-Ca  +  bj  z 

+  dj  z) 

v  mm/ys 

0 

.64 

.7 

1.0 

1.3 

1.6 

<f>  VS 

oo 

15.05 

7.85 

5.56 

4.34 

m/mo 

.51 

.62 

.73 

.78 

.82 

v,  mm/vs 

0 

.29 

.76 

1.13 

1.47 

In  Herr  and  Grabarek's  report  there  are  some  indications  that  for 
some  materials,  perhaps  tungsten  rods  versus  steel,  the  residual  mass 
might  first  increase  and  then  decrease  somewhat  as  vq  increases.  The 

data  is  not  good  enough  to  say  that  this  occurs.  If  a  case  is  found 
in  which  this  behavior  occurs,  then  the  above  model  can  be  modified 
by  the  addition  of  a  simple  ad  hoc  postulate  to  include  such  a  case. 
For  given  materials  we  might  expect  mo>  the  initial  erosion  rate,  to 

depend  somewhat  on  v  ,  the  striking  speed.  A  simple  representation  of 
this  might  be 


m  =  -  o  (v  -  v_)n  (158) 

o  v  o  E 

where  o  and  n  are  constants  and  v£  is  a  critical  striking  speed  below 

which  no  erosion  occurs.  When  Equation  (158)  is  put  into  Equation  (118), 
we  see  that  the  slope  in  an  m  versus  $  (or  t)  plot,  namely  m,  becomes 
more  negative  for  higher  v  .  This  behavior  is  shown  m  Figure  8  by  the 

two  curves  marked  vq^  (low)  and  v q2  (high).  For  the  low  striking  speed 

VqI  the  time  to  perforation,  t^,  will  be  greater  than  the  time  to 

perforation,  t 2>  corresponding  to  the  higher  striking  speed  vq2.  The 

intersections  of  vertical  lines  through  t^  and  t2  with  the  mass  curves 

give  the  residual  masses  at  perforation.  In  Figure  8  we  have  shown  the 
case  for  which  the  residual  mass  is  the  same  for  both  striking  speeds 
(that  is,  the  dashed  line  is  horizontal).  However,  it  is  easy  to  see 

that  a  more  negative  slope  on  the  v q2  curve  would  lead  to  a  lower 

residual  mass  after  perforation  than  for  the  lower  v  ^  if  remained 
the  same.  Similarly  a  less  negative  slope  on  the  vq2  curve  would  lead 

to  a  higher  residual  mass  than  for  v  ^  if  were  the  same.  Trial 

calculations  have  been  made  and  these  have  confirmed  the  ability  of 
such  a  modification  to  represent  the  occurrence  of  a  maximum  in  the 
residual  mass  versus  striking  speed  curve  with  leveling  off  to  a 
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Figure  8. 


Sketch  of  possible  mass  loss  behavior  ■ 
eros  rate  dependent  on  striking  spec 


initial 
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constant  value  at  higher  striking  speeds.  However,  these  calculations 
will  not  be  presented  here  since  we  have  no  convincing  evidence  that 
such  a  phenomenon  occurs. 
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VIII.  A  MODEL  OF  PROJECTILE  BREAKUP 


Here  we  are  thinking  of  breakup  as  a  discrete  or  discontinuous  process 
in  contrast  to  the  continuous  process  of  erosion.  In  recent  years 
microscopic  fracture  theories  have  been  advanced  and  linked  to  continuum 
fracture  mechanics  by  constructing  statistical  distribution  functions 
which  represent  molecular  bond-breaking  or  flaw  activation  rates.35 
In  this  way  more  physics  can  be  incorporated  into  a  theory  which,  however, 
remains  basically  statistical  in  nature.  This  is  partly  necessitated  by 
the  enormous  complexity  involved  in  giving  a  deterministic  description 
to  each  component  of  a  large  ensemble.  There  is  also  a  virtue  in  this 
necessity  since  the  behavior  of  large  ensembles  seems  to  follow  rather 
simple  laws  which  are  of  more  interest  than  the  behavior  of  the  indivi¬ 
dual  components.  Gredenko  has  expressed  this  very  well  in  discussing 
the  molecular  structure  of  matter: 

"The  problem  here  is  not  in  fact,  that  of  investigating  the  motion 
of  individual  particles  but  that  of  determining  the  laws  to  which  aggre¬ 
gates  of  a  large  number  of  interacting  particles  are  subject.  However, 
the  laws  that  arise  in  consequence  of  the  mass  character  of  the  parti¬ 
cipating  components  have  their  own  peculiarities  and  do  not  amount  to  a 
simple  summation  of  the  individual  motions.  Furthermore,  within  certain 
limits,  these  laws  are  found  to  be  independent  of  the  individual  propert¬ 
ies  of  the  particles  that  give  rise  to  them."36 

Early  models  of  fragmentation  such  as  that  proposed  in  England  by  Mott 
and  Linfoot37  were  also  statistical.  At  the  suggestion  of  B.  L.  Welsh 
these  authors  proposed  that  the  number  of  bomb  or  shell  fragments  with 
masses  between  m  and  m  +  dm  could  be  represented  by  the  distribution 


where  C1  and  a  are  constants  and  dn  is  the  number  in  the  mass  interval 
dm.  The  fact  that  this  distribution  fit  experimental  observations 


351>.  R.  Curran,  L.  Seaman  and  D.  A.  Shockey ,  "Dynamic  Failure  in  Solids" 
Physics  Today ,  Jan  1977 ,  pp  46-55. 

36S.  V.  Gnedenko,  The  Theory  of  Probability, Trans,  by  B.  D.  Seckler, 

NY,  Chelsea  Publ.  Co.,  1962. 

37W.  F.  Mott  and  E.  H.  Linfoot,  "A  Theory  of  Fragmentation, "  Advisory 
Council  on  Scientific  Research  and  Technical  Developments ,  A.  C. 

3348,  Jan  1943. 
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reasonably  well  suggested  the  following  idea  to  Mott  and  Linfoot.  If 
the  number  of  fragments  with  characteristic  length  x  in  the  interval 

1  “X/x 

between  x  and  x  +  dx  were  C  e  o  dx  and  the  length  x  were  proportional 
to  the  cube  root  of  the  mass  for  a  three-dimensional  body,  then 
Equation  (159)  would  result.  They  pointed  out  the  fact  that  similar 
distributions  had  been  used  to  describe  the  weight  or  diameter  distribu¬ 
tions  of  crushed  minerals  or  sand  particles38  and  discussed  the  one-  and 
two-dimensional  analogues  of  Equation  (159)  which  might  be  thought  of 
if  linear  or  planar  breakup  were  involved. 

Soon  after,  also  in  England,  Ursell39  suggested  that  the  chance  that  the 
number  of  fractures  be  n  ought  to  be  given  by  a  Poisson  distribution 

e”^  (X£)n/n!  for  a  rod  of  length  Z  with  the  expectation  of  fractures  being 
XZ.  Transatlantic  communications  must  have  been  excellent  for  in  less 
than  six  months  Thomas  published  a  report  in  the  United  States  showing 
that  the  results  of  Mott,  Linfoot  and  Ursell  could  be  obtained  from 
more  general  considerations  which  did  not  require  random  planes  traversing 
a  volume40.  He  pointed  out  that  Equation  (159)  (or  its  one-  and  two- 
dimensional  analogs)  does  not  depend  on  particular  breaking  mechanisms 
such  as  the  simultaneous  formation  of  all  fragments,  but  allows  for 
example,  smaller  fragments  to  be  formed  at  a  later  stage  from  larger 
fragments  already  formed. 

Thomas' 'observation  on  the  generality  of  Equations  of  the  type  of 
Equation  (159)  has  been  enlarged  upon  by  Molitz41  who  used  a  Weibull42 
distribution  to  describe  fragment  masses  from  shells.  This  has  a 
density  function 

dn  =  C  m*'1  e"a(mX)  dm  (160) 


which  reduces  to  Equation  (159)  when  X  =  1/3  (or  to  the  "two-  or  one 
dimensional"  analogs  when  X  =  1/2  or  1).  Of  course,  the  "one-dimensional" 


36Lienau3  J.  Franklin  Inst.3  (1935)  p  485 

3SH.  D.  Ursell 3  "Fragmentation  Data  and  Theories  of  Fragmentation 3 " 

A.C.  3817 ,  April  1943. 

40L.  ti.  Tnomas3  "Comments  on  Mott's  Theory  of  the  Fragmentation  of  Shells 
and  Bombs "3  BRL  R  3983  Sept  1943.  (ATI  #36152) 

41tf.  Molitz 3  "Einige  Bemerkungen  zu  den  Verteilungsfunktionen  der  Splitter 
grbssen3  Explosivstoffe  2_  Mar/Apr  1973s  p  33. 

h>1W.  Weibull 3  Statistical  Evaluation  of  Data  from  Fatigue  and  Creep-rupture 
Tests.  Fundamental  Concepts  and  General  Methods".  WADC  Tech. Rep t  59-400 , 
Part  I3  1959. 
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analog  with  A=1  is  the  exponential  distribution  while  the  case  A=2  is 
called  Rayleigh's  distribution. 

Gredenko43  showed  in  the  asymptotic  theory  of  extreme  values  that 
only  three  types  of  limiting  distributions  exist  with  Weihull's  formula 
being  the  third  asymptotic  distribution  of  the  smallest  values44.  In 
short.  Equation  (160)  is  a  generalization  (with  A  as  an  adjustable 
parameter)  which  does  not  depend  on  the  details  of  any  physical  model 
of  breakup.  In  fact,  it  has  been  applied  to  a  wide  variety  of  physical 
and  social  phenomena45. 

Since  practical  fragmentation  experiments  require  a  non- zero 
(measurable)  lower  mass  limit,  we  can  introduce  a  cutoff  mass,  m  , 
and  write  Equation  (160)  as  c 

dn  =  C  (m-rac)A  A  e  <it‘lu  racJ  dm  (161) 


which  is  the  usual  generalization  to  include  a  shift  in  the  origin.  Here 
A  is  an  adjustable  parameter  which  at  least  for  shell  fragmentation  by 
high  explosives  experience  indicates  is  less  than  unity  (near  1/2  or 
1/3). 

In  view  of  the  above  considerations  about  the  generality  of  Equation 
(160),  it  seems  reasonable  to  apply  it  to  projectile  breakup  under  impact 
as  well  as  to  shell  fragmentation  by  high  explosives.  When  a  projectile 
shatters  or  breaks  into  many  small  pieces  during  a  high  speed  impact, 
then  Equation  (160)  might  be  sufficient.  However,  for  impacts  at  ordnance 
speeds,  projectiles  tend  to  break  up  into  relatively  few  pieces  of  size 
comparable  to  each  other  and  the  original  projectile.  In  this  case  a 
discrete  distribution  function  related  to  Equation  (160)  should  be  more 
appropriate.  Let  us  trace  a  connection  between  certain  discrete  and 
continuous  distributions  which  will  be  suitable  for  our  purpose. 


43tf.  V.  Gredenko,  "Limit  Theroms  for  the  Maximum  Term  of  a  Variational 
Series ",  Doklady  Akad.  Nccuk,  USSR  32,  1941. 

44I?.  J.  Gumbel ,  "Statistical  Theory  of  Extreme  Values  (Main  Results)" 
o.6  in  Contributions  to  Order  Statistics ,  Ed.  by  A.  E.  Sarkan  and  E.  G. 
Greenberg  (NY:  John  Wiley  6  Sons ,  Inc.,  1962. 

45W.  Weibull,  "A  Statistical  Distribution  Function  of  Wide  Applicability" , 
J.  Appl.  Mech .,  Sept  1951 ,  p  293. 
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Consider  an  interval  y  and  divide  it  into  k  parts  each  of  size  e=y/k 
with  k  large  enough  so  that  subinterval  e  contains  on  the  average  only 
one  random  event.  Let  r  be  the  constant  average  rate  of  occurrence  of 
events  in  the  interval  y.  Then  the  probability  of  one  such  event 
occurring  in  the  sub interval  e  is  p  =  r  e  =  ry/k.  The  probability  of 
of  no  events  occurring  in  the  subinterval  e  is  (1-p).  If  these  events 
are  randomly  distributed,  the  probability  of  finding  exactly  s  of  them 
in  k  intervals  (k  trials)  is  given  by  the  binomial  distribution. 


q  (s;  k,  ry)  =  - — -  (ry/k)s  (l-ry/k)k"s  (162) 

s !  (k-s) ! 

where  s  =  0,l,2...k  and  q  =  0  for  all  other  s.  Now  consider  the 
limit  as  k+  00 ,  namely. 


lim 

k**°° 


q  = 


=  CT-LL 


=  P  (s;ry) 


[lim  - — - -J 

lim(l-ry/k)k 

|lim(l-ry/k)S  I 

lk-»-  oo  (k-s) '.  k.  J 

[k->  00 

rZ 

4 

8 

)  -  in  [ 

e“rXJ  [1] 

(163) 


which  is  the  poisson  distribution.  After  division  of  the  numerator 
by  the  demominator  the  first  limit  in  Eq  (163)  becomes 


lim  k(k-l)(k-2)...(k^s+l)  =  Um  i(iA)(l£\...(l-  (s-l)/kj=l  (164) 
k+  °°  ks  k+  *  k  k  ^ 


while  the  third  limit  is  obviously  unity.  If  we  let  z  ■ =  -(ry) /k,  the 

second  limit  becomes  the  (-ry)  power  of  the  limit  w.ac..  -e.mes  - 

of  the  natural  logarithm,  namely. 


lim 
z->  o 


(165) 


Finally  consider  the  probabily  that  one  or  more  events  will  occur  in 


the  interval  y,  namely. 


Y  p(s»ry) 

S=1 


Y  p  (s;ry)  -  P  (0;ry) 
s=0 


(166) 
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(167) 


by  equation  (163)  and  the  fact  that 

t  t  Util  -  <f  ^  = 

s=0  s>  s-u  s. 

From  Equations  (166)  and  (163)  it  is  clear  that  e-r^  is  the  probability 
that  no  events  will  occur  in  the  interval  y.  If  we  look  upon  ry  as  a 
continuous  random  variable  it  is  obvious  that  the  cumulative  probabilities 

— ry  -ry 

e  and  1-e  }  have  as  their  density  function 


d(n/N)  =  dn/N  =  e  d(ry)  =  r  e  dy 


(168) 


which  is  the  usual  exponential  distribution  giving  the  probability  that 
an  event  will  occur  in  the  interval  from  ry  to  ry+rdy.  Here  dn  is  the 
number  of  such  events  in  the  interval  dy  or  rdy  and  N  is  the  total  number  of 
such  events.  If  we  integrate  Equation  (168)  over  (ry)  from  0  to  ry  we 
obtain  Equation  (166),  while  if  we  do  the  same  from  ry  to  <*>,  we  obtain  e'rf. 
The  mean  value  of  y  is  1/r  while  the  median  value  y  is  (In  2)/r  =  . 7r. 


In  summary,  we  have  seen  that  the  Poisson  distribution  is  a  limiting 
form  of  the  binomial  distribution  for  small  p  (in  particular  for 
k>>  ry  >>  p  =  ry/k) .  We  have  also  seen  that  sums  over  the  discrete  Poisson 
distribution  which  divide  the  set  of  possible  outcomes  into  that  of  no 
event  occurring  and  its  complement,  that  of  one  or  more  events  occurring, 
have  the  form  of  a  cumulative  exponential  distribution  with  continuous 
density  function  given  by  Equation  (168) . 

The  above  analysis  has  been  deliberately  general  to  emphasize  the 
point  that  probability  distributions  of  random  variables  are  not  restricted 
to  any  particular  phenomenon  much  less  to  a  particular  mechanism  envisioned 
for  the  phenomenon.  They  are  rather  simple,  approximate  expressions  of 
laws  which  a  variety  of  large  aggregates  seem  to  follow  rather  well.  For 
example,  the  above  analysis  is  commonly  employed  in  discussions  of 
reliability  or  "time-to-first-failure"  problems.  The  interval  y  is  taken 
to  be  a  time  and  r  is  a  failure  rate  with  ry  being  the  expected  number 
°6  failures  in  time  y.  Similarly,  y  might  be  a  time  and  r  an  arrival  rate 
of  trucks  at  a  loading  dock  with  ry  the  expected  number  of  arrivals  in 
time  y.  In  another  application  the  interval  y  is  taken  to  be  a  volume. 

For  example,  if  microorganisms  are  distributed  randomly  in  a  body  of  water 
and  a  sample  drop  is  examined  under  a  microscope,  we  expect  to  find  ry 
organisms  if  r  is  the  density  or  number  per  unit  volume.  In  view  of  the 
great  variety  of  applications  for  such  formulas  it  does  not  seem  necessary 
or  appropriate  to  put  much  weight  on  discussions  of  dimensionality. 
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If  we  let  y  =  (m/P)^  in  Equation  (168)  it  becomes 

dn  =  N  e“r(-m/p-)  dTrCm/p)  X1  =  ■  -Ny--  n/-1  expj^-rmX/pX]  dm  (169) 

L  ^  p 

<which  is  Equation  (161)  with  C=NX(r/p*)  =  NXa.  Ijfhen  X=l/3  Equation  (169) 

becomes  Equation  (160)  with  C*=Na=N(p^3/r)  1  =  Ny  if  we  define  p=p/r'y. 

In  other  words,  Mott ' s  formula  is  a  particular  case  of  Equation  (161)  which 
in  turn  is  a  particular  case  of  Equation  (168)  which  is  an  exponential 
distribution  of  the  random  variable  y.  If  m  is  a  mass  and  p  is  a  mass  per 
unit  volume  then  x  can  be  interpreted  as  a  characteristic  dimension  of  the 
volume  in  question.  If  the  original  volume  is  cubical  and  breaks  into  cubes 
then  x  is  the  length  of  an  edge  and  X  is  1/3. 

If  the  volume  is  not  a  cube  or  does  not  break  into  cubes  but  the 
fragments  are  chunky  then  an  effective  value  for  X  should  not  be  far 
from  1/3.  If  the  volume  and/or  fragments  are  rods,  plates,  shells  or  some 
arbitrary  shapes,  then  the  a  priori  value  of  X  is  not  obvious.  If  Equation 
(169)  is  to  be  applied  to  such  a  volume  it  seems  better  to  treat  X  as  an 
adjustable  parameter  to  be  determined  by  sample  experiments. 

As  we  have  noted.  Equation  (168)  is  the  density  function  for  certain 
sums  over  the  Poisson  distribution.  It  is  instructive  to  return  to 
Equation  (162)  and  repeat  the  analysis  for  the  particular  application  of 
material  breakup.  In  this  case  k  is  the  very  large  number  of  defects  which 
are  found  in  real  materials,  while  s  is  the  number  which  are  activated  or 
converted  into  breaks  in  a  given  stressful  situation.  Now  r  is  the  average 
rate  of  occurrence  of  activated  defects  or  breaks  which  are  assumed  to  be 
randomly  distributed  over  the  interval  y  which  is  a  characteristic 
dimension  of  some  type.  Since  an  event  m  this  description  is  a  break, 
then  N  in  Equations  (168)  and  (169)  is  the  number  of  breaks  which  occur 
in  the  interval.  The  number  of  pieces  will  be  N+l  which  is  practically  equal 
to  N  for  N»l.  Here  r  or  ry  is  a  function  of  the  stress  applied  to  the 
body  in  question  and  will  depend  on  the  manner  in  which  it  is  applied 
as  well  as  on  its  magnitude.  In  addition,  it  will  depend  on  the  material 
properties  of  the  body  as  well  as  its  geometry.  In  the  case  of  an 
explosive-filled  shell,  N  is  usually  taken  to  be  equal  to  the  large 
number  of  fragments  collected,  X  is  taken  to  be  1/2  or  1/3,  and  y  is 
treated  as  an  adjustable  parameter.  Various  schemes  for  estimating  y 
for  certain  geometries  and  materials  have  also  been  proposed,  starting 
with  a  report  by  Mott1*6.  If  we  integrate  Eq  (169)  from  m^  to  m2  we 


^  .F’  Mott,  "Fragmentation  of  HE  Shells ;  a  theoretical  Formula  for  the 

Distribution  of  Weights  of  Fragments ",  A.C.  3642  (1943). 
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obtain  the  number  of  fragments  predicted  to  be  in  this  mass  range  for 
given  N,  A  and  y  X=^pr^j  •  If  we  let  m^=0  we  obtain  the  number  of  fragments 
which  should  have  mass  less  than  m  or  variable  less  than  w  =  ry 

A  2  2 

=r(m2/p)  ,  namely, 

™2  X  w 

N(m<ra2)  =  J  N  e"r(m/p)  d[r(m/p)X]  =  $  N  e"w  dw 

=  N[l-  exp  j  -r  ^  ~)A|]  =  N[l-e"W2]  (170) 

i 

with  N(m<m2)/N  giving  the  chance  that  an  individual  fragment  will  have 
mass  less  than  The  complementary  probability,  namely,  the  chance 

that  a  fragment  will  have  a  mass  greater  than  m2>  is  of  course 


N(m>m2)/N  =  j"  e  w  dw  =  exp  j  -r (m2/p)  X  j  =  e-W2  (171) 

w2  ' 

Here  m2  obviously  cannot  exceed  the  original  mass  of  the  body  before 

breakup  (mQ)  so  taking  the  upper  limit  of  the  integral  to  be  °°  is  an 

approximation.  However  the  error  involved  for  explosive  fragmentation 
is  usually  small  and  quite  acceptable  since  even  the  largest  fragment 
produced  is  very  much  smaller  than  m^.  The  error  is  also  acceptable  in 

view  of  the  simplicity  of  the  formula  and  the  experimental  accuracy 

which  can  be  achieved.  In  equating  N  to  the  total  number  of  fragments 

collected  we  are  effectively  defining  r  to  be  infinite  for  m^  m  so  that 

—  o 

Equation  (171)  vanishes  if  m2>m0  and  N(m<m0)  -  N  in  Equation  (170).  Other 

formulations  can  be  found  to  handle  both  a  cutoff  mass  for  particles  too 

small  to  measure  and  a  maximum  (unbroken)  mass  equal  to  m  .  These  will  be 

o 

the  subject  of  a  future  report  since  they  are  more  concerned  with  experi¬ 
ments  in  which  breakup  becomes  shatter,  that  is,  with  the  high  explosive 
fragmentation  of  shells  into  a  large  number  of  pieces  or  with  projectile 
breakup  upon  impact  in  the  hypervelocity  regime. 

% 

The  median  value  of  y,  namely  y,  is  found  by  setting  N (m<m2)/N=0 . 5  = 
N(m>m2)/N  in  Equations  (170)  and  (171)  and  solving  for  y2  =  y.  Thus 
y  =  (ln2)/r  as  before  so  the  median  value  of  the  mass  is 

111  =  P  (y)1/A  =  P  [(ln2)/rJ1/X  =(ln2)1/X  (p/r1/X)  (172) 
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The  mean  value  of  y,  namely  y,  is  found  from 

-  _  J  y  e'(ry)  d(ry)  =  r'1  (173) 

y  ~  0 

which  is  consistent  with  the  meaning  of  r  as  the  average  number  of  breaks 
per  unit  interval.  Then  the  average  value  we  expect  for  the  mass  is 

m  =  pty)1^  =  p/r1^  =  p  (174) 

by  our  previous  definition  of  p  if  X  =  1/3.  By  comparing  Equations  (172) 
and  (174)  we  see  that  the  mean  and  median  values  of  the  mass  differ  by  the 

factor  (ln2)^^  which  is  about  1/3  for  X  ■  1/3.  Having  in  <  in  is  character¬ 
istic  of  a  distribution  positively  skewed  (to  the  right) . 

In  this  report  we  are  interested  in  projectile  breakup  during  impact 
at  ordnance  speeds  where  the  number  of  breaks  is  not  much  different  from 
unity  and  the  mass  of  each  piece  is  not  very  much  smaller  than  the 
striking  mass  m^.  Consequently,  we  cannot  afford  to  set  the  number  of 

pieces  equal  to  the  number  of  breaks,  nor  can  we  let  the  upper  limit  of 
any  mass  integral  be  effectively  infinite.  We  can  however,  use  the 
Poisson  distribution  since  the  number  of  defects  or  unactivated  breaks,  k, 
is  still  very  large  (effectively  infinite)  in  any  real  material  so  that 
the  limiting  process  of  Equation  (163)  is  still  applicable. 

Equation  (163)  is  a  probability  distribution  for  the  discrete  variable 
s,  given  the  continuous  parameter  (ry) .  As  is  well-known,  the  mean  value 
of  s  is 


s  =  ry  =  Ng  (175) 

which  is  the  average  number  of  breaks  we  expect  in  the  interval  y  under 
stress  conditions  specified  by  ry.  Here  ry  is_  not  necessarily  an 
integer,  and  P  is  not  necessarily  maximum  at  s.  If  we  form  the  ratio 
P(s)/P(s-1)  =  (ry)/s  x  1  from  Equation  (163),  we  see  that  for  s<ry, 

P  increases  as  s  increases,  while  for  s>ry,  P  decreases  as  s  increases. 
If  ry  happens  to  be  an  integer,  then  P  has  two  equal  maxima  at  the  modes 

s  =  ry  and  s  =  ry-1.  Otherwise  P  has  one  maximum  at  s  equal  to  the 
largest  integer  which  is  less  than  ry  so  that  s  is  not  greatly  different 
from  ¥.  The  median  value  of  s,  namely  s  can  be  found  by  calculating 

the  cumulative  probabilities  or  partial  sums  S(£;ry)  =  I  P(s;ry)  for 

£=0,1,2, . .  until  we  find  a  sum  which  just  exceeds  0.5  with  S(£-l;ry) 
<0.5.  The  value  of  £  for  which  this  is  true  is  defined  as  5. 


Since  the  number  of  pieces  is  the  number  of  breaks  plus 
for  the  average  mass  _ 

mAV=nV(s+1)  =  m0/(ry+l)=  ra0/(NB+l) 


one  we  have 
(176) 
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from  Equation  (175) .  In  order  to  estimate  m^y,  let  us  make  some 


reasonable  assumptions  about  ry 


In  particular  let  us  assume  that 


is  proportional  to  an  energy  per  unit  volume  carried  by  the  projectile 
divided  by  the  energy  per  unit  volume  needed  to  fracture  the  projectile, 
namely, 


ry  =  NB  >  K[ip  (v0‘  -  O]  !oy  (177) 

where  the  density,  p,  the  striking  speed,  vn,  a  critical  speed,  v^, 

below  which  Nfi=0  and  the  yield  stress,  a  ,  can  all  be  determined 

experimentally.  Here  K  is  the  constant  of  proportionality  which 

measures  that  fraction  of  the  excess  striking  energy  density  which  is 

converted  to  internal  energy  of  the  projectile  and  is  responsible  for 

any  breaks  which  may  appear.  To  account  for  obliquity  it  is  reasonable 

to  assume  that  K  xs  proportional  to  the  line  of  sight  thickness 

T/cos  0  so 
o 


K  =  F/cos  0q  (178) 

where  F  is  a  fraction  dependent  on  factors  like  projectile  shape,  pitch 
and  yaw  as  well  as  target  hardness.  Instead  of  o  which  is  readily 

available  for  most  materials  of  interest,  it  might  be  more  satisfactory 
to  use  °pep»  where  <7^  and  are  the  true  stress  and  strain  at  fracture 

in  an  appropriate  experiment  which  combines  tensile  and  compressive 
forces.  However,  in  view  of  the  fact  that  ap  and  ep  are  not  readily 

available  for  all  materials  of  interest  even  in  simple  quasi-static 

tension  experiments,  we  will  use  a  ,  assuming  that  a  and  a„  are 

y  y  F  F 

related  for  particular  materials  of  interest. 

In  summary,  from  Equations^  (177)  and  (178)  we  can  calculate  the 
number  of  breaks  we  expect,  NR,  and  the  number  of  pieces  we  expect, 

Ng+l-  Then  we  can  find  the  average  mass  from  Equation  (176)  and 

estimate  the  variance  (or.  square  of  the  standard  deviation)  we  expect 
from  this  average  mass  from  the  well-known  fact  that  the  variance  of 
the  Poisson  distribution  is  equal  to  the  mean,  namely,  ry.  Alternatively, 
we  can  calculate  the  median  value  from  the  cumulative  probability 
S(£;ry)  or  even  the  mode  by  rounding  ry  to  the  nearest  integer  less  than 
or  equal  to  ry. 

Let  us  illustrate  the  above  model  by  applying  it  to  some  data 
contained  in  the  Holloway,  et  al .  report  cited  above1*.  Consider  a 
tungsten  sphere  of  mass  0.45  gm  and  diameter  0.36  9cm  striking  a  mild 


steel  plate  at  zero  obliquity.  All  of  the  zero  obliquity  data  for  this 
case  is  plotted  in  Figure  9,  a  total  of  39  points.  For  24  of  these 
points  the  target  thickness  was  0.16  cm,  for  14  of  them  it  was  0.32  cm 
and  for  one  it  was  0.64  cm.  All  but  three  involved  spheres  nominally 
0.45  gm,  while  three  were  nominally  0.26  gm.  This  information  is 
represented  by  the  different  symbols  used  in  the  figure.  From  the 

figure  we  can  take  vc=500  m/s  =  ,5xl05cm/sec  while  the  density  is  p=17.1 

gm/cm^  and  0^=6  kilobars  =  0.6x10^  dyne/cm^. 

If  we  assume  K=F=0.1_in  Equation  (178)  we  can  calculate  (Nfi)  as  a 
function  of  vq  namely,  Ng  =  .1  [±(17.1)  (vo2-.52)xl010J  /  .6*1010  = 

1.425(v  2-  .25).  The  results  are  given  in  Table  IV  with  the  ratio 


m,,,/m  from  Equation  (176)  plotted  in  Figure  9  as 
AV  o 

Table  IV 

the  dashed  line. 

vo(km/s) 

<0.5  .75 

1.00 

1.25 

1.50  1.75 

2.00 

CO 

nb+1 

0  .445 

1.069 

1.870 

2.850  4.008 

5.344 

oo 

1  1.445 

2.069 

2.870 

3.850  5.008 

6.334 

00 

mAV//mo 

1  .692 

.483 

.348 

. 260  . 200 

.158 

0 

Here  we 

are  neglecting 

mass 

loss  by 

erosion.  Below  v  =. 

0 

5  km/s  erosion 

is  quite  small  experimentally  and  from  considerations  given  in  the 
previous  section  it  should  be  even  smaller  for  higher  v^.  For 

comparison  we  have  also  plotted  in  Figure  9  the  residual  mass  curve  given 
by  the  Thor  Equation  in  the  Holloway  report4  namely. 


m=m  -10 
1  o 


4.3516,™.  .1455,,  x  .7191 


(T) 


(k) 


. .4048 r  .1. 3324,  ..9911  n ,Q. 

(sec  8q)  (mQ)  (vq)  (179) 


which  in  the  present  case  becomes  the  mass  ratio 


ml/mo  =  1  ‘  •315tVo)’"U  (180) 

with  v  expressed  in  k  m/s  where  we  have  used  k=20,  6q=0,  mQ=.45  gm 

and  T=1.6  mm.  For  T=3.2  mm  the  factor  .315  becomes  .348.  For  all 
practical  purposes  Equation  (180)  is  a  straight  line  and  predicts  that 

m  /m  vanishes  for  v  near  3.2  km/s  which  is  not  reasonable.  In 
r  o  o 
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Figure  9. 


Residua  1  mass  fract i on 

steel  plates  at  0  =  n° 

o 


versus  striking  speed  for 
:  largest  piece  found  and 


tungsten  spheres  perforating 
average  piece  expected. 


i 


contrast  Equation  (176)  vanishes  only  as  v  -►  ®.  As  noted  before.  Equation 

(179)  also  makes  the  unreasonable  prediction  that  ra^  will  decrease  and 

eventually  turn  negative  if  is  increased  while  all  other  parameters 

are  held  constant,  since  the  exponent  1.3324  exceeds  unity.  Equation 
(176)  does  not  have  this  problem  since  the  average  residual  mass  is 
always  proportional  to  the  striking  mass. 

If  we  change  one  parameter,  namely,  8  =45u,  then  we  have  the  data 

plotted  in  Figure  10.  If  we  keep  F=.l  so  K=.l/. 707=, 1414  instead  of  .1, 
then  Ng  is  multiplied  by  1.414  in  Table  IV  and  we  have  Table  V. 


TABLE  V 


Vo (km/s) 

_<  0.5 

.75 

1.00 

1.25 

1.50 

1.75 

2.00 

oo 

0 

.629 

1.512 

2.644 

4.030 

5.667 

7.556 

00 

V1 

1 

1.629 

2.512 

3.644 

5.030 

6.667 

8.556 

GO 

\v/mo 

1 

.614 

.398 

.274 

.199 

.150 

.117 

0 

the  mass  ratio  in  Table  V  is  given  by  the  dashed  curve  in  Figure  10 
which  is  lower  than  the  dashed  curve  of  Figure  9,  following  the  trend 
of  the  data  for  the  largest  pieces.  Equation  (179)  with  0q=45  leads  to 

m/m  =  1-.431(Vq) -9911  (181) 

which  is  plotted  as  the  solid  line  in  Figure  10  and  predicts  negative 
m1  beyond  v^=2.34  km/s. 

Since  we  have  a  Poisson  distribution  for  the  number  of  breaks,  we 

can  find  the  median  value  s=s  by  calculating  partial  sums.  For  example, 

if  we  use  ry=N,=1.512  for  v  =1.0  km/s  in  Table  V  we  find  the  first 
B  o 

partial  sum  (£= 0)  is  .220,  while  for  1  the  cumulative  sum  is  .553. 

Consequently,  the  median  number  of  breaks  is  s=l  and  the  median  number 
of  pieces  is  2.  Thus  half  the  time  we  should  find  fewer  than  2  pieces 
and  half  the  time  more  than  2  pieces  with  an  average  mass  fraction 

m  , ,/m  =.398.  If  we  count  the  number  of  data  points  above  and  below  the 
AV  o 

dashed  line  in  Figure  10  for  the  range  0.9<V<1.1  km/s,  we  find  ten  below 
and  ten  above  if  we  include  the  three  eroded  masses  with  mj/m^'vO.98.  This 
is  fortuitous  since  each  data  point  represents  the  mass  fraction  of  the 
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figure  10. 


Residual  mass  fraction 
steel  plates  at  0  =  45 


versus  striking  speed  for 
argest  piece  found  and 


tungsten  spheres  perforating 
average  piece  expected. 
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largest  piece  found,  and  we  cannot  expect  a  50/50  distribution  about 
the  average  mass.  If  we  consider  v  =1.5  km/s  in  Table  V  with  ry=N  =4.03, 

0  D 

the  partial  sums  are  .018,  .089,  .233,  .426  and  .621  so  the  median 

number  of  breaks  is  s  =  4  and  half  the  time  we  expect  to  find  more  or 
less  than  5  pieces  with  average  mass  fraction  .199.  All  six  of  the  data 
points  near  v_=l.5  km/s  in  Figure  10  lie  above  the  dashed  line.  This 
is  consistentuWith  the  fact  that  the  variance  of  the  Poisson  distribution 
is  also  equal  to  ry  so  the  distribution  spreads  out  as  ry  increases  and 
the  largest  pieces  are  all  likely  to  be  above  average  in  size.  Similarly 
the  smallest  pieces  are  more  likely  to  be  much  smaller  than  the  average, 
although  no  data  was  reported  to  confirm  or  deny  this  prediction. 

In  the  report  by  Herr  and  Grabarek22  an  effort  was  made  to  observe 
the  masses  of  all  the  projectile  pieces  as  well  as  target  debris  down  to 
a  cutoff  size  of  0.019  gm.  For  example,  on  pp  103-4  of  their  report 
they  give  particle  data  behind  steel  targets  6.35  mm  thick  struck  at 
zero  obliquity  by  tungsten  alloy  rods  with  mQ=3.89  gm  and  L/D=49.2  mm/ 

2.46  mm  =  20.  ForVQ=.776  km/s  the  projectile  was  intact  with  a 

residual  mass  of  1.782  gm,  indicating  an  erosion  loss  of  2.108  gm.  They 
also  reported  a  target  plug  of  0.991  gm  and  nine  fragments  with  average 
mass  of  0.075  gm  totalling  0.673  gm.  If  these  were  all  projectile 
fragments,  then  (2. 108-0. 673) =1.435  gm  was  either  out  of  the  field  of 
view  or  distributed  as  fragments  smaller  than  the  cutoff  mass.  At  vq  = 
.869  the  projectile  broke  into  three  pieces  of  .311  gm,  .739  gm 
and  1.238  gm,  totalling  2.288  gm  with  average  mass  of  0.76  gm,  indica¬ 
ting  a  possible  erosion  loss  of  1.602  gm,  less  than  the  2.108  gm 
erosion  loss  at  the  lower  as  expected.  They  also  reported  a  target 

plug  of  0.765  gm  and  15  fragments  with  average  mass  0.062  gm,  totalling 

0.925  gm.  If  these  were  all  projectile  fragments,  then  (1.602-0.925)= 

0.677  gm  was  not  measured.  At  v  =  1.384,  the  projectile  broke 

o 

into  four  pieces  of  .285  gm.  .408  gnu  .460  gm  and  2.313  gm  with  average 
mass  0.867gm,  totalling  3.466  gm,  indicating  an  erosion  loss  of  only 
0.424  gm.  They  also  reported  a  target  spall  fragment  of  0.143  gm 
and  20  other  fragments  with  average  mass  0.0476  gm  and  maximum  mass 
0.058  gm  totalling  0.953  gm.  Obviously  not  all  of  these  fragments  could 
be  projectile  fragments  since  their  total  mass  exceeds  the  erosion  loss 
of  the  projectile.  At  least  half  of  them  and  probably  most  of  them 
are  target  fragments,  which  may  also  be  true  for  the  lower  speeds.  This 
is  the  type  of  data  we  need  both  for  estimating  vulnerability/ lethality 
and  for  judging  the  validity  of  any  model  we  propose.  Of  course,  for 
both  purposes  we  need  a  great  deal  more  data. 

This  case  resembles  the  previous  case  in  that  tungsten  projectiles 
strike  steel  plate  targets  at  zero  obliquity.  However,  the  target 
here  is  thicker  by  a  factor  of  2  to  4  while  the  projectile  is  a  long, 
thin  rod  rather  than  a  sphere.  Let  us  assume  that  the  dominant 
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difference  is  the  projectile  shape.  Since  such  a  rod  striking  at  a 
slight  yaw  of  2  to  3  degrees  is  more  likely  to  fracture  than  a  sphere, 
let  us  increase  the  value  of  F  or  K  in  Equation  (178)  by  one  third  so 

that  N  =1.333  times  larger  or  N=i.9(v  -.25).  This  gives  us  the  values 

D  DO 

in  Table  VI. 


Table  VI 


V 

0 

nb 

V1 

m  /m' 
AV;  o 

m^(obs) 

mAV(calc) 

mAV(°bs) 

s  s  +1  Pieces (obs) 

.776 

.669 

1.669 

.599 

1.782 

- 

- 

0  1  1 

.869 

.960 

1.960 

.510 

2.288 

1.167 

.763 

12  3 

1.384 

3.163 

4.163 

.240 

3.466 

.832 

.867 

3  4  4 

Here  is  taken  to  be  the  sum  of  the  masses  of  the  projectile  pieces 

observed  after  perforation  and  so  does  not  include  the  mass  lost  by 
erosion.  Presumably  the  eroded  mass  survives  as  particles  too  small  to 
measure  and  breakup  of  the  eroded  mass  occurs  as  the  projectile  exits 

the  target^.  So  m^  is  taken  to  be  the  mass  to  be  broken  in  Equation  (176). 
By  using  ND=.669  for  v  =.776  we  calculate  S(0;.669)=  .512  so 

D  O 

s=0  and  the  median  number  of  pieces  predicted  is  one  in  agreement  with 
the  observation  of  an  intact  projectile.  It  is  meaningless  to  speak  of 
an  average  mass  of  the  projectile  pieces  in  this  case.  By  using  TTB= . 960 

for  VQ~- 869  we  find  S(0; .96)=. 383  and  S(l; .96)=. 751  so  s=l  and  the 
median  number  of  pieces  is  two.  If  the  two  smaller  pieces  observed, 
namely  .739  gm  and  .311  gm,  had  been  one  piece  of  mass  1.050  gm,  then 
the  observed  number  of  pieces  would  have  been  two  (as  predicted)  while 
the  observed  average  mass  would  have  been  (1.050+1. 238)/2=l. 144  gm 
which  is  close  to  the  predicted  value  of  1.167  gm.  Such  deviations  are 
expected  of  course  because  of  the  small  sample  size  and  the  large 
variance.  Surprisingly,  both  the  number  of  pieces  and  the  average  mass 
observed  at  v^=1.384  km/s  are  quite  close  to  the  values  predicted. 

This  is  encouraging.  However,  much  more  information  is  required  before 
we  can  say  that  this  model  has  been  validated  by  comparison  with  experi¬ 
ment.  This  too  will  be  the  subject  of  future  work. 
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IX.  CONCLUSION 


The  advantages  and  limitations  of  a  number  of  previous  descriptions 
of  target  penetration  by  a  projectile  have  been  reviewed,  and  new 
approaches  have  been  proposed  which  incorporate  these  advantages  and 
overcome  many  of  these  limitations.  A  deterministic  particle  model  for 
an  eroding  penetrator  has  been  proposed  and  the  simplest  versions  of  this 
model  have  been  illustrated  by  comparisons  with  simple  experiments. 

The  model  seems  capable  of  describing  three  dimensional  trajectories  in 
a  target  consisting  of  various  density  components  distributed  in  an 
arbitrary  manner.  Although  the  illustrations  given  here  have  been 
for  single  plate  or  semi-infinite  targets,  extensions  to  other  targets 
of  interest  such  as  multiple  (possibly  spaced)  plates,  convex  or  concave 
cubes,  etc.,  are  obvious.  In  addition,  a  stochastic  model  has  been 
proposed  for  projectile  breakup  and  illustrated  both  for  chunky  fragment 
projectiles  and  long  rod  penetrators.  Thus,  the  principle  features  of 
penetration  in  the  ordnance  range  can  be  described  in  a  theoretically 
sound  yet  simple  manner,  readily  adapted  to  practical  calculations  and 
suitable  for  obtaining  a  broader  fundamental  understanding  of  penetration. 

Much  work  still  remains  to  be  done.  In  order  to  consolidate  the 
gains  already  made  here,  many  more  comparisons  with  experiment  are  in 
order.  From  a  theoretical  point  of  view  it  would  be  advantageous  to 
find  closed- form  expressions  for  s(t)  and  m(t)  using  a  force  which  depends 
on  the  first  and  second  powers  of  the  speed  as  well  as  on  at  least  the 
first  power  of  the  depth  of  penetration.  In  view  of  the  long  history 
of  this  unsolved  problem,  this  may  not  be  easy  to  do  and  approximate 
solutions  may  have  to  be  employed  in  order  to  describe  hypervelocity 
penetrations  by  this  technique.  This  reminds  us  also  of  the  desirability 
of  obtaining  a  simple,  closed- form  solution  for  both  projectile  and 
shaped  charge  jet  penetration  in  a  unified  theory  as  well  as  a  unified 
model  of  projectile  breakup,  jet  particulation  and  high  explosive 
fragmentation.  Another  area  where  much  work  remains  to  be  done  involves 
the  mass  and  velocity  distributions  which  describe  the  target  debris . 
Little  experimental  information  is  available  to  guide  or  check  the 
model  maker  much  less  the  vulnerability  analyst  who  might  use  such  a 
model.  However,  now  that  a  simple  method  is  available  to  describe 
ricochet,  this  long-neglected  feature  can  be  incorporated  into  vulner¬ 
ability  codes  after  suitable  experimental  checks  have  been  made. 
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APPENDIX 


A  COMPUTER  CODE 


10 
20 
30 
40 
50 
60 
70 
80 
90 
1  00 


REM  THIS  PROGRAM  FOR  THE  HP  9830A  CALCULATOR  COMPUTES  AND  PLOTS  EXIT 
REM  SPEEDS  AND  ANGLES  VERSUS  STRIKING  SPEED  FOR  A  PENETRATOR  OF  GIVEN 
REM  MASS  M < GM >  STRIKING  A  PLATE  TARGET  OF  GIVEN  THICKNESS  T<MM>  AT  A 
REM  GIVEN  OBLIQUITY  THETA (DEG).  THE  Z  COMPONENTS  OF  THE  RICOCHET  AND 
REM  PERFORATION  LIMIT  VELOCITIES  <KM/S>  ARE  REQUESTED  AND  USED  TO  COMPUTE 
REM  HZ  AND  CZ  FROM  DZ  INSTEAD  OF  USING  THE  VALUES  GIVEhL  IF  THESE  LIMITS 
REM  ARE  UNKNOWN? SET  THEM  TO  ZERO  AND  THE  GIVEN  AZ  AND  CZ  WILL  BE  USED. 

REM  ANY  VALUES  ARE  ALLOWED  FOR  AX?  BX? CX? DX? AZ?  BZ  AND  CZ  BUT  ONLY 
REM  DZ<M*(BZ/2)t2  IS  ACCEPTED  AS  WRITTEN  TO  KEEP  BETAZ  REAL. 

D I M  V[  1 00  ]  j  CC  1 00  ]  j  NC  1 00  ]  ?  ZC  1003?  SC  1 00  ]  j  PC  28  3  ?  QC  20  3  ?  RC  20  3  ?  AC  b  1  3  ?  GL"  1 OO  3 


110 
1 20 
130 
140 
150 
160 
170 
180 
190 
200 
2 1 0 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
3  SO 
390 
400 
410 
420 
430 
440. 
450 
460 
470 
480 
490 
500 


R  =  0 

D I SP  " M  >  THETR)  T >  V0ZRL  j  VQZF'L "  5 


INPUT  M>T7>T8jV8»V9 

PR  I  NT  " MfiSS= " 5  M ! "  OBL I QU I TY= "  5  T7 5 "  TH I CKHE 


rUSF’  "  V M I N  s  V M ft X  >  V STEP"! 

INPUT  P1»P2»P3 

PRINT  "VMIN="5P1 5 "VMRX=")P2) 
BISP  "RZ? BZ> CZ) BZ' 5 


f'STEP="  5  P3 


INPUT  C5) C6 j  C7) C8 


5T8J 


02=0. 5*C6/M 

IF  <fi2't'2-C8/M><0  THEN  ISO 

B  2  =  8  U  R  fl  2  T  8  -  L  8  M 

G2=-R2+B2 

H2=-R2-B2 

IF  V8=0  THEN  290 

IF  V9=0  THEN  290 

j”;  ^ (  W  C{  /  (  V  Q  —  'y*  3  [j  j  •+•  0  Q 

C  7  =  ( -  T  8  /  ( V  9  -  V  8  >  >  *  C  8 
PRINT  "RZ="5C5J"  BZ=" »  CS5 " 
BISP  "RX) BN) CX) BX" ) 

INPUT  C1)C2)C3)C4 
R 1  =  0 . 5*C2.-'M 
W  1=0 
Bl  =  l 

IF  (Rlt2-C4/M><0  THEN  400 
B 1  =  S  Q  R  ( R 1 T  2  -  C  4  •••"  M  > 

Gl=-fil+Bl 
HI =-fll -B1 
GOTO  410 

N 1  =  S  Q  R  C  4  M  ~  fl  1 T  2 ) 

PRINT  "flX=" SCI? "  BX=" 5 C2J " 
T6=T7*PI/180 


BZ= 


;  cs 


cx=" ) C3 ) 


BX= " ; C4 


PRINT 

BISP  "THANK  YOU.  I’LL  WORK  ON  IT" 
WRIT  1000 

F'R  I  NT  TAB 2 "  VO "  TRB9 "  T I  ME "  TRB23 "  VX "  T RB 


K=0 

FOR  V0=P 1  TO  P2  STEP  P3 
K=K+ 1 
SC  K  ]  =  V0 


" V0ZRL 


■'0ZPL=  ‘ 
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519  IF  V0  =  0  THEN  10=10 

520  S 1  =V0*S  I H  <:  T 6 ) 

5  3  0  S  2 = V  0  *  C  O  S  <  T  6  > 

5  4  0  D 1  =  "  C 1  *  T  8  *  T  Fl  N  <  T  6 ')  +  f :  3  * i  j  .■■■  r  4 
550  D 2 = - <  C5*T 8 + C 7 * s p )  / r p 

560  IF  S2CV8  THEN  970 
570  IF  S2< V9  THEN  930 
580  REM***F'EEFORflT  1 fiN*** 

590  IF  R=  1  THEN  740 

609  1=1 

610  T0=0 

620  FOR  T=T0  TO  500  STEP  I 
630  IF  B2=0  THEN  660 
640  J=FNfl(T) 

650  60 TO  670 
660  J  =  FNE ( T > 

670  IF  J > T 8  THEN  690 
680  NEXT  T 

6  9  0  i f  < j - 0.0i > <  r a  the  n  7  3 0 
700  T0=T~I 

710  1=1/10 
720  60 TO  620 
730  6 L K ]  =  T 
740  T =GC K  ] 

750  IF  82=0  THEN  780 
760  V2=FNB(T) 

770  GOTO  790 
780  V2=FNF ( T ) 

790  IF  81=0  THEN  850 
800  IF  N1>0  THEN  830 
810  V 1 =FND ( T ) 

820  GOTO  860 
830  V 1 =FNL ( T  > 

840  GOTO  860 
850  V 1 =FNH ( T  > 

8  6  0  V  C  K  ]  =  S  Q  R  (  V 1 1 2  +  V  2 1 3  i 
8  7  0  C  [  K  ] = H  T  N  <  V 1  /  V  2  >  *  ( 1 8  0  /  p  I  > 

880  NR  I TE  ( 1 5  > 890 ) V 0 ?  T  5  V 1 ?  V2 

8  9  0  F  0  R  M  fl  T  F  5 . 2  5  F  9 . 2  -  F 1 5 . 3 E  1  6  3 
900  GOTO  1370 

010  IF  II 2/0  THEN  970 
9 2 0  REM*** E M 8 E D M E N T  *  *  * 

930  VI K  ]=0 
940  Cl k  ]=0 
950  GOTO  1370 

9  6  0  p  e  f'i  *  +  •¥•  r  j  f;  f, :  h  e  t  +  *  + 

970  1=1 

980  IF  R= 1  THEN  1220 
990  IF  K >2  THEN  1100 
1000  IF  82=0  THEN  1030 

1 0 10  T 1  =  L  U  6  '•  'i.  H  2  /  G  2  >  *  ( S  2'  +  G  2  *  D  2  •<  (  8  2  +  H  *  D 

1020  GOTO  1040  . 

1 S  3  0  T 1  =  <  S  2  /  fl  2 )  /  ( S  2  -  fl  ?  *  D  3 

1040  IF  T 1 >0  THEN  1100 
1050  VC K 3=0 
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CC K 3= 1 80-T7 
GC  K  3=0 

WRITE  Cl  5  ? 8 9 0 ) V 0 ?  V U  ? V U ?  V 0 
GOTO  1370 

FOR  T=T1  TO  50O  STEP  I 
IF  B2=0  THEN  1140 
J=FNA<T> 

GOTO  1150 
J=FHECT) 

IF  J < 0  THEN  1170 

NEXT  T  _ 

IF  C  J+0 . 0 1 >  >0  THEN  1210 

T 1 =T- I 
1=1/10 
GOTO  1100 
GC  K  3=T  • 

T=GC  K.  ] 

IF  B2=0  THEN  1260 
V2=FNBCT> 

GOTO  1270 
V2=FNF (T) 

IF  B 1=0  THEN  1330 
IF  W1>0  THEN  1310 
V1=FNDCT> 

GOTO  1340 
V 1 =FNL C  T  > 

GOTO  1340 
V1=FNHCT> 

VC  K  ] = S Q R  C  V 1 1 ' 2 + V 2 + 2 ) 

C  C  K  3 = A  T  N  C  V 1  /  V  2 )  *  <  1 8  0  P  l  >  +  3 *  y 
N  RITE  < 1 5 j  8 9 0 >  V 0  ?  T  ? V 1 ?  V 2 

NEXT  V0  _  -Drcri"TfiP — ci “ F I T  SPEED ” TRB40 "EX I T  ANGLE  (DEGREE 
PRINT  "IMP'Hl'T  BREED  TH&iU  tn ll  •-■rc.u. 

WRITE  <15?  1410 ) SC  L.  3?  VC  L  3?  CC  L  3 
F ORMRT  F 1 0 . 2  ?  E20 . 3 1  F20 . 

REm'  (REMOVE  REM  TO  RECYCLE!:'  GO  TO  29^ 
msp  "PLOT  THESE  VALUES''  C  tEo=1  ?  NU-O-  j 
INPUT  P5 

IF  P5=0  THEN  2710 
REM  PLOT  SPEEDS  AND  ANGLES 
DlfiP  "DATA  FILE  NUMBER?  <U=NUNE,'  » 

INPUT  N1 

IF  N1=0  THEN  1600 
LOAD  DATA  N1?A 
11=2 

FOR  1  =  1  TO  AC  1 3 
PC  I 3=AC  1 1  3 
QC I  3= AC  1 1  +  1 3 
RC  I  3= AC  11+23 
11=11+3 
NEXT  I 

GOTO  1680  . 

1600  DISP  "#  DATA  POINTS? FILE  #  ? 


1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1 190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1 3 1 0 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1 400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
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1 6  i  0 
1 620 
1630 
1640 
1  650 
1 660 
1 670 
1680 
1 6  y  0 
1  700 
1710 
1  720 
1  730 
1740 
1750 
1  760 


INPUT  AC1],N1 

FUR  1=2  TO  3*AC  1  ]+i  STEF'  8 
Dlyp  " V0, v j THETfl" ;  " 

HE^TTrHC  1  3’  FlC  1  +  1  AC  1+2  I 

STORE  DATA  N 1  <  A 
GOTO  1520 

P4=<P2-pi  )...■•  10 

SCALE  P1-P4?  P2  +  0. 5*P4, -P4,  po+fi  <=-. 

PLOT  P 1  j  R2)  1 

CP LOT  10,0 

!P_W1 >0  THEN  1758 

LABEL  ( 1 7 6 0 > ' ' E X IT  S P E E H  r' y m c  •, 

GOTO  1770 

LABEL  (17 6 0 ) " E X IT  sPFF n  i- M  ■■  ■■ 

format  2F8.3  “  M' U 


BE  I  A? 


•MEGA; 


1  *  O  01 

•  l.  u  i  r 

CPL.OT 

l  >  r  4  ,  1 
10,  -2 

1  790 

label 

(  1  800  j 

MM=:|  , 

M,  “  j". 

CT1 !  T  i  {  t  t  s  i  i .  _ _ 

1800 

format 

4F8 . 4 

i  1  ?  i j 

DL1UU1 1 r=  9 T?n  " 

r  i 
l  - 

1 8 1 0 

1  p  O 

LABEL 

•  !~!  r~  t 

(  1  800  ) 

?  C  i  *  " 

B  i  !_a  ■  { ”■  y  --  ■ 1 

1 830 

LnbtL  1 

LABEL  ' 

(  1 800  > 1 
!.  •+•  ’)  "  “ 

l,flZ  = 

,  C5»  " 

BZ= " , C6 , "  C2=" 

H  !  .’J  “ 

?  C  7  H 

1840 

CPLOT  f 

:1  ?  - ! -j 

1850 

label  < 

.  *  )  " 

lltl 

0 

j  |  i| 

1  o'  b  0 

LABEL  t 

*  ■+■  )  "  " 

1870 

FOR  1=1 

TO  K 

STEP 

lNTck, 

8 ) 

1 8  8  0 

1  PQCi 

LABEL  ( 

JIT  M  D  M  L!  T 

1 890 ) S 

T  I  ], 

VC  I  ] 

a  1  r»r 

r  ur.rin  i 

R- 1"  1  0  ■  -.j 

1  900 

NEXT  I 

1  9 1 0 

XAXIS  0 

»  0  .  d 

1  ..  C  ~. 

1  920 

LABEL  ( 

W  *  .»  f 

-+■  j 

j.  21  r 

1  930 
1940 
1950 
1 9  60 
1970 
1  980 
1  990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2078 
2080 
2090 
2100  | 
2110  i 
2120  [ 
2130  f 
2140  F 
:150  f 


Pi°nTX-P-  T°  P2  STEP  0-2 
T  L  U  T  ,  H  ,  J 

'  CP LOT  -2,-£ 

1  LABEL  (  1 970  )  v 

1  FORMAT  F4.1 

NEXT  X 

V  A  X I S  p  1 , 0 .  2 1  0 ,  p 
LABEL  (*>  “ 

FOR  y=0  TO  P2  STEP  0.2 
FLuT  0, Vi  1 
CPLOT  -5, -0.3 
LABEL  (1970 ’5V 
NEXT  V 

PLOT  P2j 0.1,1 
CPLOl  -20,0 

EPEEE  <*>:  IMPACT  SPEED  (KM. •■"S' > " 

fur  i  =  i  tu  ac i  ] 

PLOT  PC  I  ],Q[  n,  1 
CPLOT  -0. 3, -0. 3 
LABEL  ( * ':■  "X" 

NEXT  I 

FOR  L=1  TO  K 
PLOT  SC  L  ] ,  VC  L  3 


BETAZ= :: ,  E'2 
BETflZ=",B2 

DX= " , 04 


104 


2160  NEXT  L 

2180  DISP  “CHANGE  PAPER.  THEN  CONTINUE 
2190  STOP 

2200  P4=(P2~P1  10  _  ..  . 

2210  SCALE  Pi-P4»P2+P4^»-^»*10 
2220  PLOT  PI » 206?  1 
2230  CP  LOT  10»0 

2240  IF  W1>0  THEN  22."0  ••  nrrPFF1- "• 

2299  LABEL  C2280VEXIT  ANGLE  cBEuREEo. 

2 2 6H  I j U T U  2290  _  _ .  _  .•  rirrC'CC c  'i 

o:-'7L“i  I  ABEL  (2280)  "EX II  HNcLE  DELKtt-' 

2280  FORMAT  2F8.3 
ooqfi  PLOT  F'l  j  200?  1 

P?S  Sel  ««5>V.ir  OBLIQUITY-’. T7.“  T 

Hit  LABEL1  cSIetlAX*;. Cl.  :  |X!:;C2.: 

8 840  LABEL  <2320>-H2-  -LS-  6-*  ,LB’ 

2250  LABEL  (*)"  " 

2360  CF'LOT  25?  -2 
2370  LABEL  (*)" 

2380  LABEL  (*)"  "  ,JT  ,o  , 

90  FOR  1  =  1  TO  K  bTEF  INT -K.- ..... 

2400  LABEL  (24 10)  SCI  3> >-C  I  3 
2410  FORMAT  2F10.3 
2420  NEXT  I_  . 

2430  XAXIS  0j0.2jF1jF* 

2440  LABEL  <*)  p  fl 

2450  FOR  X=F'l  3  0  F c.  olLF  0 ■  - 

246  0  F'  L  U I  '  0  ’  I  _ 

2470  CF'LOT  -2? -2 
2480  LABEL  (2490)X 
2490  FORMAT  F4.1 


BETAX=‘ 

•jBIi  " 

BETHi-  » B 

OMEGAX=' 

' »  W 1  ?  " 

BETAZ  =  " j  B 

?  T8 


J  C4 

7=" , C8 


V0 


THETA- 


500 

NEXT  X 

510 

VAX  IS  F'l  i  2 

0  J 

0?  2 

00 

520 

LABEL  (*) 

1530 

FOR  Y=0  To 

•“( 

:00 

STEP 

:540 

PLOT  0>Yj1 

:550 

CF'LOT  -61- 

0. 

:560 

LABEL  (25  i'1 

0 : 

•  Y 

5:570 

FORMAT  F5. 

0 

2580  NEXT  Y 
2590  PLOT  P 2  ?  1 0 »  1 

2 b 0 0  C PLUl  —  0 '0  _ _  r. [T CTl  I-'' M  / .1 

2810  LABEL  (*) " IMPHOT  oFEEI 
2820  FOR  1  =  1  TO  AC  1 3 
2630  PLOT  PC  I 3j  RC I  3 » 1 
2640  CF'LOT  -0.3>~0.3 
2650  LABEL  (*) "X" 

2660  NEXT  I 
2670  FOR  L= 1  TO  K 
2 6 00  PLOT  SC  L  3 ?  CC  L  3 
2690  NEXT  L 
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-  (?  0  0 

PEN 

i'7  1 0 

STOP 

2720 

END 

2730 

DEF 

FNfiCT 

— S  2  +  D  2  *  H  2  )  *  E X  F'  (.  G  2  *  T  —  i.  S  2  +  D  2  * G  2 

>*ENP'- H2*T':< 

1  (.  G 2 -H2  )  +  i.i E 

2  7  40 

DEF 

FHBCT 

=  <  (  S  2  +  D  2  *  H  2  >  *  G  2  *  E  X  P  <  G  £  *  T  >  -  ( S  2  +  D  2 

*G£  )  *H2*EXF 

042*1)  •)/  CG2-H2 

2750 

DEF 

FNCCT 

= < < S 1 +D 1 *H 1 >*EXP < G1*T>~ < S 1 +D 1 *G 1 

> *EXP i  H 1  *T  > 

)■  ■  •:  G  [--Hi  )  +  D 1 

2760 

DEF 

FNIKT 

=  ( (.  S 1  +  D 1  * H 1 )  * G 1  * E X P ( G 1  * T )  -  <  S 1+D1 

*  g  1  *  Pi  l  it  t  X  F- 

H  1  *  T  ) >  •••■'  <  G  1  -  H 1 

2770 

DEF 

FNECT 

= E X P<-fi2*T>* < ( S2-R2*D2)*T-B2  > + D  2 

2780 

DEF 

FNF  <  T 

=EXP ( — R2*T >  * ( < S2-fl2*D2 ) * ( 1 -R2+T ) 

+H2*D2  > 

2790 

DEF 

FNGCT 

=EXP  ( ( ( S 1  -H 1  *B  1 )  *T~Ii  i )  +  D 1 

2800 

DEF 

FHH(T 

'  =EXP <-ft 1*T) *  < < S 1 -ft  1 *D 1 ) * ( i -H 1  *T ) 

*  h  i  *  D 1  ) 

28 1 0 

DEF 

FNKCT 

>  =  E X P «.  - H 1  * T )  *  ( (. S 1  — R 1  *D  1  >* S I N ( N 1  *  T 

/  14 1  -  D 1  *  L  U  b 

041*1)  '  +  Ui 

2820 

DEF 

FNLcT 

1  =  E  p  V  -  H  1  *T  b  1  *  t  IJ  s »,  M  1  *1  )  H  i  * b  i 

-D 1 *C4  "M )*S 

IN  041*1)  , 4 41 

28  30 

ST  OF 

2840 

END 
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